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Final  Report  Summary 

This  reports  on  a  highly  successful  two  year  program  whose  objective  was 
to  investigate  the  structural  and  optoelectronic  properties  of  a  novel  engineered 
materials  system:  vacuum  deposited  semiconducting  organic  thin  films.  The 
program  was  terminated  early  due  to  the  move  of  the  PI  from  DSC  to  Princeton. 
Hence,  this  final  report  concerns  only  the  research  done  at  DSC  over  the 
shortened  time  frame  of  the  project. 

This  research  has  led  to  the  following  breakthrough  results: 

1 .  Demonstration  of  the  first  organic  thin  film  multiple  quantum  wells  using  the 
novel  growth  technique  of  organic  molecular  beam  deposition  (OMBD).  This 
growth  technology  developed  in  our  lab  has  set  a  standard  now  being 
investigated  worldwide.  Approximately  2-4  papers  each  month  are  appearing  in 
Applied  Physics  Letters,  with  additional  papers  appearing  in  such  journals  as 
Physical  Review  Letters,  Japanese  Journal  of  Applied  Physics,  etc  Virtually  all 
groups  have  adopted  our  terminology  of  OMBD  and  "quasi-epitaxy"  (QE).  This 
latter  term  refers  to  the  growth  of  lattice-mismatched  structures  which  are  bonded 
by  the  flexible  van  der  Waals  forces.  In  addition,  a  large  percentage  of  papers 
currently  being  published  by  other  groups  refer  to  two  papers  generated  by  our 
group,  one  of  them  done  under  this  current  contract: 

"Evidence  for  Exciton  Confinement  in  Crystalline  Organic  Multiple  Quantum 
Wells",  F.  F.  So  and  S.  R.  Forrest,  Phys.  Rev.  Lett.,  66.  2649  (1991). 


2.  Discovery  and  analysis  of  exciton  confinement  in  the  first  organic  multiple 
quantum  well  structures.  This  work  has  numerous  exciting  potential  applications 
in  optical  devices  such  as  modulators,  light  emitters  and  detectors.  We 
measured,  for  the  first  time,  the  exciton  radius  and  recombination  lifetime  in 
organic  MOW  structures,  and  found  that  their  energy  was  blue  shifted,  and  their 
radiative  lifetime  reduced  due  to  quantum  effects. 

3.  Demonstration  of  monolayer  control  of  organic  thin  film  growth.  Here,  direct 
observation  via  scanning  tunnelling  microscopy  of  molecular  layers,  and  their 
relation  to  the  substrate  was  made.  The  results  lead  to  the  generation  of  a 
simple  calculational  model  which  is  currently  being  used  to  predict  the  conditions 
and  materials  which  lend  themselves  best  to  quasi-epitaxial  growth. 
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In  addition  to  these  discoveries,  we  have  fabricated  several  devices  which  take 
advantage  of  the  unique  optical  and  electronic  properties  of  the  organic  thin  films. 
Much  of  this  device  work  was  done  with  support  and  collaboration  from  Rome  Air 
Development  Center  (Hanscom)  which  has  closely  followed  the  research,  and  in 
some  cases  has  transferred  some  of  our  technology  to  their  own  labs.  As  an 
example,  RADC  issued  an  SBIR  contract  solicitation  for  a  multipurpose  ultrahigh 
vacuum  system  similar  to  our  OMBD  chamber  as  part  of  their  desire  to 
incorporate  OMBD  technology  for  their  own  use.  This  collaboration  with  RADC 
continues  even  at  Princeton  University  due  to  the  importance  of  moving  this 
strategic  technology  into  other  domestic  laboratories. 
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2.  Introduction 

Recent  work  in  many  laboratories  worldwide  [1-7]  has  demonstrated  a 
new  class  of  engineered  materials  with  application  to  photonic  devices.  These 
materials  are  variously  known  as  "van  der  Waals  solids",  "layered  materials"  or 
"quasi-epi^axial  materials".  The  property  which  these  materials  have  in 
common  is  that  the  cohesive  force  which  bonds  the  various  material  layers 
together  is  the  relatively  weak  van  der  Waals  force.  This  results  in  an  ability  to 
layer  materials  which  are  highly  lattice-mismatched  without  inducing  defects  in 
the  crystal  structure.  One  particularly  interesting  class  of  van  der  Waals  solids 
with  applications  to  a  very  broad  range  of  optical  and  electronic  devices  is 
based  on  multilayered  structures  consisting  of  crystalline  organic 
semiconductors  [1].  These  materials  have  been  found  to  have  excellent 
electrical  and  optical  characteristics  when  grown  by  the  ultrahigh  vacuum 
process  of  organic  molecular  beam  deposition  (OMBD)  onto  other  organic 
semiconductors,  or  even  onto  the  surfaces  of  inorganic  semiconductors  such  as 
Si,  GaAs  or  InP  [8,9].  For  example,  high  bandwidth  photodetectors  consisting  of 
layers  of  perylene-based  compounds  deposited  on  Si  substrates  have  been 
demonstrated  [8],  as  have  field  effect  transistors  based  on  organic/lnP 
heterojunctions  [10].  Furthermore,  waveguides  and  other  optical  devices  have 
been  demonstrated  using  these  quasi-epitaxially  grown  layers  [11]. 

In  this  work,  we  discuss  the  structural,  dielectric,  conductive,  and  optical 
properties  of  this  new  and  exciting  class  of  engineered  materials.  In  particular, 
we  have  demonstrated  the  growth  of  multiple  quantum  well  (MOW)  structures 
based  solely  on  organic  semiconductors.  These  structures  have  demonstrated 
quantum  confinement  of  excitons  [1,2],  opening  the  door  to  a  wide  range  of 
deyices  such  as  optical  modulators  and  switches  which  can  be  tailored  to  match 
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a  particular  window  in  the  visible  and  infrared  spectral  regions.  Furthermore, 
we  have  observed  giant  anisotropies  [3]  in  the  dielectric  and  conductive 
properties  of  these  epitaxial  films  which  are  a  direct  result  of  their  crystal 
structure. 

Due  to  the  ease  with  which  these  materials  can  be  grown  into  crystalline 
structures  without  regard  to  lattice  match  constraints,  as  well  as  due  to  the  wide 
variety  of  materials  which  can  be  chosen  for  a  given  application,  the  process  of 
OMBD  of  van  der  Waals  solids  opens  the  door  to  an  unprecedented  range  of 
new  challenges  and  opportunities  to  the  materials  scientist,  both  in  the 
understanding  of  materials  as  well  as  in  the  ability  to  fabricate  new  devices  for 
use  in  advanced  optoelectronic  systems. 

The  purpose  of  this  paper  is  to  review  our  recent  results  in  the  growth 
and  characterization  of  quasi-epitaxial  films  of  two  different  organic  molecules, 
and  to  discuss  the  unique  optical  and  electrical  properties  of  the  films  and 
multilayer  structures  consisting  of  two  different  materials  interleaved  in  very  thin 
(  >  10  A)  layers. 

The  paper  is  organized  as  follows:  In  the  following  section  we  discuss 
the  process  of  quasi-epitaxy.  Both  the  techniques  employed  and  a  physical 
model  for  the  growth  process  will  be  presented.  In  Sec.  3,  we  discuss  the 
electrical  properties  as  well  as  the  conductivity  mechanisms  dominant  in 
heterojunctions  consisting  of  two,  layered  organic  films.  In  Sec.  4,  the  optical 
properties  of  organic  multiple  quantum  well  (MQW)  structures  are  considered  in 
detail,  and  in  Sec.  5  we  discuss  the  future  prospects  for  quasi-epitaxially  grown 
organic  semiconductors. 

3.  Structure  of  quasi-epitaxially  grown  films 

Quasi-epitaxy  (q-e)  is  the  process  whereby  a  film  which  is  strongly  lattice- 
mismatched  to  a  substrate  can  be  grown  in  a  relatively  defect-free,  ordered 
layer.  This  seemingly  contradictory  state  of  affairs  is  a  direct  result  of  the  weak 
nature  of  the  van  der  Waals  (vdW)  interaction  which  is  the  cohesive  force 
binding  the  constituent  materials.  Here,  it  is  assumed  that  the  internal  binding 
energy  of  a  given  material  is  significantly  greater  than  the  binding  energy  of  the 
dissimilar  materials  used  in  the  layered  structures.  Although  there  is  strain 
energy  associated  with  the  mismatch  between  adjoining  layers,  this  energy  is 
insufficient  to  induce  defects  within  a  particular  layer. 
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In  our  experiments,  we  have  used  the  archetype  compound,  3,  4,  9.  10 
perylenetetracarboxylic  dianhydride  (PTCDA)  to  study  the  mechanisms  of  q-e 
growth.  The  unit  cell  of  PTCDA,  or  C24O6H8,  is  shown  in  Fig.  1.  The  molecule 
forms  tightly  bound  stacks,  with  an  intermolecular  spacing  of  only  3.21  A. 
Furthermore,  adjacent  stacks  are  oriented  (in  the  bulk  crystal)  with  the  long 
molecular  axes  at  90°  with  respect  to  each  other.  There  are  several  physical 
implications  which  arise  from  the  obvious  asymmetries  of  this  structure.  For 
example,  the  close  stacking  distance  gives  rise  to  extensive  j:-orbital  overlap 
between  adjacent  molecules  within  a  stack,  whereas  almost  no  overlap  of  the 
molecular  orbitals  exist  between  the  stacks.  This  results  in  large  conductive 
and  dielectric  asymmetries  along  various  crystalline  axes,  as  will  be  discussed 
below. 

The  first  demonstration  of  quasi-epitaxial  growth  of  an  organic 
semiconductor  deposited  onto  a  glass  substrate  was  by  Forrest,  et  al.  in  1984 
[12].  In  that  work,  it  was  found  that  when  PTCDA  was  deposited  at  a  rate 
exceeding  50  A/s  onto  a  glass  substrate  held  at  room  temperature,  single 
crystalline  thin  films  could  be  grown  which  extended  several  centimeters  across 
the  substrate  surface.  Following  that  initial  work,  Debe,  et  al.  [13]  demonstrated 
that  other  perylene  derivatives  similar  to  PTCDA  could  attain  single  crystallinity 
when  deposited  more  slowly  onto  substrates  held  well  below  room  temperature. 
More  recently,  we  have  found  that  PTCDA,  and  an  analogous  compound, 
3, 4, 7, 8  naphthalenetetracarboxylic  dianhydride  (NTCDA),  or  C14O6H8,  can 
both  be  grown  quasi-epitaxially  at  a  rate  of  1  -  3  A/s  when  deposited  onto  glass 
substrates  held  at  approximately  90K  [1].  Low  substrate  temperature  during 
growth  is  required  since  the  thermal  energy  of  the  molecules  on  the  surface 
needs  to  be  brought  well  below  the  intermolecular  bond  energy  for  ordered 
growth  to  be  established. 

Multilayer  stacks  of  alternating  layers  of  PTCDA  and  NTCDA  also  have 
recently  been  found  to  form  ordered  crystalline  structures.  This  order  is 
achieved  even  though  the  crystal  structures  of  PTCDA  and  NTCDA  are  highly 
mismatched.  In  Fig.  2  we  show  a  unit  cell  of  NTCDA  (c.f.  Fig.  1).  Here,  it  is 
seen  that  rather  than  forming  planar  stacks  as  in  the  case  of  PTCDA,  NTCDA 
molecules  stack  in  a  herringbone  habit,  with  an  intermolecular  spacing  of  3.5  A. 
In  the  herringbone  structure,  adjacent  stacking  axes  are  positioned  at  90°  with 
respect  to  each  other,  once  again  leading  to  anisotropies  in  the  conductive  and 
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dielectric  properties  of  these  crystals.  However,  the  asymmetries  are  much  less 
pronounced  than  those  found  in  PTGDA. 

In  growing  the  multilayer  stacks  of  PTCDA  and  NTCDA  onto  glass 
substrates,  pre-purified  powder  samples  of  the  materials  were  loaded  into 
separate  sublimation  cells  in  the  OMBD  growth  chamber  which  had  a  base 
pressure  of  lO-''0  Torr.  The  materials  were  then  alternately  deposited  by 
sequentially  heating  and  shuttering  each  cell  to  a  temperature  high  enough  to 
sublime  the  molecules  at  a  rate  of  1  -  3  A/s  (approximately  400  °C  for  PTCDA 
and  250  °C  for  NTCDA)  as  determined  using  a  crystal  thickness  monitor,  while 
the  temperature  was  kept  well  below  the  dissociation  temperature  for  the 
particular  molecule  being  deposited.  Typically,  from  1  -  20  layers  of  each 
material  with  individual  layer  thicknesses  ranging  from  10  A  to  1000  A  were 
grown  in  this  fashion. 

Evidence  for  crystalline  order  has  been  obtained  using  both  x-ray 
diffraction  as  well  as  optical  birefringence  measurements  [1].  In  the  case  of  x- 
ray  diffraction,  it  was  found  that  both  PTCDA  and  NTCDA  grew  into  their  normal 
crystal  stacking  habits  even  though  the  layers  alternated  between  the 
mismatched  crystals  of  PTCDA  and  NTCDA.  In  a  separate  series  of 
birefringence  experiments,  it  was  found  that  the  optic  axis  of  PTCDA  in  such 
multilayer  structures  was  oriented  in  a  single  direction  across  broad  wafer 
areas,  and  that  the  birefringence  of  each  of  the  PTCDA  layers  in  the  stacks 
added  coherently,  even  though  the  PTCDA  layers  were  separated  by 
intervening  NTCDA  layers.  The  experiment  was  accomplished  by  placing  a 
multilayer  sample  on  a  glass  substrate  between  crossed  polarizers,  and  then 
measuring  the  intensity  of  the  transmitted  light  through  the  polarizers  as  the 
sample  was  rotated  about  its  optic  axis  which  was  oriented  along  the  beam 
direction.  By  using  a  tightly  focussed  beam,  this  birefringence  measurement 
determines  the  local  orientation  of  the  optic  axis  of  PTCDA  (NTCDA  has  only 
very  weak  birefringence  and  hence  the  probe  is  not  sensitive  to  its  presence  in 
the  stacks).  This  technique,  therefore,  is  sensitive  to  the  degree  of  crystalline 
order  over  the  region  covered  by  the  beam  diameter.  The  finding  that  the  optic 
axis  of  PTCDA  is  oriented  in  the  same  direction  across  the  wafer  surface,  and 
that  the  birefringence  of  a  multilayer  stack  is  equal  in  magnitude  to  a  single 
PTCDA  film  of  equal  total  thickness,  indicates  that  the  crystalline  orientation  of 
one  PTCDA  layer  is  the  same  as  all  the  other  PTCDA  layers  in  a  multilayer 
structure.  By  inference,  the  NTCDA  layers  must  be  likewise  aligned.  This 
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ordering,  both  across  wafer  surfaces  as  well  as  between  layers  in  the  stacks  is 
evidence  for  quasi-epitaxy  in  the  multilayer  samples. 

The  process  of  q-e  can  be  quantitatively  understood  if  we  require  that  a 
monolayer  of  one  material  must  find  a  minimum  energy  configuration  with  the 
underlying  crystal  during  growth.  That  is.  when  growth  of  material  B  is  initiated 
on  the  surface  of  material  A,  the  first  few  molecules  B  align  themselves  to  A  to 
minimize  their  individual  molecular  interaction  energies.  As  the  growth  of  layer 
B  approaches  a  full  monolayer,  different  islands  of  material  B  arising  from 
nucleation  at  random  sites  across  the  surface  of  A  will  be  out  of  registry  with 
each  other,  although  they  may  all  be  approximately  aligned.  As  the  spacing 
between  islands  is  filled  in  with  additional  deposition  of  material  B,  the  islands 
must  shift  slightly  (but  in  a  rigid  manner)  to  accommodate  these  last,  interstitial 
molecules.  This  process  must  occur  without  increasing  the  overall  energy  of  the 
layer  to  the  extent  that  a  dislocation  is  generated.  Note  that  this  picture  requires 
that  the  substrate  temperature  be  high  enough  to  allow  for  some  surface 
mobility  of  the  islands  and  molecules,  yet  not  so  high  that  the  individual 
molecules  are  only  very  loosely  bound  to  the  surface.  For  this  reason, 
deposition  of  q-e  grown  materials  occurs  on  substrates  cooled  below  room 
temperature. 

One  condition  which  must  be  fulfilled  tor  the  above  scenario  to  occur  is 
that  the  binding  energy  minimum  between  molecules  A  and  B  be  very  broad  -- 
i.e.  there  is  a  range  of  positions  which  molecule  B  can  have  when  aligned  to 
molecule  A  which  does  not  significantly  change  the  total  energy  of  the  system. 
In  this  way,  small  shifts  in  position  can  be  accommodated  without  significantly 
increasing  the  total  crystal  energy.  We  have  tested  this  hypothesis  by 
calculating  the  vdW  bond  energy  between  PTCDA  (molecule  A)  and  NTCDA 
(molecule  B).  These  calculations  are  done  using  the  atom-atom  potential 
method  [14].  In  this  method,  we  calculate  (and  minimize)  the  vdW  energy 
between  each  pair  of  atoms  in  both  the  PTCDA  and  NTCDA  molecules.  The 
total  bond  energy  is  simply: 

0  =  L<J)ij  (1) 

where  0ij  is  the  atom-atom  potential  between  the  ith  and  jth  atoms  in  molecules 
A  and  B,  respectively.  This,  in  turn,  is  calculated  using  the  Buckingham 
potential: 

bij  =  -a/rij6+  pexp(-Yrij).  (2) 
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Here,  Hj  is  the  distance  between  atoms  i  and  j  in  molecule  B,  and  a,  p  ar.d  y  are 
vdW  atom-atom  potential  constants  for  the  constituent  atoms  in  the  two 
molecules  (C,  O  and  H  for  PTCDA  and  NTCDA,  for  example).  The  constants 
used  in  the  case  of  PTCDA  and  NTCDA  are  given  in  Table  1 . 

The  accuracy  of  this  method  has  been  tested  by  calculating  the  minimum 
energy  configuration  of  bulk  PTCDA.  and  comparing  the  results  to  the 
measured  crystal  structure.  Under  this  test,  the  model  works  surpri'  igly  well. 
For  example,  the  calculated  intermolecular  stacking  distance  is  3.26  A  as 
compared  with  the  actual  value  of  3.21  A.  Furthermore,  the  model  predicts  that 
two  nolecules  stack  with  their  long  axes  aligned,  but  with  their  short  (b)  axes 
offset  by  0.9  A,  as  compared  to  l.i  A  obtained  from  measurements.  Finally, 
molecules  in  adjacent  stacks  favor  a  perpendicular  to  a  parallel  relative  stack 
orientation  by  40  meV,  also  consistent  with  observation  (Fig.  1).  Thus,  we  have 
some  confidence  that  the  atom-atom  method  is  valid  for  predicting  the  stacking 
configuration  of  two  different  molecules  such  as  PTCDA  on  NTCDA. 

The  results  of  this  calculation  for  PTCDA/NTCDA  multilayers  are  shown 
in  Fig.  3.  Here,  the  diagram  shows  contours  of  constant  total  energy,  O,  as  a 
function  of  displacement  of  the  two  molecules  placed  in  a  planar  stacking 
configuration.  The  contours  arr  at  5  meV  intervals,  and  position  (0,0) 
corresponds  to  the  center  of  an  NTCDA  molecule  positioned  directly  over  the 
center  of  a  PTCDA  m.olecule.  Note  that  the  intermolecular  z-spacing  of  3.2  A 
was  determined  from  a  similar  minimization  procedure  as  that  outlined  above. 

From  the  energy  contour  plot,  we  see  that  the  energy  minimum  between 
PTCDA  and  NTCDA  is  indeed  broad,  allowing  for  nearly  a  2  A  translation  of  one 
molecule  with  respect  to  another  without  any  significant  increase  in  O.  Thus, 
we  can  conclude  from  this  first  order,  somewhat  simplified  calculation,  that  our 
"picture"  of  the  process  of  q-e  is  plausible.  That  is,  we  expect  that  once  the 
islands  are  nucleated  around  the  surface  of  the  substrate,  they  are  still  free  to 
rigidly  translate  to  accommodate  the  arrival  of  the  final  molecules  which  fill  in 
the  interstices  between  neichboring  islands.  This  shift  occurs  without  inducing 
strain  or  polycrystalline  growth. 

While  this  model  is  still  unsophisticated  in  many  ways,  it  suggests  the 
properties  which  molecules  rriust  possess  before  they  can  be  grown  into  q-e 
structures.  From  these  initial  results,  we  infer  that  the  process  is,  in  fact,  a  very 
general  property  of  vdW  solid  growth. 


4.  Organic  Heterojunctions 

It  is  well  known  that  multiple  quantum  wells  can  be  fabricated  by 
alternately  layering  two  inorganic  semiconductors  with  different  bandgaps. 
Thus,  the  question  arises:  Is  it  also  possible  to  observe  quantum  confinement  in 
fully  organic  multiple  heterojunction  structures?  This  question  can  only  be 
answered  once  we  understand  the  nature  of  a  heterojunction  between  two 
organic  semiconductors,  and  the  role  it  plays  in  governing  charge  transport 
between  the  contacting  layers.  While  there  have  been  reports  on  the  fabrication 
of  organic  heterojunctions  (HJs)  and  their  rectification  characteristics  [17,18], 
the  charge  transport  mechanisms  across  such  HJs  are  still  nc*  well  understood. 
In  this  section,  therefore,  we  describe  the  rectification  characteristics  of  organic 
heterojunctions,  and  demonstrate  that  they  can  be  understood  in  terms  of  a 
band  model  analogous  to  that  used  for  inorganic  semiconductors.  That  is,  we 
show  that  charge  transport  over  the  energy  barrier  at  the  heterointerface  is 
limited  by  the  band-like  thermionic  emission  of  charge  carriers.  From  the  results 
of  our  experiments,  we  infer  that,  indeed,  fully  organic  semiconductor  multiple 
quantum  well  (MQW)  structures  have  the  potential  for  exhibiting  many  of  the 
same  physical  properties  as  do  inorganic  semiconductor  MQWs. 

The  inset  of  Fig.  4  shows  a  fully  organic  heterojunction  device  consisting 
of  copper  phthalocyanine  (CuPc)  and  PTCDA.  In  this  structure,  indium-tin-oxide 
and  indium  are  used  as  contacts  to  CuPc  and  PTCDA,  respectively,  since  these 
materials  have  been  found  to  have  very  low  contact  resistance  to  the 
corresponding  organic  layers.  The  rectification  due  to  the  heterojunction  is 
evident  from  the  bipolar  current-voltage  (J-V)  characteristics  shown  in  Fig.  4. 
Here,  forward  bias  is  achieved  when  PTCDA  is  negative  with  respect  to  CuPc. 

We  can  understand  the  observation  of  rectification  in  the  CuPc/PTCDA 
heterojunction  in  terms  of  the  conventional  theory  for  inorganic  semiconductor 
HJs  [19].  It  should  be  noted,  however,  that  charge  transport  in  crystalline 
organic  materials  might  alternatively  be  described  by  the  hopping  model,  where 
the  transport  of  localized  carriers  between  adjacent  molecules  is  a  thermally 
activated  process.  In  this  case,  therefore,  the  carrier  mobility  should  have  a 
strong  temperature  dependence.  However,  such  a  temperature  dependence  is 
not  observed  in  organic  crystals,  therefore  suggesting  that  charge  transport  is 
inconsistent  with  the  hopping  model.  On  the  other  hand,  the  use  of  band  theory 
for  some  crystalline  organic  materials  is  justified  by  previous  work  on  hole 
mobility  measurements  on  ultra-high  purity  naphthalene  and  perylene  crystals 


[20].  In  that  work,  Warta  and  Karl  found  that  the  hole  mobility  reaches  a 
maximum  value  of  400  cm^/V-s  at  4.2  K  in  some  materials.  Furthermore,  the 
hole  mobility  was  found  to  increase  linearly  at  low  electric  field,  and  saturates  at 
high  fields.  The  high  field  carrier  velocity  saturation  is  due  to  phonon  emission, 
consistent  with  band  theory.  Further,  the  long  carrier  mean  free  path  (as 
compared  with  the  intermolecular  stacking  distance)  in  these  crystalline  organic 
materials  estimated  from  the  low  temperature  mobility  data  also  supports  the 
band  description.  Provided  that  we  can  use  this  model  to  describe  the  charge 
transport  in  PTCDA  and  CuPc,  the  transport  theory  for  inorganic  semiconductor 
heterojunctions  can  also  be  applied  (with  caution)  to  organic  heterojunctions. 

Using  the  treatment  of  Chang  for  isotype  heterojunctions  [21],  the 
saturation  current  density  for  an  HJ  is  given  by; 

Js  =  BTi/2exp[  -  qVD2/kT]  (4) 

where  B  is  a  pre-factor  which  depends  on  the  density  of  acceptors  and  the  hole 
effective  mass  in  CuPc,  kT  is  the  Boltzmann  energy  at  temperature  T,  Vd2  is  the 
diffusion  potential  on  the  depleted,  CuPc  side  of  the  heterojunction,  and  q  is  the 
electron  charge.  The  valence  band  discontinuity  energy,  AEy ,  is  then  given  by: 

AEv  =  qVoi  +  qVD2.  (5) 

In  the  above  equation,  Vqi  is  the  diffusion  potential  on  the  PTCDA  side  of  the 
heterojunction.  For  simplicity,  we  assume  that  there  is  no  significant  difference 
in  the  effective  hole  densities  of  states  and  in  the  acceptor  concentrations 
between  CuPc  and  PTCDA. 

From  Eq.  4,  it  can  be  seen  that  Js  is  thermally  activated  with  energy  qVD2- 
In  Fig.  5,  we  show  the  dependence  of  log(Js)  on  1/T  for  a  CuPc/PTCDA 
heterojunction.  Here,  Js  is  obtained  by  extrapolating  both  the  forward  and 
reverse  bias  currents  to  an  applied  voltage  of  V  =  0.  The  activation  energies  for 
both  forward  and  reverse  bias  are  the  same,  yielding  Vd2  =  0.51  ±  0.05  V.  This 
value  is  consistent  with  that  obtained  from  capacitance-voltage  (C-V) 
measurements,  where  we  find  that  Vd2  =  0.45  ±  0.15  V.  Following  the  treatment 
of  Chang,  we  obtain  Vqi  =  0.01  V,  and  hence  the  valence-band  discontinuity 
energy  for  CuPc/PTCDA  heterojunctions  is  AEy  =  Vqi  +  Vd2  =  0.52  ±  0.05  V. 

Based  on  the  electrical  characteristics  of  the  CuPc/PTCDA 
heterojunctions,  a  band  diagram  near  the  heterointerface  is  proposed  in  Fig.  6. 
Here  the  "bandgap"  of  CuPc  is  taken  as  1.7  eV,  and  the  "bandgap"  of  PTCDA  is 
2.2  eV.  For  the  purposes  of  this  discussion,  the  highest  occupied  molecular 
orbital  (HOMO)  is  treated  as  the  valence  band,  and  the  lowest  unoccupied 


molecular  orbital  (LUMO)  is  treated  as  the  conduction  band.  The  band  offset, 
AEv,  is  also  shown  in  the  figure.  Note  that,  although  the  HOMO  and  LUMO  are 
somewhat  similar  to  bands  in  conventional  semiconductors,  the  bandwidth 
typical  of  organic  semiconductors  is  ~  1  -  10  meV  [22],  whereas  for  inorganic 
semiconductors  the  bandwidth  is  ~  1  eV.  These  factors  have  profound  effects 
on  the  charge  transport  since  both  the  density  of  states  and  carrier  effective 
masses  are  significantly  different  for  organic  and  inorganic  semiconductors.  We 
note  too,  that  the  band  picture  cannot  be  accurately  applied  to  organic  materials 
whose  carrier  mobility  is  significantly  lower  than  PTCDA  (where  a  hole  mobility 
from  0.1  cm2A/-s  to  1  cm2A/-s  at  room  temperature  is  typical). 

An  alternative  picture  for  the  CuPc/PTCDA  HJ  based  on  molecular 
spectral  levels  is  shown  in  Fig.  7.  This  diagram  phenomenologically  describes 
the  processes  of  heterojunction  transport  using  the  hopping  model,  without 
invoking  band  structure.  In  the  figure,  the  HOMO  and  LUMO  bands  are  shown 
as  functions  of  the  configuration  coordinate,  Q  (where  the  subscripts  C  and  P 
refer  to  CuPc  and  PTCDA,  respectively).  Both  the  ground  states  (So)  and  first 
excited  singlet  states  (Si)  are  resolved  into  different  vibronic  levels  (labeled  0, 
1,  2, ...).  As  indicated  in  Fig.  7,  the  ground  states  (So)  of  CuPc  and  PTCDA  are 
offset  by  an  energy  of  AEv  This  is  the  energy  required  to  transfer  a  hole  from  a 
CuPc  to  a  PTCDA  molecule  across  the  heterojunction  interface,  and  hence  is 
just  the  band  offset  energy  obtained  from  the  C*V  and  J-V  transport 
measurements.  While  this  molecular  level  diagram  is  useful  for  understanding 
transport  across  the  heterointerface,  it  fails  to  fully  describe  the  transport  within 
the  layers.  In  that  case,  the  band  model  is  more  useful  for  interpreting  the 
macroscopic  electrical  characteristics  of  organic  HJs.  Whichever  model  is  used, 
however,  our  results  suggest  that  the  organic  HJ  provides  an  energy  barrier  to 
charge  transport  similar  to  that  observed  in  inorganic  semiconductor  HJs. 

5.  Optical  properties  organic  multiple  quantum  wells 

In  the  previous  section,  we  demonstrated  that  semiconductor 
heterojunctions  can  be  formed  by  layering  two  different  crystalline  organic  thin 
films,  and  that  the  energy  barrier  at  the  interface  controls  the  flow  of  charge 
across  the  HJ.  By  sandwiching  an  ultra-thin,  small  bandgap  organic  layer 
between  two  large  bandgap  organic  layers,  organic  quantum  well  structures 
can  then  be  realized.  Due  to  the  reduced  dimensionality  in  such  structures, 
charge  carriers  are  confined  between  the  barriers,  hence  strongly  affecting  the 
optical  properties  of  the  resulting  materials.  In  the  following,  we  describe  the 


optical  properties  of  crystalline  organic  MQWs  consisting  of  alternating,  ultra- 
thin  layers  (between  10  A  and  200  A  thick)  of  PTCDA  and  NTCDA  grown  by 
OMBD.  A  schematic  diagram  of  such  a  quantum  well  structure  is  shown  in  Fig. 
8  where  PTCDA  is  the  well  layer  sandwiched  between  adjacent  NTCDA  (with  a 
bandgap  energy  of  3.1  eV)  barrier  layers. 

The  room  temperature  optical  absorption  and  low  temperature  (20  K) 
photoluminescence  bands  of  a  5  period,  symmetric  PTCDA/NTCDA  organic 
MOW  sample  with  an  individual  layer  thickness  of  40  A  is  shown  in  Fig.  9.  The 
characteristic  absorption  spectra  of  both  PTCDA  and  NTCDA  are  indicated  in 
the  figure.  The  absorption  bands  are  due  to  exciton  generation  in  the 
corresponding  crystalline  layers,  and  are  relatively  broad  due  to  the  strong 
exciton-phonon  interaction  characteristic  of  organic  crystals  [22].  The  lowest 
energy  singlet  exciton  absorption  peak  in  PTCDA  (indicated  by  the  left-hand 
arrow)  shifts  to  higher  energy  as  the  layer  thickness  is  decreased,  as  shown  by 
the  data  points  in  Fig.  10.  No  apparent  shifts  in  the  higher  energy  exciton  lines 
were  observed,  possibly  as  a  result  of  the  broad  nature  of  these  absorption 
bands. 

Since  the  electronic  polarization  energy  is  large  (~  1  eV)  in  crystalline 
organic  materials,  it  plays  a  vital  role  in  determining  the  exciton  energy.  Thus 
we  might  expect  the  presence  of  NTCDA  molecules  to  influence  the  PTCDA 
exciton  energy,  as  observed  in  our  MOW  samples.  An  alternative  explanation 
for  the  blue  shift  is  due  to  quantum  confinement  of  excitons  in  a  quantum  well 
structure.  We  now  discuss  both  models  in  detail. 

The  change  of  polarization  energy  of  PTCDA  due  to  the  proximity  of 
NTCDA  molecules  in  closely  spaced  layers  can  be  calculated  using 
perturbation  theory.  A  derivation  of  the  perturbation  to  the  exciton  energy  due  to 
polarization  effects  is  given  in  Appendix  C.  The  results  of  the  calculation  show 
that  the  first  order  energy  shift  to  the  polarization  energy  increases  linearly  with 
the  layer  thickness,  t.  This  effect,  while  being  very  short  range,  can  actually 
increase  the  exciton  energy  with  increasing  well  width,  contrary  to  our 
observations.  Hence,  we  can  rule  out  this  first  order  "energy  offset"  effect.  On 
the  other  hand,  the  second  order  polarization  energy  correction  is  proportional 
to  t"^.  As  shown  in  Fig.  10,  the  data  are  not  consistent  with  the  results  of  the 
calculation  since  their  energy  dependence  is  considerably  smaller  than  even 
t*2.  Hence,  we  conclude  that  the  polarization  model  does  not  adequately 
explain  the  observed  blue  shift  in  the  exciton  line. 


Quantum  confinement  provides  an  alternative  explanation  for  the 
observed  blue  shift  in  the  absorption  spectrum  [1,2].  Here,  we  use  the 
variational  method  to  evaluate  the  binding  energy  of  an  exciton  in  a  quantum 
well.  The  Schrodinger  equation  for  an  exciton  can  be  solved  using  the  trial 
wavefu  notion: 

Vtot  =  Ve(Ze)Vh(Zh)Vl  s(p  2)  (6) 

where  \j/is(p.z)  is  the  hydrogenic  Is  wavefunction,  p  =  x2  +  y2,  zq  and  Zh  are 
coordinates  in  the  z-direction  for  the  electron  and  hole,  respectively,  Ve{ze)  and 
\|/h(zh)  are  the  exact  solutions  to  the  finite  square  well  problem,  and  z  =  Ze  -  Zh- 
In  Eq.  6,  the  hydrogenic  state  has  the  form: 

Vis  =  exp  [  -  (p2/a2  +  z2/(52)1/2]  (7) 

Here,  a  is  a  variational  parameter  in  the  trial  solution,  and  (3  is  set  equal  to  ao, 
where  ao  is  the  exciton  Bohr  radius  in  the  bulk  crystal.  For  crystalline  organic 
MQW  structures,  typically  ao  «  Lz,  where  Lz  is  the  quantum  well  width.  In  this 
case,  therefore,  it  is  reasonable  to  assume  P  =  ao,  with  a  being  the  only 
variation  parameter. 

The  parameters  used  in  this  model  to  fit  the  data  in  Fig.  10  are  listed  in 
Table  2.  Here,  the  sum  of  Ve  and  Vh  is  equal  to  the  difference  in  energy  gaps  of 
PTCDA  and  NTCDA.  Of  all  the  parameters  used  in  the  calculation,  the  results 
are  most  sensitive  to  the  choice  of  mj,  and  the  ratio,  Ve/Vh-  Using  these 
parameters,  a  good  fit  to  the  experimental  data  is  obtained  for  mh  =  0.1 8mo, 
where  mo  is  the  electron  rest  mass.  A  small  value  of  mh  and  a  relatively  large 
value  of  me  were  chosen  in  the  calculation  since  PTCDA  is  predominantly  a 
hole-transporting  material.  The  effect  of  using  different  values  of  mh  are  also 
shown  for  comparison  in  Fig.  10  indicating  the  sensitivity  of  the  fit  to  this 
parameter.  Note  that  our  calculation  implies  that  Ve  »  Vh.  This  is  consistent 
with  electrical  measurements  made  for  PTCDA/NTCDA  heterojunctions  which 
do  not  exhif*  significant  rectification  of  the  hole  current  at  temperatures  ranging 
from  room  temperature  to  about  90  K.  This  indicates  that  the  energy  barrier  to 
hole  transport  is  less  than  100  meV.  A  best  fit  to  the  data  is  obtained  assuming 
Vh  =  50  meV,  in  accordance  with  these  findings.  The  calculated  values  of  a  are 
also  shown  in  Fig.  10.  Here,  a  increases  monotonically  with  increasing  well 
width,  and  asymptotically  approaches  a©  =  12  A. 

The  effect  of  quantum  confinement  of  excitons  in  organic  MQWs  is  further 
evident  from  exciton  lifetime  measurements,  where  we  have  found  that  the 
exciton  lifetime  decreases  with  decreasing  layer  thickness.  For  example,  at  Lz  = 


200  A,  the  radiative  lifetime  is  10.8  ±  0.5  ns,  decreasing  to  5.7  ±  0.5  ns  for  Lz  = 
14  A.  In  previous  work,  the  decrease  in  lifetime  in  PTCDA/NTCDA  MQWs  has 
been  attributed  to  the  shrinkage  of  exciton  volume  due  to  quantum  confinement 
[2].  This  phenomenon,  which  is  often  observed  in  inorganic  semiconductor 
heterostructures,  results  from  an  increased  overlap  of  the  electron  and  hole 
wavefunctions,  thus  enhancing  the  probability  for  exciton  recombination.  It  can 
be  shown  [2]  that  the  shrinkage  in  exciton  volume  can  be  used  to  quantitatively 
explain  the  exciton  lifetime  data  using  the  parameters  given  in  Table  2,  and 
taking  the  wavefunction  given  in  Eq.  6  and  7. 

The  most  important  conclusion  drawn  from  the  above  experiments  is  that 
we  have  observed  effects  of  exciton  quantum  confinement  in  crystalline  organic 
multiple  quantum  wells  for  the  first  time.  From  the  results  of  variational 
calculations,  the  exciton  Bohr  radius  in  PTCDA  is  determined  to  be  ~12  A, 
implying  that  the  exciton  wavefunction  extends  to  3  -  4  nearest  neighbors.  We 
can  therefore  conclude  that  excitons  in  the  crystalline  organic  materials  studied 
here  are  Wannier-like  [16]  rather  than  Frenkel-like. 

6.  Giant  anistropies  in  crystalline  organic  films  deposited  by 
OMBD 

We  have  measured  the  indices  of  refraction  and  dielectric  constants 
along  different  directions  in  thin  films  of  the  crystalline  organic  semiconductor; 
3,4,  9,  10  perylenetetracarboxylic  dianhydride.  The  films  were  deposited  via 
organic  molecular  beam  deposition,  resulting  in  highly  ordered,  "quasi- 
epitaxial"  films.  Due  to  inherent  asymmetries  in  the  molecular  crystal  structure, 
the  ordering  of  the  films  results  in  giant  anistropies  in  their  dielectric  properties. 
For  example,  the  index  of  refraction  measured  at  a  wavelength  of  1  =  1 .064  pm 
in  the  direction  perpendicular  to  the  substrate  plane  is  ni  =  1.36  ±  0.01, 
whereas  parallel  to  the  plane,  ny  =  2.017  ±  0.005,  resulting  in  an  index 
difference  of  An  =  0.66.  Furthormore,  the  dielectric  constant  of  the  films 
measured  between  100  Hz  and  10  fvlhz  is  ei=  1.9  ±  0.1  and  £||  =  4.5  ±  0.2.  To 
our  knowledge,  these  are  the  largest  index  anisotropies  ever  measured  for  thin 
films.  In  this  section,  we  also  discuss  a  guided  wave  polarization-selective 
device  which  takes  advantage  of  the  large  dielectric  anisotropies  characteristic 
of  the  thin  organic  films. 

One  would  expect  that  the  dielectric  properties  of  highly  ordered  films 
would  also  have  a  considerable  degree  of  anisotropy  along  different  crystalline 
directions.  The  molecules  form  planar  stacks,  with  a  distance  of  3.21  A  between 


molecules  in  the  stack.  Typically,  the  films  consist  of  such  stacks  tilted  at  p  =  1 1° 
from  the  substrate  normal.  The  long  axes  of  molecules  in  adjacent  stacks  are 
rotated  by  90°  to  the  long  axis  of  a  molecule  in  the  stack  at  the  center  of  the  unit 
cell.  The  off-resonance  dielectric  constant  of  a  material  at  frequency,  to,  is 
given  by: 

>  (8) 

where  N  is  the  molecular  number  density,  coj  is  the  dipole  transition  frequency, 
q  is  the  electronic  charge,  and  m*  is  the  electron  effective  mass.  The  oscillator 
strength  of  the  dipole  transition  from  level  o  to  j  is  given  by  fj  =  2m*hQ)|xoj|2/h2, 
where  Xoj  is  the  electron  position  expectation  value,  and  h  is  Planck's  constant 
divided  by  2k.  Assuming  that  the  largest  contribution  to  the  dipole  moment  of 
the  molecule  is  due  to  the  p-orbitals,  and  that  the  electron  within  a  given  orbital 
is  completely  delocallized,  then  we  can  make  the  approximation  that  fji/fjn  ~ 
|xojil^/|xoj||l^  ~  d2/L2,  Here,  d  =  3.2  A  is  the  extent  of  the  7t-orbital  system 
perpendicular  to  the  molecular  plane,  and  L  is  its  extent  in  the  plane,  which  is 
approximately  equal  to  the  length  of  the  perylene  molecular  core  of  PTCDA  (i.e. 
L  =  6.9  A).  Hence,  provided  that  there  is  perfect  crystalline  alignment 
throughout  the  thin  film,  we  can  expect  an  anisotropy  in  e  between  directions 
perpendicular  and  parallel  to  the  thin  film  plane  (and  hence  approximately 
perpendicular  and  parallel  to  the  molecular  stacking  axis)  of  (ei-1)/(E||-1)  "• 
d2/L2  ~  0.22.  This  value  is  only  an  approximation,  since  we  are  assuming  that 
m*  is  isotropic  [22],  and  that  the  dipole  moment  is  due  to  a  completely 
delocallized  electron  in  the  extended  p-system.  Nevertheless,  it  is  shown  below 
that  this  value  is  close  to  the  measured  anisotropy  in  e  for  highly  ordered,  quasi- 
epitaxial  films  of  PTCDA. 

To  measure  the  dielectric  constant  along  different  film  directions, 
capacitors  oriented  along  the  various  film  axes  were  fabricated  on  quartz 
substrates.  For  measuring  the  capacitance  along  the  in-plane  direction, 
interdigitated  patterns  of  Cr-Au  contacts  were  employed.  Here,  a  thin  layer  (~ 
150  A)  of  Cr  was  deposited,  followed  by  a  0.7  pm  thick  Au  layer  which  forms  a 
current-blocking  contact  to  PTCDA  [15].  Next,  the  contacts  were  patterned  into 
2  mm  long  interdigitated  electrodes  with  2pm  spacing  and  finger  width.  In 
addition,  some  electrode  sets  were  positioned  perpendicular  to  other  sets  to 


determine  if  there  were  significant  asymmetries  in  the  capacitance  in  the  two  in¬ 
plane  directions.  Next,  purified  PTCDA  was  deposited  via  organic  molecular 
beam  deposition  (OMBD)  onto  a  quartz  substrate  which  was  maintained  at  a 
temperature  of  approximately  90K  during  growth.  The  PTCDA  thickness  was 
0.9  pm,  filling  the  spaces  between  the  electrode  fingers.  As  has  been 
discussed  previously  [1],  deposition  under  these  conditions  ensures  that  the 
resulting  film  is  single  crystalline,  independent  of  the  substrate  material 
employed. 

The  capacitance  measured  between  100  Hz  and  10  MHz  for  one  set  of 
electrode  pairs  using  this  in-plane  geometry  is  shown  in  Fig.  11.  The 
capacitance  is  independent  of  frequency  over  the  measurement  range, 
indicating  that  the  films  are  relatively  free  of  traps.  Traps  tend  to  affect  the 
capacitance  when  their  emission  rate  is  comparable  to  the  test  frequency. 

The  capacitance  of  a  co-planar  strip  line  lying  on  a  quartz  substrate  with 
dielectric  constant  es  is  given  by  [23]:  C  =  (es  +  eF)eo[FK(k)]l/2,  where  £o  is  the 
permittivity  of  free  space,  I  is  the  finger  length,  F  =  16  is  the  number  of  finger 
pairs,  K  is  the  ratio  of  two  complete  elliptic  integrals  of  the  first  kind,  and  k  is  a 
geometrical  parameter  related  to  the  finger  width  and  spacing.  Here,  we 
assume  quartz  fills  one  half  plane,  and  a  material  with  dielectric  constant  ep  the 
other.  The  ratio  of  capacitances  of  the  pattern  measured  before  (Cwo)  after 
(Cw)  deposition  of  the  PTCDA  thin  film  is  then  simply: 

Cwo^Cw  =  ^  V  (^s  +  £||)  (9)- 

We  note  that  the  small  amount  of  electric  field  penetration  above  the  PTCDA 
film  can  lead  to  deviations  from  Eq.  (9).  However,  since  the  PTCDA  thickness 
was  larger  than  that  of  the  contact  fingers,  this  error  leads  to  an  insignificant 
underestimate  of  e||.  We  obtain  in-plane  dielectric  constants  for  PTCDA  of  £((  = 
4.6±  0.2  and  4.3±  0.2  for  the  two  perpendicular  orientations  of  the  patterns. 
Since  these  values  are  within  each  other’s  error  limits,  we  assume  that  the 
dielectric  constants  along  the  two  directions  is  not  significantly  different,  leading 
to  e||  =  4.5±  0.2. 

To  measure  ex,  a  parallel  plate  geometry  was  used,  where  the  back 
contact  was  formed  using  a  broad  Cr-Au  contact  deposited  on  glass.  This  was 
followed  by  deposition  of  a  0.9  pm  thick  layer  of  PTCDA.  Top  contact  to  the 
PTCDA  was  made  using  circular  Au  pads  of  area  2.5  x  10"^  cm^.  Using  this 
geometry,  it  is  found  that  ex  =  1.9  ±  0.1  over  the  same  frequency  range  as  that 


used  to  measure  e||.  From  these  data  we  obtain  (ej_-  l)/(e||  -  1)  =  0.26  which  is 
close  to  our  theoretical  estimate  of  0.22  for  perfectly  oriented  crystalline  films. 

In  earlier  measurements  of  ex  for  PTCDA  deposited  under  conditions 
leading  to  polycrystalline  growth,  it  was  found  that  [15]  ex=  3.6  ±  0.4.  Assuming 
a  homogeneous  medium  was  being  analyzed  in  that  case,  then  the  spatially 
averaged  dielectric  constant  is  e  =  (2eii  +  ex)/3.  Using  the  values  of  e  obtained 
in  this  work,  we  find  e  =  3.6.  This  is  in  agreement  with  the  former  value, 
indicating  that  the  films  deposited  on  room  temperature  substrates  using 
standard  vacuum  sublimation  techniques  are  structurally  randomized,  whereas 
deposition  via  OMBD  on  cold  substrates  results  in  films  which  have  nearly 
perfect  crystalline  structure  across  broad  areas. 

Since  £«  =  n2  ,  where  n  is  the  refractive  index  of  the  film,  and  £«  is  the  thin 
film  dielectric  constant  at  optical  frequencies,  we  also  expect  to  observe 
asymmetries  in  n  as  a  result  of  the  ordering  of  the  quasi-epitaxial  films  of 
PTCDA.  In  earlier  work  [11],  we  measured  an  in-plane  refractive  index  of  n||  = 
2.017  ±  0.005.  In  this  work,  therefore,  we  measure  only  nA. 

Determination  of  ni  proceeds  by  measuring  the  reflectivity  from  the  thin 
film  as  a  function  of  beam  incident  angle  using  several  film  thicknesses  and 
light  wavelengths.  Considering  an  isotropic  homogeneous  thin  film 
sandwiched  between  two  semi-infinite  media,  the  reflectivity  is  given  by  [24]: 

rie4-r23f2ri2r2;POs2|3 

1+ri2r23+2ri2r23pos2|3 

where  p  =  27cn2tcos92/^.  Here,  n2  is  the  refractive  index  measured  along  the 
refractive  angle  02  in  the  thin  film  of  thickness,  t.  Also,  ri2  (r23)  is  the  reflectivity 
at  the  interface  between  medium  1  (3)  and  medium  2.  The  reflectivity  maxima 
and  minima  for  a  given  index,  n2,  occur  at  angles  02  which  satisfy  n2t  = 
m^4cos02,  where  m  is  the  integer  order  of  the  extremum.  If  the  film  thickness  is 
accurately  determined,  the  refractive  index  can  thence  be  obtained  by 
measuring  the  angles  corresponding  to  reflectivity  extrema  as  predicted  by 
these  equations. 

To  extend  this  technique  to  anisotropic  thin  films,  the  incident  light 
polarization  is  made  parallel  to  the  plane  of  incidence;  i.e.  it  is  a  TM-polarized 
wave.  Furthermore,  using  birefringence  measurements  discussed  previously 
[1],  the  optical  b-axis  is  located,  and  is  also  oriented  parallel  to  the  incident  light 
polarization  vector.  When  the  incident  angle  (0i)  of  the  beam  to  the  film  is 


changed,  this  is  equivalent  to  a  rotation  of  the  angle  62  in  an  index  ellipsoid  with 
axes  of  length  ny  =  nb  and  nj.,  as  shown  in  the  inset  in  Fig.  12.  Here,  nb  is  the 
index  along  the  b-axis.  From  this  figure,  it  can  be  shown  that  the  index 
r\2^  =  n|i2  +  (1  -  n||2/nj_2)sin20i  from  which  nj_is  easily  extracted. 

A  typical  measurement  result  made  at  X  =  1 .064  pm  for  a  1 .25  pm  thick 
PTCDA  film  together  with  a  theoretical  calculation  using  Eq.  (10)  is  shown  in 
Fig.  1 2.  The  experimental  setup  consists  of  a  YAG  laser  providing  the  incident 
beam.  A  large  area  detector  is  positioned  near  the  top  film  surface  to  measure 
the  reflected  light  intensity.  The  PTCDA  is  deposited  on  a  quartz  substrate  at 
low  temperatures  to  achieve  crystalline  order.  Also,  the  back  surface  of  the 
quartz  is  frosted  to  minimize  reflections  from  the  quartz/air  interface.  It  can  be 
seen  that  the  positions  of  the  theoretical  extrema  match  those  obtained 
experimentally,  although  the  magnitude  of  the  reflected  signal  is  sometimes 
different  than  predicted  due  to  unwanted  reflections  from  surface  imperfections 
and  the  quartz/air  interface.  Similar  fits  are  obtained  for  a  film  thickness  of  1 .13 
pm,  and  at  X  =  1 .3  pm,  and  the  results  are  listed  in  Table  3.  From  these  data,  we 
obtain  n^  =  1.36  ±  0.01  at  X  =  1.064  pm,  giving  An  =  n||  -  nj.  =  0.66  which 
apparently  is  the  largest  value  of  An  reported  for  thin  films  far  from  their 
absorption  edge  (which  for  PTCDA  lies  at  X  =  6000  A  [1 5]).  Very  little  dispersion 
is  observed  in  n  measured  at  X  =1 .064  pm  and  1 .3  pm,  as  indicated  by  the  data 
in  Table  3. 

It  is  interesting  to  compare  the  asymmetries  measured  for  the  dielectric 
constant  and  the  index  of  refraction.  Ignoring  dispersion  between  the 
frequencies  at  which  these  parameters  were  determined,  we  find  that  e||/ei  =  2.4 
=  (n|i/ni)2  =  2.2.  This  is  truly  a  remarkable  agreement  given  the  widely  different 
conditions  under  which  the  various  measurements  were  made.  Indeed,  we  see 
that  the  low  frequency  dielectric  constant  is  approximately  equal  to  n2  for  each 
of  the  various  film  directions. 

A  waveguide  polarizer  was  fabricated  to  take  advantage  of  the  large 
birefringence  characteristic  of  OMBD-grown  PTCDA  thin  films.  Thus,  2  pm 
wide  PTCDA  rib  waveguides  were  made  of  lengths  of  1.5  and  10  mm.  The 
guides  were  fabricated  by  first  spinning  a  1  pm  thick  film  of  the  photoresist  AZ 
1 400  (with  npR  =  1.61)  onto  the  (1 00)  surface  of  an  InP  wafer  (with  ns  =  3.27  at  X 
=  1 .064  pm).  Next,  2  pm  wide  ridges  were  photolithographically  patterned  in 
the  AZ  1400,  followed  by  deposition  of  approximately  1  pm  of  PTCDA.  The 
regions  where  PTCDA  was  deposited  onto  the  photoresist  ridges  form  the 


guides  to  TE  modes  (where  n))  >  npR),  but  not  to  TM  modes  (where  <  npp). 
Elsewhere,  PTCDA  is  directly  deposited  on  the  large  index  InP  substrate  where 
guiding  does  not  occur.  The  substrates  were  cleaved  along  (110)  planes,  thus 
allowing  for  end-fire  coupling  of  light  into  the  cleaved  facets.  Similar  guides 
have  previously  been  measured  to  have  a  loss  of  <2.5  dB/cm  [6]. 

Due  to  the  selective  guiding  of  TE  over  TM  modes,  the  PTCDA 
waveguides  form  a  strongly  polarizing  medium.  To  measure  the  TM/TE 
extinction  ratio,  X  =  1.064  pm  light  from  a  diode-pumped  YAG  laser  was  coupled 
into  and  out  of  the  guides  using  microscope  objective  lenses.  The  input  light 
polarization  was  continuously  rotated  through  360^  using  a  quarter  wave  plate 
placed  between  the  laser  and  the  sample.  In  addition,  a  Glan-Thomson 
polarizer  (GTP)  with  an  extinction  ratio  >  50  dB  was  placed  in  front  of  the  laser. 
This  set  up,  shown  in  the  inset  of  Fig.  13,  ensures  a  constant  coupling  efficiency 
of  a  highly  polarized  beam  along  any  direction.  The  out-coupled  light  was 
passed  through  a  second  GTP  which  analyzed  the  light  prior  to  being  detected 
using  a  Si  CCD  camera. 

Measurement  results  for  the  10  mm  long  device  are  shown  in  Fig.  13. 
Here,  both  the  TE  and  TM  components  of  the  output  light  beam  are  plotted 
versus  the  input  light  polarization  angle.  An  extinction  ratio  of  TM/TE  <  -48  dB 
was  obtained,  where  the  measurement  accuracy  was  limited  by  the  GTP. 
Similar  results  were  obtained  for  the  1.5  mm  guide,  in  which  an  extinction  ratio 
of  TM/TE  «  -  25  dB  was  measured.  The  accuracy  of  this  measurement  was 
limited  by  light  scattered  over  the  top  of  the  short  guide.  The  propagation  loss 
of  the  short  guide  was  found  to  be  <  0.4  dB. 

In  summary  of  this  part  of  the  work,  extremely  large  anisotropies  in 
dielectric  properties  have  been  observed  in  crystalline  organic  thin  films  of  the 
compound,  PTCDA.  The  anisotropies  in  both  the  dielectric  constant  and  index 
of  refraction  for  these  films  is  a  result  of  their  nearly  perfect  crystalline  order 
achieved  via  OMBD  on  cold  substrates.  Similar  results  are  expected  for  other 
anisotropic  organic  semiconductor  crystals  deposited  using  the  OMBD 
technique.  Furthermore,  a  TE-pass  waveguide  polarizer  with  a  high  TM 
polarization  extinction  and  low  TE  propagation  loss  was  demonstrated.  Since 
such  polarizing  guides  can  be  deposited  on  conventional  semiconductor 
substrates,  these  films  have  uses  such  as  optical  isolators  integrated  with 
semiconductor  lasers  or  amplifiers. 

7.  Resonant  nonlinear  effects  in  crystalline  organic  films 
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In  this  work,  we  demonstrated  all-optical  modulation  in  crystalline  organic 
semiconductor  waveguides  grown  by  the  ultrahigh  vacuum  process  of  organic 
molecular  beam  deposition.  Two  light  beams  with  wavelengths  of  1.06  pm  and 
0.514  pm  were  used  as  the  guided  and  the  pump  light  sources,  respectively.  A 
refractive  index  change  of  5.4  x  lO’^  at  1.06  pm  was  observed  at  a  pump  intensity 
of  1.0  W/cm2.  This  large  nonlinear  effect  is  attributed  to  free  electron-hole  pairs 
produced  by  the  dissociation  of  excitons  generated  by  the  short  wavelength 
beam.  A  carrier  lifetime  of  (17  ±  1)  ps  which  determines  the  modulator  switching 
time  is  in  good  agreement  with  theoretical  predictions.  To  our  knowledge,  this  is 
the  first  observation  of  free  carrier  induced  index  modulation  in  crystalline  organic 
waveguides. 

Charge  transport  in  molecular  semiconductors  is  thought  to  be  initiated  via 
excitonic  transitions  [25].  Here,  free  electron-hole  pairs  are  generated  in  a 
second  order  process  via  dissociation  of  photogenerated  excitons  at  surfaces  or 
bulk  material  impurities.  That  is,  given  a  molecule  in  ground  state  So,  a  dominant 
channel  for  free  carrier  production  due  to  absorption  of  light  of  energy  hw,  is  So  + 
hw  ->  S',  followed  by  S*  +  M  ->  e*  +  h+.  Here,  S*  is  the  excited  singlet  exciton 
state  of  the  molecule,  M  represents  an  impurity  or  interface,  and  e*  and  h+  are  the 
free  electron  and  holes,  respectively.  The  absorbed  light  must  have  a  photon 
energy  greater  than  the  exciton  transition  energy  (typically  1.5-3  eV).  Since  the 
process  is  second  order,  the  quantum  efficiency  (h)  for  electron-hole  generation 
is  usually  <1%.  Assuming  that  there  is  a  refractive  index  change  (An)  associated 
with  the  optical  generation  of  free  carriers  (via  excitons),  the  Drude  model  gives: 

-  e2N 

=  - - ; - 2  (^1) 

2nm  eoo)-^ 

where  e  is  the  electronic  charge,  N  is  the  free  carrier  density  induced  by  the 
optical  beam,  n  is  the  refractive  index  in  the  dark,  m*  is  the  effective  carrier  mass, 
Eo  is  the  permittivity  in  vacuum,  and  w  is  the  light  frequency.  A  steady  state 
solution  of  the  continuity  equation  for  N  gives: 


An  =  n2l  = 


-e^axTiI 

2nm*eo 


(12) 


Here,  a  is  the  absorption  coefficient,  I  is  the  optical  intensity,  x  is  the  free  carrier 
lifetime,  and  h  is  Planck's  constant  divided  by  2k.  Note  that  the  NLO  coefficient, 
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n2,  given  by  Eq.  (12)  is  only  due  to  the  free  carrier  concentration  (not  the  exciton 
population),  and  is  intensity  dependent. 

To  fabricate  the  waveguide  devices  in  which  n2  was  measured,  a  ipm- 
thick,  AZ  1400  photoresist  layer  ,'with  np  =  1.61)  was  first  spun  onto  the  surface  of 
a  cleaned  and  polished  (100)  mP  substrate.  Next,  a  series  of  2  pm  wide  strips 
were  patterned  onto  the  phote.-^esist  along  the  (110)  direction  by  standard 
photolithographic  techniques.  A  1  pm  thick  crystalline  PTCDA  layer  was  then 
deposited  onto  both  the  InP  substrate  and  the  photoresist  ridges  using  the 
ultrahigh  vacuum  (~10'9  Torr)  process  of  organic  molecular  be^m  deposition  [2]. 
Deposition  proceeds  by  heating  a  pre-purified  source  of  PTCDA  to  450  °C  to 
achieve  a  deposition  rate  of  approximately  3  A/s,  while  the  substrate  temperature 
is  maintained  at  90  K.  It  has  previously  been  found  that  thin  film  deposition  under 
these  conditions  results  in  single  crystalline  growth  of  PTCDA  on  the  photoresist 
strips  into  low  loss  (<  2.5  dB),  TE-mode  waveguides  [3].  To  form  waveguide 
facets,  the  InP  substrate  was  cleaved  along  the  (110)  direction.  The  facet 
surfaces  are  sufficiently  smooth  to  allow  for  observation  of  Fabry-Perot 
resonances  as  the  index  of  the  guide  is  thermally  or  optically  varied. 

In  the  experiments,  1.06  pm  wavelength  light  from  a  YAG  laser  was 
coupled  into  and  out  of  the  waveguide  using  two  microscope  objective  lenses 
(see  inset.  Fig.  14).  A  high  speed  lno.53Gao.47As  avalanche  photodetector  (APD) 
was  used  to  detect  the  light  signal.  To  ensure  a  constant  coupling  efficiency,  a  TV 
camera  monitored  the  beam  spot  position  using  a  beam  splitter  placed  in  front  of 
the  APD.  The  waveguide  was  placed  on  a  copper  block  mounted  on  a 
thermoelectric  cooler.  Using  a  thermocouple  and  a  temperature  controller,  the 
temperature  of  the  waveguide  device  was  controlled  to  within  ±0.05  °C. 

The  2pm  wide  PTCDA  rib  waveguide  allowed  for  propagation  of  a  single 
TEqo  mode.  The  waveguide,  with  its  cleaved  facets,  acts  as  a  lossy  Fabry-Perot 
resonator  with  a  transmission  function  of: 


P(5)  = 


(l.R)2e-aL 


( 1  -Re'“T)2+4Re-aLsjn 


(13) 


where  a  is  the  absorption  coefficient,  L  is  the  device  length  (1.62  mm),  R  is  the 
reflection  coefficient  at  the  waveguide  facets  which  is  calculated  to  equal  0.11 
considering  only  Fresnel  reflections,  and  d  is  the  phase  shift  defined  as:  d  =  do  + 
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A5  =  4:tL/X  (n+An).  When  the  temperature.!,  is  varied,  the  transmitted  intensity 
varied  periodically  (see  Fig.  14)  as  is  characteristic  of  a  Fabry-Perot  resonator.  A 
complete  fringe  (two  adjacent  maxir  i)  corresponding  to  a  phase  difference  of  Dd 
=  2k  was  observed  for  a  temperature  change  of  2  °C,  from  which  we  obtain 
(1/n)An/AT  =  (8.1  ±  0.5)  x  Using  Eq.  (13),  however,  a  maximum 

transmission  modulation  of  -50%  was  expected,  compared  to  an  observed 
modulation  of  only  10%.  This  smaller  than  expected  modulation  is  possibly  due 
to  imperfect  waveguide  facets  which  lower  the  finesse  of  the  resonator. 

To  measure  the  all-optical  modulation  characteristics  of  the  guide,  a  pump 
light  beam  of  0.514  pm  wavelength  with  a  diameter  ._‘=1  mm  directly  illuminated 
the  surface  of  the  PTCDA  waveguide,  while  the  1.06  pm  beam  was  end-fire 
coupled  into  the  waveguide,  as  in  the  case  of  An/AT  measurements.  In  order  to 
avoid  thermally  induced  intensity  changes  from  the  pump,  the  device  temperature 
was  precisely  maintained  at  22.1  °C  using  the  thermoelectric  cooler.  The  pump- 
induced  intensity  changes  are  shown  in  Fig.  15.  Here,  the  transmitted  intensity  is 
decreased  nearly  linearly  proportional  to  the  pump  intensity  at  a  temperature  of 
22.1  °C.  In  contrast,  the  transmitted  intensity  due  to  a  temperature  rise  at  22.1  °C 
is  increased  (Fig.  14). 

While  we  suggest  that  the  pump  power  dependence  of  P(5)  is  due  to  an 
optically  induced  index  change,  it  is  nevertheless  important  to  consider  other 
competing  processes.  The  most  significant  contributions  to  An  are  thermally 
induced  index  changes  and  optically  induced  absorption  effects.  To  differentiate 
between  these  various  mechanisms,  we  have  investigated  the  modulation 
behavior  at  very  low  pump  light  intensifies  (<  1  W/cm^)  over  a  wide  temperature 
range  (16  °C  -37  °C).  Three  observations  exclude  the  possibility  of  thermal 
effects:  Firstly,  the  changes  in  transmitted  intensity  under  a  constant  pump  power 
could  either  be  increased  or  decreased,  depending  on  temperature.  However, 
the  sign  of  the  op^'cally  induced  change  of  transmitted  intensity  was  always 
opposite  to  the  thermally  induced  changes  over  the  entire  temperature  range 
explored.  For  instance,  at  a  certain  temperature,  if  the  transmitted  intensity  was 
decreased  due  to  pump  light  illumination,  it  was  increased  due  n  a  temperature 
rise.  This  is  similar  to  observations  for  inorganic  semiconductors  [26]  in  which 
nonlinear  effects  were  attributed  to  index  changes  due  to  excitons. 

A  further  observation  ruling  out  thermal  effects  is  that  the  temperature 
increase  is  calculated  to  be  -0.03  °C  under  a  maximum  pump  light  intensity  of  1 
W/cm2.  This  calculation  assumes  a  steady  ctate  thermal  conductivity  of  0.1 , 0.007 


and  0.8  W/cm-K  for  PTCDA,  AZ1400  and  InP,  respectively.  Such  a  small 
temperature  rise  does  not  significantly  contribute  to  the  intensity  modulation,  as 
inferred  from  Fig.  14.  Finally,  the  transient  thermal  relaxation  time  was  calculated 
to  be  >13  ms,  which  is  three  orders  of  magnitude  larger  than  the  measured 
modulation  time  of  17  ps,  as  shown  below.  A  very  low  amplitude  "tail"  in  the 
modulation  response  with  a  time  constant  of  ~6  ms  was  observed  at  the  highest 
pump  intensities.  We  attribute  this  latter  time  constant  to  be  characteristic  of  very 
small  thermal  effects,  whereas  the  shorter  time  constant  results  from  significantly 
larger  carrier  lifetime  effects. 

The  second  "parasitic"  effect  to  be  considered  is  pump  induced  absorption 
changes,  Aa.  Using  Eq.  (13),  it  can  be  shown  that  the  change  of  transmitted  light 
intensity,  AP,  in  the  limit  of  aL«1 ,  is  given  by:  AP  =  -A  Aa,  where  A  is  a  constant. 
Since  the  waveguide  used  in  our  experiment  had  aL  =  0.092,  this  expression  is 
valid.  Thus,  if  the  optically  induced  intensity  change  is  due  to  Aa,  we  should  see 
a  constant  decrease  of  the  transmitted  intensity  under  a  constant  pump  light 
illumination,  independent  of  temperature.  However,  as  noted  above,  the  sign  of 
the  observed  change  varied,  dependent  on  temperature.  In  addition,  the 
absorption  change  due  to  the  optically  induced  plasma  is  given  by  Aa/An  =  2/ct, 
where  i  is  the  carrier  lifetime  and  c  is  the  velocity  of  light.  For  PTCDA,  Aa/An  = 
4.2  X  10'6  cm-^and  absorption  changes  are  insignificant.  We  conclude, 
therefore,  that  the  modulation  is  due  to  an  intensity  dependent  refractive  index, 
rather  than  a  change  in  the  absorption  coefficient. 

Since  both  optically  and  the  thermally  induced  modulation  is  caused  by 
the  changes  in  refractive  index,  the  thermal  measurements  can  be  used  to 
determine  An  =  n2l  [27].  That  is,  since  the  thermal  measurements  indicated  that 
an  output  intensity  change  of  10%  corresponds  to  a  p-phase  change  (Fig.  14), 
then  the  2%  decrease  in  transmitted  intensity  induced  by  a  pump  power  of  I  = 
0.98  W/cm2  corresponds  to  A5  =  n/5.  Using  the  relationship:  A5  =47t(n2l)L7>.  we 
get  n2  =  5.4  x  10'5  cm2/w,  where  L'  =  1  mm  is  the  interaction  length.  Here,  22.1 
°C  was  chosen  as  the  measurement  temperature  since  the  optically  induced 
intensity  change  is  nearly  linear  over  a  small  range  of  pump  intensity  (<1  W/cm2). 
Of  course,  the  measurement  can  be  made  equally  well  at  other  temperatures. 

That  n2  can  be  explained  by  the  presence  of  optically  generated  free 
carriers  was  verified  By  comparing  the  measured  carrier  lifetime  to  theoretical 
predictions  we  verified  that  n2  is  indeed  due  to  an  optically  generated  plasma. 
Here,  t  was  directly  measured  from  the  transient  response  time  of  the  1.06  pm 


wavelength  beam  due  to  a  pulsed  0.514  pm  pump  beam.  Using  a  measurement 
system  with  a  response  time  of  <  1  ps,  the  rise  and  fall  times  of  the  1.06  pm  beam 
were  found  to  be  18  ps  and  16  ps,  respectively  (Fig.  16).  For  comparison,  from 
the  steady  state  measurements  of  na,  using  Eq.  (2)  we  obtain  t  =  (14  ±  1)  ps, 
which  is  very  close  to  the  measurement.  The  constants  used  in  the  calculation  for 
PTCDA  are:  a  =  4.6  x  10^  cm*'*  at  X  =  0.514  pm  as  obtained  from  absorbance 
measurements,  a  hole  effective  mass  for  PTCDA  of  m*  =  0.18  me  (where  me 
electron  rest  mass)  [2],  and  a  quantum  efficiency  of  ti  =  0.01  [28].  The  presence  of 
electrons  was  ignored  since  their  effective  mass  is  at  least  an  order  of  magnitude 
larger  than  that  of  holes  [2].  In  this  case,  Eq.  (12)  implies  that  their  effect  on  An  is 
negligible. 

In  Table  4  we  list  the  resonant  nonlinear  properties  of  some  typical 
inorganic  semiconductors,  along  with  one  organic  polymer,  and  PTCDA.  From 
this  table,  r\2  of  PTCDA  is  one  of  the  highest  observed  to  date.  However,  the 
response  time  (t)  is  comparatively  long.  To  reduce  the  carrier  lifetime  for  high 
speed  applications,  organic  multiple  quantum  well  structures  [12]  are  potentially 
useful.  However,  decreasing  the  carrier  lifetime  might  also  reduce  h,  thereby 
ultimately  decreasing  r\2.  Alternatively,  one  can  reduce  x  by  electrically  injecting 
holes  from  contacts,  and  then  sweep  them  out  using  an  applied  electric  field. 

In  the  conclusion,  we  report  the  apparently  first  observation  of  free  carrier- 
induced  nonlinear  effects  in  crystalline  organic  semiconductors.  An  intensity- 
dependent  refractive  index  of  na  =  5.4  x  l0-5cm2/w  has  been  measured,  and  is 
among  the  highest  values  reported  for  practical  optical  device  materials. 

8.  Future  Prospects  for  Quasi-Epitaxial  Organic  Structures 

Our  results  for  the  charge  transport  and  optical  properties  of  these  novel 
layered  structures  provide  an  unusually  large  range  of  possibilities  for  materials 
and  device  physicists.  That  is,  crystalline  organic  compounds  have  proven  to 
be  an  exciting  class  of  materials  from  the  standpoint  of  providing  a  model 
system  for  understanding  fundamental  transport  properties  in  van  der  Waals 
solids.  In  addition,  their  environmental  stability,  compositional  versatility  and 
combination  of  physical  properties  which  bridge  the  materials  spectrum 
between  semiconductors  and  insulators  opens  the  prospect  for  a  broad  new 
class  of  engineered  materials  for  use  in  future  generation  optoelectronic  device 
applications.  To  summarize  our  findings  to  date,  crystalline  organic  compounds 
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have  several  properties  which,  taken  together,  make  them  unique  in  the  known 
spectrum  of  materials: 

1 .  Many  crystalline  materials  bonded  by  vdW  forces  can  be  grown  into  ordered 
"quasi-epitaxial"  films  on  substrates  or  other  films  without  regard  to  lattice 
matching  requirements. 

2.  The  anisotropic  crystalline  structure  of  the  films  results  in  giant  anisotropies 
in  their  conductive  and  dielectric  properties.  As  an  example,  PTCDA  films  have 
a  conductivity  anisotropy  of  ~  10®,  and  a  birefringence  of  An  =  0.66.  These 
properties  are  advantageous  for  a  broad  range  of  devices.  Recently, 
polarization-selective  waveguides  employing  the  large  birefringence  of  PTCDA 
have  been  demonstrated  [3]. 

3.  Heterojunctions  consisting  of  an  organic  film  in  contact  with  an  inorganic 
semiconductor  substrate  exhibit  electrical  and  optical  characteristics  which  are 
often  the  optimal  combination  of  these  properties  for  the  particular  HJ  materials. 

4.  Heterojunctions  consisting  of  two  organic  films  in  contact  with  each  other 
often  exhibit  electrical  and  optical  properties  analogous  in  many  ways  to  fully 
inorganic  heterojunctions. 

5.  Evidence  for  quantum  confinement  in  very  thin  layer  organic  multiple 
quantum  well  structures  has  furthered  our  understanding  of  the  nature  of 
excitons  in  crystalline  organic  materials.  Given  that  excitons  in  these  materials 
have  been  the  focus  of  study  for  over  40  years  [33,34],  in  many  ways  they  are 
the  key  to  our  understanding  the  physics  of  vdW  solids  in  general.  In  addition, 
the  observation  of  quantum  confinement  suggests  that  many  devices  based  on 
organic  MQWs  are  becoming  a  realistic  possibility.  One  example  of  such  a 
device  would  be  an  organic  MOW  optical  modulator  which  could  be  made  to 
cover  all  regions  of  the  UV,  visible  or  near  IR  spectral  regions,  depending  on  the 
choice  of  organic  compounds  employed. 

6.  The  large  exciton-phonon  coupling  strength  characteristic  of  organic 
semiconductors  can  lead  to  exciting,  very  large  nonlinear  optical  effects.  For 
example,  Lam  and  co-workers  [35]  predicted  that  organic  MQWs  have  large 
^(3)  effects,  and  are  therefore  potentially  useful  for  a  broad  range  of 
applications  which  are  unattainable  using  the  polymers  and  photorefractive 
materials  studied  today.  They  have  shown  that  quasi-epitaxially  grown 
crystalline  organic  thin  films  and  MQWs  with  large  Franck-Condon  (FC)  shift 
energies  exhibit  an  intensity-dependent  absorption  and  index  of  refraction  at 
relatively  low  input  powers.  Here,  the  FC  shift  is  the  energy  difference  between 
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the  lowest  energy  absorption  and  the  highest  energy  luminescence  spectral 
peaks,  and  is  proportional  to  the  square  of  the  exciton-phonon  coupling 
strength.  As  shown  by  the  arrows  in  Fig.  10,  the  PTCDA/NTCDA  MQW  structure 
with  40  A  layer  thickness  has  FC  =  4050  cm-"*,  which  to  our  knowledge  is  the 
largest  yet  observed  for  organic  semiconductors. 

To  understand  the  nonlinear  optical  properties  of  organic  MQWs,  the 
Hamiltonian  for  organic  MQWs  under  illumination  can  be  written  [35]: 

H  =  h((0x  -  coi)a+a  -  hXQa+a  +  hcoob+b  -  pa+E  -  paE*  (14) 

where  cox,  toi,  and  coq  are  the  exciton,  incident  light  and  phonon  frequencies, 
respectively,  and  X  is  the  exciton-phonon  coupling  constant,  p  =  qao  is  the 
exciton  dipole  moment,  E  is  the  external  electrical  field,  and  a+a  and  b+b  are  the 
exciton  and  phonon  populations,  respectively.  Also,  Q  =  b+  +  b  is  the  pnonon 
amplitude. 

The  origin  of  the  optical  nonlinearities  in  organic  MQWs  lies  in  the 
second  term  of  Eq.  14,  which  is  proportional  to  the  exciton-phonon  coupling 
strength,  X,  given  by: 

X  =  V  ci}q(FC) 

The  second  term  in  Eq.  14  "renormalizes"  the  exciton  energy,  (Ox,  by  -hXQ.  As 
the  incident  optical  field  intensity  increases,  the  exciton  population  increases, 
resulting  in  a  decrease  in  the  energy  of  the  exciton  absorption  peak.  Such  an 
effect  leads  to  an  intensity  dependent  absorption,  and  hence  an  intensity 
dependent  index  of  refraction  through  the  Kramers-Kroenig  effect.  Note  that  if 
the  incident  light  beam  is  detuned  from  the  exciton  peak  on  the  low  side  of  cox  by 
a  factor  of  coi  -  cox,  then  the  absorption  peak  energy  is  shifted  toward  the  pump 
light  energy,  which  decreases  the  detuning.  This  increases  the  absorption  with 
increasing  incident  light  intensity,  leading  to  optical  bistability. 

The  two-wave  mixing  gain  for  a  weak  probe  beam  interacting  with  a 
strong  pump  beam  for  PTCDA/NTCDA  MQW  structures  predicted  by  the  time 
dependent  solution  of  Schrodinger's  Equation  and  using  Eq.  14  is  shown  in 
Fig.  17.  The  horizontal  axis  corresponds  to  the  (Rabi  frequency)^  which  is 
proportional  to  the  incident  optical  power  in  the  pump  beam,  and  the  absorption 
or  gain  in  the  low  intensity  beam  is  plotted  on  the  vertical  axis.  In  the  plot,  the 
absorption  becomes  negative  (i.e.  gain)  for  a  Rabi  frequency  of  V60, 
corresponding  to  a  power  density  of  ~  100  W/cm^  for  PTCDA.  Also,  a  steep 


increase  in  exciton  population  (hence  absorption)  occurs  at  a  Rabi  frequency  as 
low  as  20.  Thus,  strong  nonlinear  absorption  is  expected  at  power  densities 
well  below  50  W/cm2.  The  results  of  the  calculation  indicate  that  crystalline 
organic  materials  with  a  large  FC  shift  exhibit  unusual  optical  nonlinearities 
which  are  unattainable  with  other  classes  of  materials.  With  the  ability  to  grow 
MQW  structures  and  tailor  their  optical  properties,  the  use  of  this  new  class  of 
engineered  materials  for  photonic  applications  can  thus  be  widened. 
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APPENDIX  A 


TRIP  REPORT 


FRONTIER  FORUM  ON  ORGANIC  MBE  AND  STM 

Riken  Institute  of  Physical  and  Chemical  Research 
Wako,  Saitama,  Japan 
March  11  -  13,  1991 

S.  R.  Forrest 

University  of  Southern  California 


This  topical  conference,  or  ''forum”  is  held  semi-annually  by  the 
Nonlinear  Optics  Group  at  Riken.  This  group  is  directed  by 
Professor  A.  Garito  at  the  University  of  Pennsylvania,  with  local 
direction  provided  by  Professor  Sasabe. 

Riken  is  a  very  forward-looking  and  respected  basic  research 
organization  in  Japan.  Research  carried  out  at  Riken  is  done  by 
both  permanent  and  temporary  Japanese  staff,  and  by  foreign 
visitors.  Indeed,  even  some  of  the  groups  are  directed  by  non- 
Japanese,  such  as  Professor  Garito.  Many  of  the  temporary  Japanese 
researchers  are  on  loan  from  universities  and  industry,  and  hence 
Riken  works  very  closely  with  both  sectors  of  the  Japanese  research 
establishment.  In  many  ways,  what  happens  at  Riken  is  picked  up 
and  followed  closely  by  Japanese  industry.  There  does  not  appear 
to  be  any  organization  in  the  U.S.  which  is  comparable  to  Riken. 

CNRS  in  France  is  somewhat  similar  in  structure  and  goals. 

The  forum  was  attended  by  over  100  visitors  mostly  from  Japanese 
industry.  Papers  on  both  MBE  and  STM  were  presented  from  Japanese 
researchers  at  Riken,  in  academia  and  industry.  Invited  papers  by 
foreign  researchers  from  DuPont,  the  University  of  Munich,  the 
University  of  Pennsylvania,  and  the  University  of  Southern 
California  were  also  presented.  A  few  invited  papers  from  Europe 
and  the  U.S.  were  also  cancelled  due  to  the  Persian  Gulf  war.  A 
full  program  of  the  conference  with  abstracts  is  attached. 

Taken  together,  a  quality  of  the  presentations  was  superb,  with  a 
great  deal  happening  in  the  field  of  organic  MBE  and  STM.  I  will 
give  a  few  of  the  highlights  below: 

Organic  MBE 

Virtually  all  of  the  Japanese  OMBE  research  centers  around  the 
growth  of  the  phthalocyanines  (Pc's)  on  alkali  halide  or  metal 
dichalcogenide  (mdc)  (e.g.  MoSe2)  substrates.  The  work  basically 
follows  the  lead  of  Professor  Koma  at  University  of  Tokyo  who,  for 
the  last  three  years,  has  concentrated  on  growing  VOPc  and  PbPc 
epitaxially  onto  MoSe2  substrates.  Also,  he  has  successfully  grown 
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a  wide  range  of  mdc's  and  alkali  halides  onto  each  other. 

His  primary  analysis  tool  is  RHEED,  with  some  x-ray  diffraction  and 
STM.  To  achieve  epitaxy,  extraordinarily  low  growth  rates  of  from 
0.01  -  O.lA/s  are  required  on  low  temperature  substrates.  No 
electrical  or  optical  characterization  of  the  films  has  been 
accomplished. 

Surprisingly,  Tanigaki  at  NEC  deposits  GaPcCl  and  Ga-Tetraphenyl 
porphorin  -  Cl  on  substrates  held  at  200°c  to  achieve  epitaxy. 
Tanigaki  has  also  shown  a  clear  dependence  of  single  crystal  growth 
on  base  pressure,  with  "perfect”  growth  only  occurring  at  3x10 
Torr . 

One  exciting  development  is  that  Koma  claims  to  have  grown  NaCl  on 
top  of  Pc's,  opening  up  the  door  to  heteroepitaxial  superlattices 
bonded  by  Van  der  Waals'  forces. 

The  DuPont  effort  is  also  quite  sizeable.  Here  the  emphasis  is 
almost  completely  on  growing  and  microscopically  analyzing  arious 
epitaxial  combinations  of  the  mdc's.  Again,  no  electiical  or 
optical  characterization  of  the  materials  has  been  done,  although 
the  DuPont  group  is  doing  some  rudimentary  calculational  analysis 
of  the  basic  reasons  behind  the  observation  of  Van  der  Waals 
epitaxy. 

The  use  group's  work  is  qualitatively  different  than  all  of  the 
other  work  presented.  Its  primary  focus  is  on  growth  of  fully 
organic  MQW  structures  and  organic/inorganic  heterojunctions  based 
on  the  polyacenes  (e.g.  perylene  and  napthalene) .  The  growth  rates 
are  high  (1-3A/S)  and  the  substrate  temperatures  are  very  low 
(90K) .  This  work  appeared  to  present  the  first  exciton  size 
measurements  using  organic  MQWs,  and  also  concerned  itself  with 
predicting  very  large  non  linear  optical  effects  in  crystalline 
organic  MQW  structures.  The  work  was  clearly  distinct  from  the 
Japanese  and  DuPont  research  in  that  it  subjected  the  samples  to 
both  optical  and  electrical  analysis,  and  overall  was  somewhat  more 
applications-oriented  than  the  other  work  presented.  Another 
difference  was  that  the  structural  techniques  emphasized  use  of  x- 
ray  and  optical  information  as  opposed  to  RHEED  analysis  of  the 
grown  structures. 

In  addition  to  the  papers,  I  learned  a  considerable  amount  about 
the  OMBE  efforts  both  at  Riken  in  particular,  and  worldwide  in 
general.  The  Riken  OMBE  group  involves  approximately  10  PhD  level 
researchers.  Two  or  three  actually  work  on  materials  growth,  while 
the  balance  work  on  materials  characterization.  Here,  particular 
emphasis  is  on  measuring  the  NLO  properties  of  the  films  using 
lasers  operating  over  a  wide  range  of  wavelengths  (YAG,  HeNe,  Ar, 
COj,  etc.).  Also,  a  femtosecond  laser  setup  is  in  operation  using 
a  colliding  beam  technique.  Thus  far,  the  optics  group  is  awaiting 
samples  from  the  OMBE  group,  whose  films  are  still  too  thin  to  be 
useful. 


Furthermore,  Dr.  Hara  at  Riken  indicates  that  every  major 
electronics  company  in  Japan  has  either  already  installed  or  has 
purchased  an  OMBE  system.  These  companies  include  Hitachi,  Sony, 
NEC,  KDD,  Mitsubishi,  Sharp,  Sanyo,  etc.  The  common  system 
supplier  is  Vieetech  —  a  Japanese  MBE  company  associated  with 
Thermionics  in  the  U.S.  During  FV  1990  alone,  Vieetech  took  13 
orders  from  Japanese  companies  for  OMBD  systems!  According  to 
Hara,  orders  were  still  coming  in,  but  were  being  refused  since 
Vieetech  was  unable  to  keep  up  with  production  requirements.  All 
systems  ordered  used  the  somewhat  flawed  but  nevertheless 
functional  Riken  design. 

According  to  Dr.  Hara,  each  OMBE  industrial  group  consisted  of 
approximately  5  researchers  at  the  Ph.D.  (or  equivalent)  level. 
Discussions  with  Dr.  Koma  indicated  that  a  similar  number  were 
involved  in  his  group.  At  least  two  other  professors  from  Japanese 
educational  institutions  are  involved  in  OMBE  since  they  approached 
me  at  the  conference.  If  we  presume  the  number  of  members  of  these 
latter  groups  is  approximately  three,  then  the  total  number  of 
Japanese  researchers  now  getting  involved  in  OMBE  is  estimated  at: 


75  Industrial 

10  Riken 

11  University 


96  Total 

This  compares  with  -  5-10  in  Germany  (a  new  group  at  Tubingen  is 
complementing  the  Stuttgart  work)  and  -  6-7  in  the  U.S.  (4  at  USC 
and  ~  2  at  University  of  Arizona) . 

The  focus  of  this  massive  Japanese  effort  is  still  unclear.  For 
the  time  being,  most  appear  to  be  interested  in  following  Riken 's 
lead  to  "see  what  the  future  holds".  All  are  aware  of  the  vast 
potential  of  Van  der  Waals  epitaxy  in  being  able  to  engineer 
materials  without  the  constraints  of  lattice  match. 

Organic  STM 

Riken  (Dr.  Hara)  has  taken  a  lead  in  recognizing  the  importance  of 
STM  in  understanding  OMBE  films.  Their  work  at  elucidating  liquid 
crystal  structures  using  this  technique  is  clear  from  the  paper 
attached.  Hara  has  4  STMs  in  his  laboratory  and  1  atomic  force 
microscope.  One  of  his  STMs  is  connected  through  a  UHV  mobile 
"trolley"  to  his  OMBE  system.  The  STM  connected  in  this  way  is 
also  mounted  in  a  UHV  system,  although  it  has  not  been  successfully 
used  at  this  time. 


Typically,  layers,  analyzed  were  -  1-3  monolayers  thick,  and  were 
grown  on  graphite  or  MoSj  substrates.  Effects  on  image  quality  due 
to  tip  composition,  tip  geometry  and  observational  parameters  (e,g. 
voltage  and  current)  were  all  extensively  discussed. 

This  part  of  the  conference  was  important  to  the  USC  work  in  that 
we  are  currently  attempting  to  interpret  PTCDA  images  generated 
using  the  Digital  Instruments  Nanoscope  II  STM.  Indeed  a  large 
proportion  of  the  data  presented  at  the  conference  were  obtained 
using  this  instrument. 

Perhaps  the  most  beautiful  STM  work  was  presented  by  Dr.  Heckl  at 
the  University  of  Munich.  He  presented  extensive  studies  of  DNA, 
liquid  crystals  and  other  organic  molecules. 

It  appears  that  the  work  of  the  Riken  group  is  making  an  important 
and  long  lasting  linkage  between  STM  and  OMBE.  They  are  interested 
in  this  linkage  both  for  its  application  to  resolving  structural 
detail,  and  also  for  its  use  in  "machining"  microstructures 
consisting  of  organic  molecules  manipulated  by  the  STM  tip. 

In  summary,  I  feel  it  is  of  utmost  importance  to  maintain  contact 
between  the  Riken  group  (and  hence  the  broader  Japanese  research 
establishment)  and  our  group  at  USC.  To  date,  we  have 
independently  followed  somewhat  different  "evolutionary  trends"  in 
OMBE,  both  of  which  have  their  separate  strengths  and  weaknesses. 
By  maintaining  contact,  we  can  learn  a  considerable  amount  from  the 
Japanese  effort.  To  this  end,  I  have  invited  Dr.  Kara  to  visit  USC 
during  his  visit  to  the  Digital  Instruments  STM  Users  meeting  in 
June.  Further,  I  hope  to  attend  the  NLO  meeting  at  the  University 
of  Pennsylvania  (co-sponsored  by  Riken  and  University  of 
Pennsylvania)  in  April.  Finally,  I  have  invited  Dr.  Paul  Burrows, 
a  Riken  post-doc  from  the  U.K.,  to  visit  and  work  in  our  laboratory 
once  he  finishes  his  "tour  of  duty"  at  Riken  in  the  summer  of  1991. 

It  was  a  very  valuable  experience  for  me  to  attend  this  forum,  and 
hope  to  do  so  again  at  this  same  time  next  year. 
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The  purpose  of  this  "Frontier  Forum"  is  to  bring  together  recent 
advances  in  " Nanoscopic"  science  and  technology  for  fabrication  and 
fundamental  study  of  low-dimensional  materials.  The  topics  will  focus  on 
novel  application  of  Molecular  Beam  Epitaxy  (MBE)  and  Scanning 
Tunneling  Microscopy  (STM)  techniques,  especially  for  organic 
molecular  systems.  The  official  language  is  English. 
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Van  der  Waals  Epitaxy  and  Organic  MBE 
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Heteroepitaxial  growth  of  ultrathin  films  of  organic  materials  has  potential  appli¬ 
cations  in  such  new  fields  as  nonlinear  optics,  molecular  electronics  and  molecular 
biology.  But  there  have  been  no  good  epitaxial  growth  method  which  is  applicable 
to  organic  materials.  Here  we  will  introduce  new  concepts  in  heteroepitaxial 
growth  of  organic  material  films. 

Usually  there  appear  dangling  bonds  on  a  clean  surface  of  a  substrate.  This 
makes  it  difficult  to  grow  good  heteroepitaxial  film  without  good  lattice  match  in  the 
constituent  materials,  resulting  in  very  limited  combinations  of  substrate  and  over¬ 
grown  materials.  There  are,  however,  materials  having  no  dangling  bonds  on  their 
clean  surfaces,  on  which  epitaxial  growth  proceeds  via  van  der  Waals  forces.  We 
call  that  kind  of  epitaxy  van  der  Waals  epitaxy  [1].  Good  heterostructures  can  be 
grown  even  between  materials  having  large  lattice  mismatch  in  the  van  der  Waals 
epitaxy.  We  have  proven  almost  all  layered  materials  can  be  heteroepitaxially 
grown  on  another  layered  materials  regardless  of  their  lattice  matching  conditions. 
This  idea  is  especially  useful  in  heteroepitaxial  growth  of  organic  materials,  since 
many  organic  materials  crystallize  in  forms  of  van  der  Waals  crystals.  Actually  we 
have  demonstrated  ordered  films  of  various  metal  phthalocyanines  and  coronene  can 
be  grown  on  a  cleaved  face  of  M0S2,  which  has  typical  van  der  Waals  nature.  This 
method  seems  to  be  applicable  to  many  other  organic  materials,  although  the  sub¬ 
strate  should  be  cooled  for  tlie  materials  with  high  vapor  pressures. 

In  some  cases  electrostatic  interactions  help  the  epitaxial  growth.  This  is  espe¬ 
cially  true  in  the  preparation  of  films  of  polar  organic  materials.  For  example  good 
VO-phthalocyanine  film  can  be  grown  on  various  alkali  halide  substrates.  Here 
electrostatic  interactions  witli  alkaline  or  halogen  ions  help  to  orient  polar  VO-ph¬ 
thalocyanine  molecules  on  alkali  halide  substrates  [2j. 


[1]  A.  Koma,  K.  Sunouchi  and  T.  Miyajima,  J.  Vac.  Sci.  Technol.,  B3  (1985)  724. 

[2]  H.  Tada,  K.  Saiki  and  A.  Koma,  Jpn.  J.  Appl.  Phys.,  30  (1991)  L306. 
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structural  Analysis  of  Ultrathin  Films  of  Metal-Phthalocyanines 
by  Reflection  High  Energy  Electron  Diffraction 

H.  TADA,  M.  SAKURAI.  T.  MORIOKA  and  A.  KOMA 
Dept,  of  Clieni.,  Fac.  of  Sci.,  Univ.  of  Tokyo 

Ultrathin  films  of  some  metal-phtlialocyanines  (MPc,  M=VO,  Pb,  AlCl,  and 
Cu,  Pc=C32Hi6N8)  have  been  grown  heteroepitaxially  on  cleaved  faces  of  alkali  hal¬ 
ides  (KBr,  KCl,  and  NaCl)  and  M0S2  by  molecular  beam  epitaxy  (MBE).  Their 
structures  were  investigated  in  situ  by  reflection  high  energy  electron  diffraction 
(RHEED).  Figure  1  shows  tlie  RHEED  pattern  of  VOPc  ultratliin  film  grown  on 
KBr.  The  incident  beam  is  parallel  to  the  [luO]  direction  of  the  substrate.  Sharp 
streaks  originating  from  ordered  VOPc  film  are  observed.  The  interval  is  one  third 
of  tliat  between  substrate  streaks.  Taking  into  account  the  RHEED  patterns  observed 
from  various  direcdons,  VOPc  molecules  are  found  to  form  a  single  domain  square 
lattice  as  shown  in  Fig.  2.  On  the  other  hand,  tlicy  form  two  equivalent  square  lat¬ 
tices  on  KCl  and  NaCl.  From  tJiese  results  tliat  VOPc  molecules  are  arranged  differ¬ 
ently  depending  upon  the  lattice  constants  of  alkali  halides,  it  is  considered  that  the 
arrangement  is  primarily  governed  by  the  electrostatic  interaction  between  VOPc 
molecule  and  the  substrate.  Van  der  Waals  interaction  between  neighboring  mole¬ 
cules  seem  to  play  a  role  to  determine  molecular  distances  and  arrangement  direc¬ 
tions.  Ultrathin  films  of  PbPc  and  AJClPc  molecules  grown  on  alkali  halides  gave 
almost  tlie  same  RHEED  patterns,  whereas  CuPc  films  gave  ordered  dot  patterns  due 
to  the  formation  of  a-type  crystalline  islands.  Thus  metal-phllialocyanines  with  py¬ 
ramidal  molecular  structures  seem  to  be  preferable  to  form  smooth  heteroepitaxial 
films  on  alkali  halide  substrates.  They  are  apt  to  line  up  along  [210]  direction  of  the 
substrate,  and  tliey  are  considered  to  be  stable  when  the  molecular  distances  are 
about  14  A. 

CuPc  ultrathin  film  grown  on  M0S2  gave  streaky  RHEED  patterns,  and  tlien 
CuPc  molecules  are  found  to  form  square  lattices  with  molecular  distance  of  13.7  A. 
Most  organic  ultratliin  films  prepared  on  M0S2  so  far  , (MPc,  coronene,  perylene, 
MNA)  gave  streaky  RHEED  patterns.  The  molecules-  seem  to  be  closely  packed  in 
the  first  layer  on  M0S2. 


Fig.l.  RHEED  pattern  of  VOPc  ultratliin  Fig.2.  Possible  arrangement  of  VOPc  on 
film  grown  on  KJBr.  KBr :  (•)  VO^*,  (o)  K*.  (O)  Br ' . 
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Thin  Film  Crystals  of  Phthalocyanines  and  Porphyrins 

Katsumi  Tanipaki.  Sadanori  Kiiroshima  and  Thomas  IV.  Ebbcsen 
Fundamental  Research  Laboratories,  NEC  Corporation, 

34  Miyukigaoka  Tsukuba,  Ibaraki  305,  Japan 

TeI;81-298-50-1138.  Fax:81-298-56-2326 

Phthalocyanine  and  porphyrin  families  are  promising  molecules 
for  electronics.  Especially  thin  film  crystals  of  these  molecules  have 
possibilities  to  open  the  way  for  their  advanced  applications  as 
nanoscopic  and  low-dimensional  materials  [1].  Recent  progress  in 
the  technology  of  thin  film  crystal  growth  by  molecular  beam 
epitaxy  (MBE)  is  now  being  applied  to  organic  molecules. 

We  have  applied  MBE  technique  for  thin  film  crystal  growth 
focusing  on  phthalocyanines  and  porphyrins  containing  Ill-V 
elements  in  the  periodic  table.  Tliis  is  because  d2sp3/sp3d2  hybrid 
orbitals  formed  with  tetragonally-positioned  four  nitrogen  atoms  and 
a  central  element  are  thought  to  play  an  important  role  for  tlie 
artificial  molecular  assembly  in  the  direction  perpendicular  to  the 
molecular  plane.  In  this  presentation,  we  report  thin  film  crystals  of 
a  phthalocyanine  (Pc)  and  a  mesotetraphenylporphyrin  (TPPr) 
containing  gallium  (typical  III  element)  as  a  central  element  and 
chlorine  as  a  ligand,  monochloro(phthalocyaninato)gallium  (GaPcCl) 
and  monochloro(mesotetraphenylporphyrinato)gallium  (GaTPPrCl). 

Precise  control  of  the  growth  condition  by  organic  MBE  gives 
thin  film  single  crystals  of  GaPcCl  and  GaTPPrCl  on  a  KBr  substrate 
over  2mmx2mm  area.  TEM  and  X-ray  analyses  of  these  thin  film 
crystals  reveal  tliat  their  structure  is  body-centered  tetragonal  (BCT) 
with  lattice  parameters  of  a=9.9  A  and  c=13.4  A.  Crystal  analysis 
data  of  a  bulk  single  crystal  of  GaTPPrCl,  grown  with  a  vapor-phase 
apparatus  support  the  determined  BCT  structure,  being  tetragonal, 
space  group  I4/m,  with  a=13.5  A,  c=9.9  A  for  Z=2.  Tlie  GaTPPrCl 
thin  film  density  (p=1.324)  measured  by  buoyancy  is  also  in  good 
agreement  with  the  calculated  one  from  the  BCT  structure.  Direct 
TEM  lattice  image  of  GaTPPrCl  is  demonstrated  with  a  computer 
simulation.  These  results  indicate  that  Organic  MBE  will  be  useful 
in  controlling  an  architecture  of  molecular  assembly  and  should  open 
a  new  field  from  the  leclinological  point  of  view.  .• 

Reference 

[1]  C.  W.  Dirk,  T.  Inabc,  K.  F.  Schoch,  and  T.  J.  Marks,  J.  Am.  Cliem.  Soc., 
105,1539  (1983). 
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A  New  Concept  of  Molecular  Electronics:  From  Inchworm 
Memory  To  Organic  Quantum  Well. 
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Frontier  Research  Program,  Rikcn  Institute,  Hirosawa  2-1,  Wako-shi,  Japan. 

E.G.  Wilson. 

Queen  Mary  &.  WestHeld  College,  Mile  End  Road,  London  El  4NS,  England. 

If  the  present  trend  continues,  tlie  size  of  electronic  devices  will 
decrease  to  the  size  of  organic  molecules  by  early  next  centur}'.!  This  has 
prompted  a  plethora  of  device  concepts  utilising  single  molecules  or 
polymer  chains,  but  tliere  has  been  little  experimental  evidence  that 
manipulation  and  external  contacting  of  such  devices  is  feasible,  or  that  the 
reliability  of  single  molecule  devices  could  be  made  sufficiently  high. 

A  new  device  concept  is  presented,  the  inchworm  memory^^  which 
stores  information  as  charge  on  a  three-dimensional  stack  of  ultra-thin, 
high  electron  affinity  sheets  of  conjugated  organic  rnolecules  alternated 
with  low  electron  affinity  sheets  of  "insulator"  molecules.  This  device 
offers  information  storage  on  a  molecular  scale  without  the  need  for 
lithography  at  molecular  dimensions.  A  possible  design  for  such  a  meraor}', 
utilising  thin  films  of  phthalocyanine  molecules,  is  presented. 

Organic  Molecular  Beam  Epitaxy  (OMBE)  is  a  possible  route  to  the 
production  of  such  devices.  Results  are  presented  from  raetallo- 
phthalocyanine  layers,  grown  heteroepitaxially  on  graphite  under  a  base 
pressure  of  2x10'*°  torr  by  the  OMBE  technique  with  in  siru  RHEED 
monitoring.  The  orientational  order  of  the  initial  stages  was  clearly 
determined  by  STM  even  in  air.  The  combination  of  surface-sensitive 
diffraction  methods  with  real  space  imaging  by  STM  allows  direct  studies 
of  epitaxial  layers  at  the  atomic  scale.  Preliminary  evidence  of  organic 
RHEED  intensity  oscillations  during  film  growth  is  presented. 

*  F.L.  Carter  (Ed.),  Molecular  Electronic  Devices  Marcel  Dekkcr,  New  York 

(1982,1987). 

2  P.E.  Burrows  &  E.G.  Wilson,  J.  Mol.  Elec.  209  (1990). 
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Quasi'Epitaxial  Growth  of  Organic  Multiple  Quantum  Well 
Structures  by  Organic  Molecular  Beam  Deposition 
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Abstract 


We  report  on  the  gro\A4h  and  characterization  of  thin  films  and  multiple  quantum 
wells  consisting  of  alternating  layers  of  crystalline  organic  semiconductors.  We  have 
found  that  two  compounds,  namely  PTCDA  (a  "narrow  band-gap"  perylene-based 
material)  and  NTCDA  (a  "wide  band  gap"  naphthalene  derivative)  can  be  deposited  in 
crystalline  form  on  a  wide  variety  of  substrates  using  the  UHV  process  of  organic 
molecular  beam  deposition  (OMBD).  Furthermore,  stacks  of  alternating  layers  of 
PTCDA  and  NTCDA  can  be  grown  in  a  completely  ordered  manner,  thereby  forming 
the  first  crystalline  organic  MOW  structures. 

The  optical  properties  of  the  thin  films  and  the  MQWs  have  both  been  studied. 
Measurements  of  the  PTCDA  exciton  spectrum  suggest  that  the  charge-transfer  (CT) 
exciton  can  be  quantum  confined  in  PTCDA/NTCDA  MQWs.  Solution  of 
Schrodinger's  Equation  for  an  exciton  confined  to  a  quantum  well,  both  the  blue  shift 
of  the  exciton  spectrum,  and  the  decrease  in  its  radiative  lifetime  with  decreasing 
PTCDA  layer  thickness  (from  10A  -  200  A)^can  be  explained  assuming  that  the  CT 
exciton  is  "Wannier-like"  with  a  radius  of  12  A.  Furthermore,  solutions  of  the  Davydov 
Hamiltonian  for  these  MQWs  predict  large  nonlinear  optical  effects  in  the  PTCDA 
absorption  coefficient  resulting  from  the  large  exciton-phonon  coupling  characteristic 
of  organic  molecular  semiconductors.  This  opens  up  the  possibility  that  organic 
MQWs  can  be  useful  as  light  modulators,  or  for  energy  transfer  in  two-wave  mixing 
experiments. 

Due  to  the  perfection  of  the  deposited  thin  films,  we  have  successfully  fabricated 
very  low  loss  optical  waveguides  to  TE  mode  propagation.  On  the  other  hand,  due  to 
giant  anisotropies  discovered  in  the  dielectric  properties  of  these  films  (i.e.  the 
birefringence  of  PTCDA  is  0.66).  TM  modes  are  not  guided,  leading  to  naturally 
polarizing  thin  film  waveguide  devices. 

This  work  was  supported  in  part  by  the  Air  Force  Office  of  Scientific  Research, 
Rome  Air  Development  Center  and  the  3M  Corp. 
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TUNNELING  MICROSCOPIES 
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A  field  emission  microscope  (FEM),  a  field  ion  microscope 
(FIM)  and  a  scanning  tunneling  microscope  (STM)  utilize  electron 
tunneling  and  have  realized  the  atomically  high  resolution  imag¬ 
ing.  The  high  resolution  has  been  attained  by  confining  the 
tunneling  area  to  the  size  of  an  atom.  Unique  features  of  the 
microscopes  are; 

L  Direct  projection  of  atomic  airrangement  of  a  tip  apex  by 
FIM  and  topographic  depiction  of  individual  atoms  on  a  flat 
specimen  surface  by  STM. 

Z  Electron  spectroscopy  of  a  small  apex  area  with  FEM  (field 
emission  electron  energy  spectroscopy:  FEES)  and  individual 
surface  atoms  with  STM  (scanning  tunneling  spectroscopy:  STS). 
3.  Mass  analysis  of  individual  apex  atoms  by  an  atom-probe(A-P), 
the  combined  instrument  of  an  KiM  and  a  mass  spectrometer. 
Although  the  performance  of  the  microscopes  appears  peerless, 
there  are  many  restrictions  such  as: 

L  Extremely  small,  specimens  of  FEM  and  FIM:  Shperical  apex  of 
a  sharp  tip. 

Z  Highly  sensitivity  of  STM  images  and  STS  spectra  to  the  apex 
profile  and  composition  of  a  scanning  tip. 

3.  Incomprehensible  tunneling  characteristics  of  polymers  and 
biomolecules. 

Unique  features  and  imposed  restrictions  of  the  microscopies 
will  be  discussed  based  on  the  experimental  results. 
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Organic  MBE  and  Organic  STM  for 
"  Nanoscopic "  Materials 
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‘Permanent  Address:  Department  of  Physics,  University  of 
Pennsylvania,  Philadelphia,  PA  19101,  USA. 

Application  of  molecular  beam  epitaxy  (MBE)  and  scanning  tunneling 
microscopy  (STM),  especially  for  organic  molecular  systems,  has  been  drawing 
our  attention  as  a  novel  approach  to  realizing  novel  material  structures  which 
exhibit  important  electronic  and  photonic  properties.  The  following  is  an  outline 
of  our  work  in  progress,  including  an  overview  of  the  new  concept  embraced  by 
"nanoscopic"  science  and  technology. 

1 .  Nanoscopic  Fa  b  ric  atio  n 

MBE  has  long  been  successful  as  an  ultrahigh  vacuum  (UHV)  deposition 
technique  for  fabricating  multilayer  thin  film  structures,  even  in  the  organic 
molecular  field.  In  an  extreme  case,  however,  heteroepitaxial  growth  of  organic 
ultrathin  films  on  single  crystal  substrates  can  be  attained  by  usual  adsorption, 
even  in  air,  like  a  single  crystal  growth  by  nature.  From  the  numerous 
observations  of  the  organic  layer  growth,  it  has  been  realized  that  the  anchoring 
sites  at  the  interface  are  the  most  important  factors  to  manipulate  the  adsorbates 
and  control  the  surface  structures.  From  this  point  of  view,  the  combination  of 
the  MBE  and/or  adsorption  techniques  with  surface  modification  by  STM  will 
provide  a  new  technique  for  "nanoscopic"  fabrication  of  low-dimensional 
materials,  which  allows  "molecular  planting"  at  the  atomic  scale. 

2.  Nanoscopic  Modification 

Since  the  single  atom  positioning  was  reported  by  Eigler  and  Schweizer, 
considerable  interest  has  centered  on  utilizing  STM  for  pinning  individual 
molecules.  It  has  been  well-known,  however,  that  strict  conditions  such  as  low- 
temperature  STM  under  UHV  is  required  to  obtain  reproducibility.  It  is  still 
more  realistic,  to  some  extent,  to  utilize  the  STM  for  lithography  rather  than 
manipulation.  By  introducing  the  micro-focusing  teclmique,  the  STM  controlled 
sub-nanoscopic  electron  beam  will  develop  the  possibilities,  in  principle,  for  the 
ultimate  in  "nanoscopic"  modification. 

3 .  N  aiio  s  c  0  p  i  c  ’  Char  ac  t  c  rizati  0  11 

The  STM  has  opened  a  completely  new  approach  to  organic  molecules  with 
ultrahigh  resolution.  More  recently,  the  STM  technique  has  been  shown  to 
provide  outstanding  capability  not  only  in  structural  analysis  but  also  in 
visualization  of  functional  parts  such  as  chiral  centers  and  dipoles  in  individual 
molecules.  This  kind  of  "nanoscopic"  characterization  is  indeed  encouraging  to 
the  possibilities  for  the  novel  application  of  the  STM  as  an  input/output  part  in 
next-generation  molecular  devices. 


Direct  Observation  of  Anchoring  Structures  of  Liquid 
Crystal  Molecules  by  Scanning  Tunneling  Microscopy 
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Sasabctand  Y.  Yamadaf 

Hitachi  Research  Laboratory,  Hitachi  Ltd.,  Hitachi,  Ibaraki  319-12 
tFrontier  Research  Program,  The  Institute  of  Physical  and  Chemical 
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Scanning  tunneling  microscopy  (STM)  has  been  attracting  much 
interest  as  a  useful  tool  for  the  study  of  the  molecular  alignment  on 
substrate  at  a  atomic  scale.  The  anchoring  structures  of  liquid  ciy'stals, 
n-  alkylcyanobiphenyl  (mCB=8-12),  on  substrate  are  directly  observed  by 
STM  in  order  to  analyze  the  alignment  mechanism  of  liquid  crystals. 

1.  The  substrate  dependence  of  the  anchoring  structure  of  liquid 
crystal 

The  anchoring  structures  of  smectic  liquid  cry'stal  8CB  are  observed 
condensed  on  graphite  or  MoS2  (IJ.  In  contrast  to  the  bilayer  structure  on 
graphite,  the  liquid  crystal  8CB  on  M0S2  exhibits  a  periodic  monolayer 
structure  which  suggests  that  the  anchoring  structure  can  be  varied  by  the 
substrate  structure. 

2.  The  molecular  structure  dependence  of  the  anchoring 

structure  of  liquid  crystal 

The  anchoring  structures  of  smectic  liquid  crystals  8CB,  lOCB,  12CB 

are  observed  condensed  on  M0S2  [2].  For  8CB,  the  anchoring  structure  is  of 

the  periodic  monolayer  type,  while  lOCB  and  12CB  take  a  bilayer  structure. 
These  structures  on  M0S2  depend  on  the  length  of  the  alkyl  group  in 
contrast  to  those  on  graphite. 

3.  The  relationship  between  bulk  alignment  and  anchoring 

structure  of  liquid  crystal 

The  anchoring  structure  of  the  liquid  cry'stal  is  governed,  to  a  large 
extent,  by  the  balance  of  molecular-molecular  and  substrate-molecular 
interaction.  The  balance  depend  on  the  degree  of  freedom  for  the  anchoring 
of  alkyl  chains. 

^Permanent  Address:  Department  of  Physics,  University  of  Pennsylvania, 
Philadelphia,  PA  19104,  USA 

[1]  M.  Hara  et  ai:  Nature  344  (1990)  228 

[2]  Y.  Iwakabe  el  al.  :  Jpn.  J.  Appl.  Phys.  29  (1990)  L2243 
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A  Model  of  Imaging  Mechanism 
and  Its  Applications 

Wataru  Mizutani 

Electrotechnical  Laboratory,  Umezono  1-1-4.  Tsukuba  305,  Japan 

The  interpretation  of  scanning-  tunneling  microscopy  (STM) 
Images  of  adsorbed  molecules  is  still  controversial.  Recently,  a 
resonant  tunneling  model  was  proposed  to  interpret  the  experimen¬ 
tal  results  of  STM  observation  of  liquid  crystals^^.  Tunneling 
current  comparable  with  the  experiments  can  be  calculated  based 
on  the  model. 

The  following  physical  aspect  of  the  model  is  now  under  in¬ 
vestigation.  One  of  the  remarkable  effects  which  should  be 
studied  is  charging  of  molecules  induced  by  electron  tunneling. 
Provided  that  the  Individual  molecule  is  observed,  the  energy  in¬ 
crease  by  the  single  electron  tunneling  is  estimated  to  be  about 
1  eV.  On  the  other  hand,  the  life  time  of  the  extra  electron  in 
the  molecule  is  calculated  to  be  10“^®  seconds.  Therefore,  the 
product  A  E- r  ~  10“^^  is  less  than  n—lO"^'*,  which  means  that  the 
process  is  allowed  by  the  indeterminacy  principle.  However,  the 
energy  levels  of  tlie  molecule  become  uncertain  during  that  time, 
and  it  is  not  so  explicit  that  the  resonant  tunneling  should  oc¬ 
cur  under  such  an  extreme  condition.  Probably  the  quantum  many- 
particle  effect  should  be  considered  to  analyze  the  behavior  of 
the  electrons  in  the  molecule,  but  such  a  theory  applicable  to 
STM  is  not  proposed  so  far. 

The  model  is  applicable  to  other  experimental  results,  i.e., 
a  change  in  the  molecular  arrangement  of  a  liquid  crystal  due  to 
the  bias  voltage^),  and  the  observation  of  other  adsorbed 
materials.  There  remain,  however,  some  factors  incomprehensible 
and  needs  for  further  modifications,  especially  taking. into  ac¬ 
count  of  quantum  chemistry. 

1)  W.  Mizutani  et  al . :  Appl.  Phys.  Lett.  56  (1990)  1974. 

2)  W.  Mizutani  et  al . :  J . Vac . Sci .Technol . A.  April/May  (1991). 
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We  present  a  wide  variety  of  organic  and  biological  molecules  which  have 
been  imaged  with  STM  in  air. 

1.  Hetero epitaxially  Grown  Phthalo cyanine  Monolayers 

The  combination  of  diffraction  methods  with  real  apace  imaging  of  STM 
allows  direct  studies  of  epitaxial  layers  at  the  atomic  scale.  Metailophthalocyanine 
monolayers  were  grown  heteroepitaxially  on  graphite  under  2  x  10  ’°  torr  by  the 
organic  MBE  technique  with  in  situ  RHEED  monitoring.  From  STM  images  of  the 
phthalocyanine  layer  and  the  substrate  beneath  the  organic  layer,  the  orientational 
correlation  between  the  organic  layer  and  the  substrate  surface  was  clearly 
determined  in  real  space. 

2.  Novel  Anchoring:  Structures  and  Phase  Transitions  in 

Liquid  Crystals'^ 

While  STM  images  of  liquid  crystal  molecules  have  been  reported  with  high 
resolution,  only  graphite  has  been  used  as  a  substrate  so  far.  We  replaced  the 
graphite  with  M0S2  and  observed  by  STM  an  entirely  new  "2  x  4"  anchoring 
structure  for  n-alkyl  cyanobiphenyls.  Furthermore,  the  results  suggest  that  a 
possibility  exists  to  study  the  phase  transition  mechanisms  using  the  STM  imaging. 

3.  Visualization  of  Functional  Parts  in  Ferroelectric 

Liquid  Crystals*'^' 

The  STM  technique  has  been  shown  to  provide  outstanding  results  not  only  in 
structural  analysis  but  also  in  visualization  of  functional  parts  such  as  chiral  centers 
and  dipoles  in  organic  molecules.  By  carefully  controlling  deposition  parameters  in 
the  sample  preparation  step,  it  is  possible  to  distinguish  between  methyl  and  ethyl 
moiety  around  chiral  carbons  and  orientations  of  dipoles  in  individual  molecules. 

4.  Novel  Guest- Ho  St  Anchoring  Method  for  STM 

Visualhation  of  Biological  Macroinolecules'^** 

Serious  problems  have  been  encountered  in  trying  to  obtain  unambiguous 
STM  images  of  biological  macromolecules,  mainly  because  of  lack  of  reliable 
deposition  methods  on  atomically  flat  substrates.  The  molecules  must  not  be  packed 
into  dense  aggregates,  but  must  be  distributed  in  a  two-dimensional  plane,  i.e.  sub¬ 
monolayer.  Thus  during  deposition,  high  molecular  mobility  on  the  surface  is 
desirable,  while  moderate  anchoring  is  required  during  scanning.  Here  we 
introduce  a  novel  and  easy  deposition  method,  ’’Guest-Host  Anchoring"  of 
biological  macromolecules  with  clear  discrimination  against  artifacts. 
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Surface  phases  of  small  organic  molecules 
observed  by  tunneling  microscopy 


W.M.  Heckl  *■  G.  Binnig  D.P.E.  Smilh+  and  H.  Klaggcs  * 

*  University  Munich,  Schcllingstr.  4.  D-8000  Miinchcn  40 
+  IBM  Research,  Physics  Group  Munich 

The  two  dimensional  order  formed  by  epitaxial  growth  of  small 
organic  molecules  on  the  surface  of  graphite  and  molybdenum  disulfide 
can  be  investigated  with  nearly  atomic  resolution  by  scanning  tunneling 
microscopy  (STM).  We  will  discuss  the  structure,  the  registry,  the 
adsorption  and  imaging  mechanism  as  well  as  lattice  defects  and 
dynamical  aspects  of  the  interfacial  layer.  We  will  concentrate  on  DNA 
bases  and  show  some  additional  examples  of  amphiphilic  molecules  and 
hexadecylpyrene  as  a  fluoropiiore. 

Several  recent  studies  have  reported  on  the  possibility  of  imaging 
DNA  by  Scanning  Tunneling  Microscopy.  However  tJie  ultimate  goal,  that 
is  to  read  the  code  of  the  strands,  has  not  yet  been  achieved.  In  order  to 
do  so  it  must  be  possible  to  clearly  recognize  and  distinguish  between  the 
four  different  bases  which  comprise  the  genetic  code.  One  essential 
requirement  for  imaging  small  organic  molecules  by  STM  is  to  find 
experimental  preparation  conditions  where  the  molecules  stick  firmly  to 
the  substrate  and  form  a  layer  with  great  stability  in  order  to  withstand 
the  forces  of  tJie  STM  tip  during  imaging.  Compared  to  the  binding  of 
complete  DNA  strands  to  the  basal  planes  of  M0S2  or  graphite  the 
adsorption  of  "naked"  nucleotide  bases  is  favoured  by  their  greater 
hydrophobicity  and  by  their  ability  to  register  with  the  substrate  fonning 
a  two  dimensional  ordered  array.  We  found  that  the  two-dimensional 
lattice  formed  by  adenine,  guanine  and  thymine,  three  of  the  four  bases 
which  comprise  the  genetic  code  of  DNA,  is  nearly  identical  on  the 
surface  of  the  graphite  and  the  molybdenum  difsulfide  substrate.  The 
resolution  of  molecular  details  is  superior  in  the  case  of  the  graphite 
substrate.  By  the  aid  of  molecular  orbital  calculations,  which  allow  for 
the  comparison  of  predicted  electron  densities  with  the  observed  STM 
contrast,  we  propose  the  interfacial  structure  of  guanine  to  be  composed 
of  linear  strands  of  guanine  molecules  with  alternating  orientation  linked 
together  by  hydrogen  bonds. 

Lit.  W.M.  Heckl  ct  al.:  Nature,  submitted 

D.P.E.  Smith  and  W.M.  Heckl,  Nature,  Vol.  346,  No.  6285,  616  -  617,  1990 
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o  n 

Surface  Chemistry  and  Physics  of  Low  Dimensional  Layered 
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Fumio  S.  Ohuchi 

E.l. DuPont  de  Nemours  and  Company 
Central  Research  and  Development  Department 
Experimental  Station 
Wilmington,  Delaware  19880-0356 

Abstract 

The  potential  for  new  materials  with  tailored  bulk,  surface  and  interface 
properties  has  motivated  many  recent  investigations  concerned  with  fabricating 
composite  systems  comprised  of  dissimilar  materials.  We  have  been 
developing  processes  involved  in  layered  materials  through  tightly  controlled 
atomic  and  electronic  structure  at  'he  interfaces. 

A  route  to  nev/  materials  by  exploiting  atomically  abrupt  nature  of  layered 
compound  materials  and  their  reactivity  to  metals  is  described  first,  then  a 
variation  of  molecular  beam  epitaxy,  called  van  der  Waals  epitaxy,  will  be 
followed.  Both  methods  have  utilized  unique  natures  of  low  dimensionality  of 
the  materials,  and  are  rich  in  surface  physics  and  chemistry. 

The  interactions  of  group  IB  metals  (Cu,  Ag  and  Au)  with  van  der  Waals 
surfaces  of  layered  semiconducting  materials  have  led  new  classes  of  material 
systems  that  have  not  been  obtained  by  other  methods.  The  peculiarity  of 
interaction  lies  in  the  fact  that  the  host  structure  is  forced  into  a  lattice  distortion 
by  electronic  driving  force  from  surface  reaction  and  intercalation.  Difference  in 
the  surface  symmetry  also  influences  the  nucleation  and  growth  of  the  metal 
atoms,  which  in  turn  demonstrates  the  decoration  of  atoms  over  the  surface. 

Lattice  mismatching  difficulties  that  limit  the  choice  of  materials  in  conventional 
MBE  of  three  dimensional  systems  are  circumvented  in  the  van  der  Waals 
epitaxial  process.  Demonstrated  are  epitaxial  growth  of  layered  materials  with 
lattice  mismatching  exceeding  more  than  +18%.  STM  provided  real  space 
images  of  the  morphology  of  the  epitaxial  layer  and  showed  novel  structures 
resulting  from' the  large  lattice  mismatch  where  the  epilayer  atoms  are 
commensulated.  The  flexibility  for  fabricating  multilayer  structures  containing 
different  materials  selected  for  their  specific  properties,  rather  than  their  lattice 
matching  requirement,  has  been  discussed. 
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Appendix  C:  Polarization  effects  in  organic  MQWs 


Consider  an  exciton  with  a  Bohr  radius  a©  generated  at  point  P  in  a 
PTCDA  layer  of  thickness  t  sandwiched  between  two  NTCDA  layers  with  the 
same  thickness,  as  shown  in  Fig.  18.  The  change  in  exci+on  potential  due  to  the 
presence  of  an  NTCDA  molecule  at  point  N  on  the  NTCDA  surface  is; 

A(p  =  A|iCOS0'/eplR-r|2  (A-1) 

where  Ap  =  pp  -  pN,  with  pp  (pn)  the  dipole  moment  per  unit  volume  of  PTCDA 
(NTCDA),  Ep  is  the  dielectric  constant  of  PTCDA,  R  is  a  vector  between  points  P 
and  N,  0'  is  the  angle  between  R  and  the  normal  to  the  NTCDA  plane,  a  :d  the 
charge  density  of  the  exciton  is  described  by  the  spherical  coordinates:  r,  0, 
and  p  .  Using  Taylor's  series  expansion  and  dropping  the  higher  order  terms  in 
Eq.  A-1,  the  change  in  potential  is: 

A(|j  =  Apcos0'/(epR2)[1  +  (2r/R)  cos0cos0’]  (A-2) 

Here  we  assume  the  ground  state  exciton  wavefunction  to  be  a  hydrogenic  Is 
state,  and  the  first  term  to  be  proportional  to  the  first  order  perturbation 
Hamiltonian.  Therefore,  the  first  order  polarization  energy  correction  is: 

E(1)  =  <1s  I  H(1)  I  1s>  =  qApcos0V(ep  R2)  (A-S) 

On  the  other  hand,  assuming  \|/n  is  the  excited  state  wavefunction  and  H(2)  is  the 
second  term  in  Eq.  A-2,  then  the  second  order  energy  is: 

E(2)  =  In  I  <n|H(2)l  1  s>  I  2  /  (En  -  Eis)  (A-4) 

Here,  n  denotes  the  nth  excited  state  wavefunction,  and  En  and  Eis  af©  fhe 
excited  and  ground  state  energies,  respectively. 

The  first  order  polarization  energy  can  be  evaluated  by  integrating  EC)  in 
Eq.  A-3  over  all  space  to  obtain: 

<E(i)tot>  =  A' (Ap/ep)t  {A-5) 

where  A'  is  a  numerical  constant.  Similarly,  ignoring  the  higher  order  term,  the 
second  order  polarization  energy  can  be  evaluated: 

<E(2)tot>  =  A"Ap2ao25/(ept4)  (A-6) 

where  A"  is  also  a  numerical  constant,  and  5  is  the  thickness  of  the 
PTCDA/NTCDA  heterointerface  region.  Therefore,  the  first  order  energy 
correction  to  the  polarization  energy  due  to  the  presence  of  the  NTCDA  layers 
increases  linearly  with  t,  and  the  second  order  term  is  proportional  to  t''^- 
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Table  1  -  van  der  Waals  Coefficients(^) 


Bond  Interaction 

a  (meV-A6) 

p  (meV) 

1 

•< 

H-H 

1.96 

432 

4.54 

C-H 

5.51 

3735 

4.57 

C-C 

15.7 

39400 

4.59 

0-0 

9.41 

5860 

4.59 

c=o-o=c 

16.1 

9110 

4.59 

H-0 

3.92 

1160 

4.57 

H-0=C 

5.29 

1500 

4.57 

C-0 

12.1 

15200 

4.59 

c-o=c 

15.9 

19000 

4.59 

C)Values  for  a,  (3  and  y  taken  from  A.  Abe,  R.  L  Jernigan  and  P.J.  Flory,  J.  Am. 
Chem.  Soc.,  fiS.  631  (1966);  and  R.  A.  Scott  and  H.  A.  Scheraga,  J.  Chem. 
Phys.,42,  2209  (1965). 


Table  2.  Variational  calculation  parameters 


Parameter 

Symbol 

Unit 

Value 

Hole  mass 

mh 

m© 

0.18 

Election  mas 

me 

mo 

10.0 

Electron  potential 

Ve 

meV 

900 

Hole  potential 

Vh 

meV 

50 

Relative  dielectric 
constant  of  PTCDA 

3.6 
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Table  3 :  Reflectivity  data  for  PTCDA  thin  films 


X(|Am) 

Sample 

Thickness 

Incident 

Angle* 

nj_  (Cal.) 

1.064 

1.13 

290 

1.354  (m=8) 

1.319 

1.25 

470 

1.364  (m=8) 

1.064 

1.13 

410 

1.343  (m=6) 

1.319 

1.25 

55.50 

1.346  (m=6) 

*  Corresponding  to  the  reflectivity  minimum. 


Table  4:  Resonant  Nonlinear  Properties  of  Some  Typical  Materials 


Material 

Energy 
Gap  (eV) 

>.(pm) 

Response 

Time 

Temp. 

(K) 

n2 

(cm2A/V) 

Refs. 

InSb 

0.18 

5.5 

0.3  ps 

80 

10-3 

[29] 

In  As 

0.40 

3.1 

0.2  ps 

80 

2x10-5 

[30] 

Si 

1.12 

1.06 

— 

300 

1.2x10-''0 

[31] 

GaAs 

1.42 

0.81 

30  ps 

80 

o.exio-io 

[32] 

SINC* 

0.81 

— 

300 

10-7 

[33] 

PTCDA 

2.2 

0.514 

16ps 

295 

5.4x10-5 

This  work 

*  SINC  is  silicon  naphthalocyanine  oligomer,  which  is  a  random  glassy  polymer. 


Figure  Captions 

Fig.  1:  Two  views  of  a  unit  cell  of  PTCDA. 

Fig.  2:  Two  views  of  a  unit  cell  of  NTCDA. 

Fig.  3:  Plot  of  contours  of  constant  van  der  Waals  energy  between  a  PTCDA 
and  NTCDA  molecule.  Each  contour  corresponds  to  an  increase  in  energy  of  5 
meV  as  one  moves  from  the  central  position  (0,0)  to  the  edges.  Each  axis 
spans  a  distance  of  ±  2k  from  the  position  where  an  NTCDA  molecule  is 
centered  directly  over  the  center  of  a  PTCDA  molecule,  with  an  interplanar 
stacking  distance  of  3.2  A. 

Fig.  4:  Bipolar  current  voltage  characteristics  of  a  CuPc/PTCDA  hetero¬ 
junction  sample.  Inset:  Cross-sectional  view  of  the  CuPc/PTCDA  test 
structure. 

Fig.  5:  Saturation  current  versus  temperature  for  a  CuPc/PTCDA 
heterojunction. 

Fig.  6:  Proposed  band  diagram  of  a  CuPc/PTCDA  heterojunction  in  the 
region  of  the  heterointerface. 

Fig,  7:  Proposed  molecular  level  scheme  of  a  CuPc/PTCDA  heterojunction. 
Fig.  8:  Schematic  view  of  a  PTCDA/NTCDA  multiple  quantum  well  structure. 
Fig.  9:  Room  temperature  absorption  and  20K  luminescence  spectra  for  a  5 
period,  40  A  symmetric  PTCDA/NTCDA  MOW.  The  luminescence  spectrum  is 
obtained  by  pumping  the  MOW  sample  using  554  nm  light  from  an  Ar  ion  laser. 
Fig.  10:  Singlet  exciton  energy  as  a  function  of  PTCDA  well  width  for 
PTCD/VNTCDA  MQWs.  Fits  to  the  data  are  obtained  using  parameters  given  in 
Table  2,  and  are  described  in  the  text.  Also  shown  is  the  exciton  radius,  a,  as  a 
function  of  well  width. 

Fig.  11:  Capacitance  measured  versus  frequency  for  a  horizontal, 
interdigitated  capacitor  with  PTCDA  serving  as  the  dielectric. 

Inset:  Two  views  of  a  unit  cell  of  PTCDA. 

Fig.  12:  Reflectance  at  1  =  1.064  pm  versus  incident  beam  angle  for  a  1.25  pm 
thick  PTCDA  film  on  a  quartz  substrate.  Data  is  shown  as  closed  circles,  and 
theory  is  shown  as  a  solid  line. 

Inset:  Index  of  refraction  ellipsoid  for  PTCDA  showing  the  relationship 
between  the  indices  along  different  thin  film  directions. 

Fig.  13:  Transmitted  light  intensity  for  TE  and  TM  waves  in  a  10  mm  long,  2 
pm  wide  PTCDA  rib  waveguide  as  a  function  of  incident  angle  of  polarization. 
Measurements  are  made  at  I  =  1 .064  pm. 
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Inset:  Experimental  set-up  used  to  characterize  the  PTCDA  waveguide- 
polarizer  under  study. 

Fig.  14:  Temperature  dependent  intensity  measured  in  a  PTCDA  Fabry-Perot 
waveguide  resonator  at  X  =  1.06  pm  with  a  length  of  1.62  mm  and  a 
width  of  2  pm. 

Inset:  The  experimental  set-up,  1:  Argon  laser  (X  =  0.514  pm),  2:  AO 

modulator,  3:  Ring  YAG  laser  (X  =  1.06  pm),  4:  Glan-Thompson 
Polarizer,  5:  Thermoelectric  cooler  and  thermocouple,  6:  TV  Camera, 
7:  Monitor,  8:  Iris,  9:  APD,  10:  Oscilloscope. 

Fig.  15:  Intensity  of  guided  light  (X  =  1.06  pm)  versus  pump  light  intensity  (X  = 
0.514  pm)  in  a  PTCDA  waveguide  with  the  same  parameters  in  Fig.  1 . 
Fig.  16  :  Rise  and  fall  times  of  the  1.06  pm  wavelength  beam  due  to 

modulation  of  the  X=  0.514  pm  beam. 

Fig.  17:  Nonlinear  absorption  in  a  weak  probe  beam  of  several 
PTCDA/NTCDA  MOW  structures  as  a  function  of  pump  beam  intensity  (which  is 
proportional  to  the  square  of  the  Rabi  Frequency). 

Fig.  18:  Geometry  used  for  calculating  the  polarization  energy  of  a  PTCDA 
exciton  sandwiched  between  two  NTCDA  layers  separated  by  distance  t. 


71 
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Evidence  for  Exciton  Conrinement  in  Crystalline  Organic  Multiple  Quantum  Wells 

F.  F.  So  and  S.  R.  Forrest 

Center  for  Photonic  Technology.  Departments  of  Electrical  Engineering! Electrophysics  and  Materials  Science. 

University  of  Southern  California.  Los  Angeles.  California  90089-0241 
(Received  31  January  1991) 

Mulliple-quantum-well  structures  based  on  two  crystalline  organic  semiconductors,  namely,  3,4,9,10 
perylenetetracarboxylic  dianhydride  and  3.4,7,8  naphthalenetetracarboxylic  dianhydride,  have  been 
grown  by  organic  molecular-beam  deposition.  Both  optical-absorption  and  time-resolved  photolumines¬ 
cence  measurements  reveal  a  significant  effect  on  the  binding  energy  and  the  radiative  recombination 
probability  of  excitons  due  to  localization  of  carriers.  Variational  calculations  of  the  ground-state  exci¬ 
ton  energy  in  quantum  wells  have  been  done,  and  the  results  agree  with  the  experimental  data.  This  pro¬ 
vides  evidence  for  exciton  confinement  in  organic  quantum-well  structures. 


PACS  numbers:  72.80.  Le 

Ittorganic  semiconductor  muhipic-quantum-well 
(MQW)  structures  have  been  a  subject  of  intense  study 
in  the  past  decade  due  to  their  interesting  electronic  and 
optical  properties  which  have  many  promising  applica¬ 
tions  in  optoelectronic  device  technology. '  However,  ep¬ 
itaxial  growth  of  semiconductor  heterostructures  is  limit¬ 
ed  to  materials  with  small  lattice  mismatch.^  In  con¬ 
trast,  organic  crystals  arc  bonded  together  by  the  rela¬ 
tively  weak  van  dcr  Waals  forcc,^  which  allows  for  the 
layering  of  materials  of  widely  differing  crystal  lattices 
whose  strain  energy  is  insufficient  to  create  lattice  de¬ 
fects.  In  principle,  therefore,  it  is  possible  to  grow  de¬ 
fect-free  heterostructures  with  a  wide  range  of  crystal¬ 
line  organic  semiconductors.  Recently,  highly  ordered 
organic  MQW  structures  based  on  3,4,9,10  perylenet¬ 
etracarboxylic  dianhydridc  (PTCDA)  and  3,4,7,8  na- 
phthalenetetracarboxylic  dianhydride  (NTCDA)  have 
been  grown  by  the  ultrahigh-vacuum  process  of  organic 
molecular-beam  deposition  (OMBD).^’  The  so-called 
“quasiepitaxial”  growth  of  organic  thin  films  has  also 
been  demonstrated  on  a  variety  of  substrate  materials  in¬ 
cluding  inorganic  semiconductors  such  as  InP  and  Si. 
More  surprisingly,  highly  ordered  organic  thin  films  can 
be  grown  on  noncrystalline  substrates  such  as  fused 
quartz  and  glass.  Thus,  new  opportunities  for  engineered 
heterostructures  consisting  of  a  wide  range  of  organic 
material  combinations  layered  without  regard  to  lattice 
match  provide  exciting  possibilities  in  the  field  of  op¬ 
toelectronics. 

Both  the  optical  and  electronic  properties  of  organic 
semiconductors  are  governed  by  the  nature  of  excitons 
in  the  crystal.  Excitons  in  most  organic  crystals  are 
thought  to  be  Frenkel-like,^  although  there  arc  a  few  re¬ 
ports  suggesting  that  excitons  in  anthracene  crystals  are 
Wannicr*  ’  in  nature.  Therefore,  it  is  important  to  de¬ 
termine  the  characteristics  of  excitons  in  the  archetype 
compound  PTCDA  in  order  to  understand  its  optical  and 
electronic  properties.  Crystalline  organic  semiconductor 
heterojunctions  (HJs)  have  been  grown  previous  to  this 
work,  and  their  electrical  and  optical  properties  indicate 
that  the  HJ  band  offsets  control  the  flow  of  charge  be¬ 


tween  the  contacting  semiconductors  in  a  manner  similar 
to  inorganic  hetcrostructurcs.*"'®  Thus.  MQW  struc¬ 
tures  based  on  PTCDA  and  NTCDA  provide  the  ideal 
•ools  for  a  study  of  the  nature  of  excitons  in  crystalline 
organic  semiconductors. 

In  this  paper  we  report  the  measurement  and  analy¬ 
sis  of  optical-absorption  spectra  of  organic  multiple- 
quantum-well  structures  grown  by  OMBD.  An  increase 
of  the  singlet  ground-state  exciton  energy  in  PTCDA 
with  decreasing  layer  thickness  has  been  observed,  and 
the  results  can  be  understood  in  the  context  of  a  change 
in  exciton  binding  energy  due  to  quantum  confinement  in 
a  potential  well"  '*  formed  by  the  energy-band  offsets 
between  the  contacting  materials.  The  results  of  our 
variational  calculations  on  the  well-width  dependence  of 
exciton  binding  energy  are  in  good  agreement  with  ex¬ 
perimental  data.  Time-resolved  photoluminescence  mea¬ 
surements  were  also  done  for  organic  MQW  structures, 
and  the  exciton  lifetimes  were  found  to  decrease  with  de¬ 
creasing  PTCDA/NTCDA  layer  thickness.  The  de¬ 
crease  in  exciton  lifetime  is  explained  in  terms  of  a 
shrinkage  in  exciton  volume  in  the  MQW’s, "  and  is  also 
a  result  of  quantum  confinement  of  charge  carriers. 

The  organic  MQW  structures  investigated  in  this  work 
consist  of  ultrathin,  alternating  layers  of  PTCDA  and 
NTCDA  of  equal  thickness,  with  layer  thicknesses  in  the 
several  samples  ranging  from  10  to  2(X)  A.  Layer 
thicknesses  are  estimated  to  be  uniform  to  within  about 
two  molecular  layers  (i.e.,  =6-7  A).  Organic  MQW 
structures  reported  in  this  paper  were  grown  on  glass 
substrates  by  organic  molecular-beam  deposition ‘ — a 
process  analogous  to  conventional  molecular-beam  epi¬ 
taxy  used  to  grow  inorganic  semiconductor  heterostruc- 
tures.  The  key  to  achieving  good  surface  morphology 
and  sharp  interfaces  was  to  keep  the  substrates  at  about 
90  K  during  growth.  Low-temperature  quasiepitaxial 
growth  of  crystalline  organic  thin  films  on  inorganic  sub¬ 
strates  has  been  demonstrated  for  a  number  of  crystal¬ 
line  organic  materials. ^  Under  such  conditions, 
highly  ordered  molecular  stacks  are  observed  for  our 
PTCDA/NTCDA  MQWs  using  both  birefringence  and 
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FIG.  I.  Typical  room-temperature  optical-absorption  and 
low-temperature  (20  K)  photoluminescence  spectra  of  a 
PTCDA/NTCDA  mulliple-quanium-well  sample. 


x-ray  measurements  in  the  multilayer  structure  consist¬ 
ing  of  PTCDA  and  NTCDA.*  “ 

Optical-absorption  measurements  were  done  at  room 
temperature  using  a  spectrophotometer.  Time-resolved 
photoluminescence  measurements  of  organic  quantum- 
well  samples  were  done  at  20  K  using  a  multiline  cw 
Ar-ion  laser  with  an  output  of  400  mW.  At  this  energy, 
excitons  were  generated  only  in  the  PTCDA  layers  since 
the  NTCDA  layers  arc  transparent  to  the  incident  light. 
Laser  pulses  were  generated  using  an  acousto-optic 
modulator  with  a  rise  and  fall  time  of  less  than  4  ns. 
The  laser-pulse  width  was  10  ns.  The  time  decay  of  the 
monochromatic  luminescence  signal  was  measured  using 
a  double-pass  0.75-m  Spex  monochromator  in  conjunc¬ 
tion  with  a  multichannel  photon  counting  averager.  De- 
convolution  of  the  measured  luminescence  signal  from 
the  instrumental  temporal  response  at  the  shortest  times 
measured  (--5  ns)  has  been  carried  out  in  all  data 
analysis. 

The  room-temperature  optical-absorption  and  low- 
temperature  (20  K)  photoluminescence  spectra  of  a  five- 
period  symmetric  PTCDA/NTCDA  organic  MQW  sam¬ 
ple  with  an  individual  layer  thickness  of  40  A  is  shown  in 
Fig.  I.  The  characteristic  absorption  spectra  of  both 
PTCDA  and  NTCDA  are  present,  and  the  absorption 
bands  are  due  to  exciton  generation  n  the  corresponding 
crystalline  layers.  The  absorption  bands  are  relatively 
broad  because  of  the  strong  exciton-phonon  interaction 
characteristic  of  organic  crystals.^  From  the  lowest- 
energy  cutoffs  of  the  spectra,  the  “energy  gaps”  are 
determined  to  be  2.2  and  3.1  eV  for  PTCDA  and 
NTCDA,  respectively.  We  have  found  that  the  lowest- 
energy  (0-0)  singlet-exciton  absorption  peak  in  PTCDA 
(indicated  by  the  arrow)  shifts  to  higher  energy  as  the 
layer  thickness  is  decreased.  No  apparent  shifts  in  the 
higher-energy  exciton  lines  were  observed,  possibly  as  a 
result  of  the  broad  nature  of  these  lines.  The  dependence 
of  the  ground-state-exciton  energy  shift  with  well  width 
is  shown  by  the  data  points  in  Fig.  2. 

The  blueshift  in  the  exciton  energy  in  PTCDA  can  be 
due  to  molecular  polarization  in  the  presence  of  an  exci- 
lon  which  extends  into  the  NTCDA  in  the  multilayer 
structure.^  It  has  been  shown  that  the  first-order  exciton 


FIG.  2.  Exciton  line  shift  as  a  function  of  well  width  The 
solid  line  is  the  result  of  variational  calculations.  The  varia¬ 
tional  parameter  a  is  also  shown. 

energy  correction  due  to  polarization  arising  from  the 
multilayer  structures  decreases  the  total  exciton  energy 
by  a  quadratic  function  of  increasing  layer  thickness. 
However,  the  observed  dependence  on  layer  thickness  is 
much  weaker  than  that  predicted  from  the  polarization 
model. 

Alternatively,  the  exciton  energy  shift  can  be  due  to 
the  change  in  exciton  binding  energy  as  a  result  of  quan¬ 
tum  confinement.  As  the  layer  thickness  decreases,  the 
exciton  motion  becomes  two  dimensional.  The  exciton  is 
“squeezed”  in  the  potential  well  formed  in  the  PTCDA 
layer,  resulting  in  an  increase  in  exciton  binding  energy 
In  this  model,  the  small-band-gap  PTCDA  layer  is  treat¬ 
ed  as  a  potential  well  bounded  by  the  large-band-gap 
NTCDA  layers  in  a  type-I  superlatticc  configuration. 
This  picture  is  consistent  with  observations  made  for  sin¬ 
gle  PTCDA/NTCDA  heterojunctions  using  current- 
voltage  measurements  similar  to  those  used  to  determine 
the  properties  of  other  PTCDA-based  heterojunctions.* 
Thus,  electrons  and  holes  are  confined  in  PTCDA  layers 
of  width  L:.  The  total  Hamiltonian"  for  an  exciton  in  a 
quantum  well  is 

//,oi-//,  +  ///,  +  //,.  ,  (I) 

where  H,  and  //*  describe  the  motion  of  electrons  and 
holes  along  the  z  axis  in  the  quantum  well,  and  //,,  de¬ 
scribes  the  two-dimensional  motion  of  excitons  in  the  x-y 
plane  (parallel  to  the  heterointerfaces).  The  relative 
dielectric  constant  of  PTCDA  used  in  this  calculation  is 
3.6,'*  assuming  the  crystal  is  an  isotropic  dielectric  con¬ 
tinuum.  While  it  is  known  that  dielectric  anisotropies  in 
such  films  arc  large,  the  isotropic  continuum  model  nev¬ 
ertheless  has  been  found  to  provide  a  very  accurate 
description  of  charge-transfer  exciton  dynamics  in  simi¬ 
lar  organic  crystals  such  as  anthracene.'’ 

The  Schrodinger  equation  can  then  be  solved  by  the 
variational  principle  using  the  following  trial  wave  func¬ 
tion: 

'•'ioi"“^r'P*exp(-(rVa--Fz  V^*)''']  .  (2) 


2650 


Volume  66,  Numepk  20 


PHYSICAL  REVIEW  LETTERS 


20  Ma'i  19*^! 


Here,  z  “r,  -zh  and  r  —  r,  -r*,  where  r,,  r,,  r»,  are 

the  positions  for  electrons  (e)  and  holes  (h)  in  the  z 
direction  and  the  x-y  plane.  In  Eq.  (2),  the  total  wave 
function  if  tot  consists  of  the  hydrogen  U-like  wave  func¬ 
tion  ('i'l,),  modulated  by  the  harmonic  part  of  the  wave 
function,  i.e.,  the  electron  ('!',)  and  hole  (if/,)  wave 
functions.  Also,  a  is  a  variational  parameter  in  the  trial 
solution,  and  p  is  set  equal  to  ao,  where  oq  is  the  exciton 
Bohr  radius  in  the  bulk  crystal. 

It  should  be  noted  that  the  trial  wave  function  used  is 
assumed  to  be  ellipsoidal  with  p’^ao.  Such  an  ellipsoid¬ 
al  trial  wave  function  is  consistent  with  the  results  of 
the  theoretical  calculations  done  by  Shinozuka  and 
Matsuura. For  very  small  well  widths,  the  harmonic 
part  of  Eq.  (2)  dominates  the  exciton  :  dimension. 

To  solve  the  Schrodinger  equation,  the  binding  energv 
is  minimized  as  a  function  of  a.  The  results  of  the  calcu¬ 
lation  are  shown  by  the  solid  line  in  Fig.  2.  The  parame¬ 
ters  used  in  this  model  are  listed  in  Table  I.  Here,  the 
sum  of  the  electron  potential  (K, )  and  the  hole  potential 
(F/,)  is  equal  to  the  difference  in  energy  gaps  of  PTCDA 
and  NTCDA.  Of  all  the  parameters  used  in  the  calcula¬ 
tion,  the  results  are  only  sensitive  to  the  choice  of  the 
effective  hole  mass  m*.  Using  these  parameters,  a  good 
fit  to  the  experimental  data  is  obtained  for  m*  “0.1 8mo. 
where  mo  is  the  electron  rest  mass.  Note  that  a  small 
value  of  mi,  and  a  relatively  large  value  of  effective  elec¬ 
tron  mass  m,  were  chosen  in  the  calculation.  A  large 
difference  in  electron  and  hole  effective  masses  is  com¬ 
mon  in  organic  crystals  such  as  naphthalene  and  anthra¬ 
cene. Indeed,  in  previous  work,  it  has  been  shown  that 
PTCDA  is  a  preferentially  hole-transporting  material, 
indicating  that  is  substantially  less  than  m,.*  The 
effect  of  different  values  of  mi,  are  shown  for  comparison 
in  Fig,  2.  Note  that  the  fit  to  the  data  is  good  except  for 
the  10- A  sample.  This  can  have  several  causes,  includ¬ 
ing  nonuniform  layer  thickness  and  penetration  of  the 
exciton  wave  function  into  adjacent  PTCDA  layers. 

The  calculated  values  of  a,  as  shown  in  Fig.  2,  increase 
monotonically  with  increasing  well  width,  and  asymptoti¬ 
cally  approach  the  free-exciton  Bohr  radius  of  12  A,  in¬ 
dicating  that  the  exciton  is  confined  in  the  quantum 
wells.  The  dependence  of  a  on  Z,.-  for  the  ellipsoidal  ex¬ 
citon  can  be  understood  as  follows:  As  well  width  de¬ 
creases,  'F|,  becomes  a  2D  wave  function  resulting  in  a 
decrease  in  its  extent  in  the  z  direction.  At  the  same 


TABLE  1.  Parameters  used  in  the  variational  calculations. 


Parameters 

Symbol 

Unit 

Value 

Hole  mass 

mt, 

mo 

0  18 

Electron  mass 

mo 

lO.O 

Electron  potential 

y. 

meV 

900 

Hole  potential 

Relative  dielectric 

F/, 

meV 

50 

constant  of  PTCDA 

er 

3.6 

lime,  the  effective  Bohr  radius  (a)  in  the  x->  direction 
also  decreases  with  decreasing  Z..  in  order  for  the  wave 
function  to  retain  spherical  symmetry  In  the  case  of  an 
extremely  narrow  well,  the  exciton  is  squeezed  to  the  ex¬ 
tent  that  its  wave  function  is  forced  to  deviate  from 
spherical  symmetry 

To  test  this  conclusion  of  quantum  confinement,  pho¬ 
toluminescence  (PL)  measurements  were  carried  out  by 
optically  pumping  the  MQW’s  with  an  Ar-ion  laser. 
From  the  absorption  spectrum  in  Fig.  1,  it  is  apparent 
that  NTCDA  is  transparent  to  the  excitation  from  the 
Ar-ion  laser,  and  therefore  the  PL  spectrum  for  a  MQW 
sample  shown  in  the  figure  is  due  to  PTCDA  alone 
Also,  no  PL  signal  is  obtained  when  NTCDA  is  optically 
pumped  by  the  Ar-ion  laser.  From  the  PL  spectrum,  we 
can  identify  two  major  emission  bands  at  720  and  800 
nm.  and  a  shoulder  at  about  680  nm.  From  the  energy 
difference  between  the  absorption  and  luminescence 
spectra  (indicated  by  arrows),  we  obtain  a  Franck- 
Condon  shift  of  4050  cm”'.  Indeed,  such  a  large 
Franck-Condon  shift  in  PTCDA  is  observed  in  MQ\N 
samples  as  well  as  in  the  bulk  PTCDA  crystal  To  our 
knowledge,  this  is  the  largest  value  reported  for  a 
Franck-Condon  shift  in  organic  crystals,  and  is  indicative 
of  an  extremely  strong  exciion-phonon  coupling  along 
the  tightly  packed  PTCDA  molecular  stacks. 

Time-resolved  PL  measurements  were  made  by  moni¬ 
toring  the  time  decay  of  the  luminescence  signal  at  720 
nm  (the  1-0  singlet  transition).  The  PL-signal  decay 
transients  for  two  MQW  samples  obtained  at  20  K  are 
shown  in  the  inset  of  Fig.  3.  The  inset  shows  the  re¬ 
sponse  for  the  two  samples  under  equal  excitation  inten¬ 
sity.  Here,  the  exponential  decay  over  several  decades  in 
intensity  of  the  transient  is  apparent.  In  addition,  the  PL 
decay  time  is  found  to  decrease  as  the  layer  thickness  de¬ 
creases.  The  dependence  of  the  exciton  lifetime  on  well 
width  is  shown  in  Fig.  3.  The  exciton  lifetime  ij,  is 
found  to  be  I0.8±0.5  ns  f^r  Z-”20  nm,  decreasing  to 
5.7  ±0.5  ns  for  Z..- “  1  nm.  The  change  in  exciton  life¬ 
time  in  the  MQW  samples  could  be  due  to  traps  at  the 
interface.  The  overall  recombination  rate  can  be  written 
as  r  r,” '  +  Tn7 where  rr  is  the  exciton  radiative  life¬ 
time  and  T„,  is  the  trapping  time.  If  the  nonradiative 
recombination  rate  is  such  that  Tn,  ^  r„  the  observed 
lifetime  should  decrease  due  to  traps.  The  results  of  our 
PL  measurements  show  that  the  time-integrated  lumi¬ 
nescence  intensity  was  found  to  be  constant  with  sample 
layer  thickness,  indicating  that  nonradiative  recombina¬ 
tion  is  not  important.  Further,  if  the  exciton  lifetime  in 
molecular  crystals  is  dominated  by  traps,  it  should  de¬ 
crease  with  increasing  temperature.  '*  However,  our 
time-resolved  measurements  show  that  the  radiative  de¬ 
cay  time  is  independent  of  temperature  from  20  to  295 
K,  indicating  again  that  the  PL  signal  here  is  not  due 
traps. 

The  decrease  of  t„  with  decreasing  Z..-,  therefore,  can 
result  from  an  increase  in  the  spontaneous  recombination 
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FIG.  3.  Exciton  liretime  r,,  as  a  function  of  well  width  The 
solid  line  is  the  fit  to  the  data  based  on  the  theory  described  in 
the  text.  Inset:  PL-signal  time-decay  transients  for  two  MQW 
samples.  Note  that  the  transient  for  the  L:  “200  A  sample  has 
been  shifted  upward  by  10%  of  a  decade  in  this  plot. 

rate  due  to  an  increase  in  the  overlap  of  the  electron  and 
hole  wave  functions  in  the  PTCDA  potential  wells.  A 
similar  explanation  of  the  observed  dependence  of  the 
exciton  lifetime  on  well  width  has  been  observed  in 
GaAs/GaAIAs  quantum  wells. For  the  hydrogenic  |j 
state,  /?,p  can  be  shown  to  be  proportional  to  the  exciton 
volume ’’  (given  by  which  differs  signifi¬ 

cantly  from  adL:  used  in  previous  work’^).  Using  the 
results  of  the  calculation  of  a  discussed  above  (Fig.  2), 
the  exciton  volume  in  a  MQW  is  observed  to  decrease 
with  decreasing  well  width.  The  shrinkage  of  the  exciton 
due  to  quantum  confinement  can  thus  account  for  the  de¬ 
crease  in  exciton  lifetime.  Based  on  the  above  analysis 
and  using  the  parameters  given  in  Table  I,  the  depen¬ 
dence  of  exciton  volume,  and  hence  r,,,  on  well  width  is 
shown  by  the  solid  line  in  Fig.  3.  The  results  of  the  cal¬ 
culation  are  in  good  agreement  with  the  exciton  lifetime 
data,  indicating  the  exciton  quantum-confinement  model 
can  quantitatively  account  for  the  change  in  both  exciton 
binding  energy  and  exciton  lifetime  using  the  model 
given  above  and  the  parameters  in  Table  I. 

In  summary,  we  have  characterized  PTCDA/NTCDA 
organic  quantum-well  structures  by  both  optical-ab¬ 
sorption  and  time-resolved  photoluminescence  measure¬ 
ments.  The  exciton  binding  energy  was  found  to  in¬ 
crease  with  decreasing  well  width,  whereas  its  lifetime 
decreased  with  decreasing  well  width.  The  dependence 


of  exciton  binding  energy  and  lifetime  on  can  both  be 
understood  ii  terms  of  confinement  of  Wannier-like 
molecular  ex'itons  in  quantum  wells.  To  our  knowledge, 
these  results  provide  the  first  evidence  for  exciton  quan¬ 
tum  confinement  in  organic  MQW  structures.  From  the 
results  of  our  calculation,  the  effective  exciton  Bohr  ra¬ 
dius  is  12  A,  indicating  that  the  exciton  wave  function 
extends  to  three  or  four  molecular  layers  in  PTCDA. 
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We  demonstrate  all-optical  modulation  in  crystalline  organic  semiconductor  waveguides 
grown  by  the  ultrahigh  vacuum  process  of  organic  molecular-beam  deposition.  Two  light 
beams  with  wavelengths  of  1.06  and  0.514  /am  were  used  as  the  guided  and  the  pump 
light  sources,  respectively.  A  refractive  inde.x  change  of  5.4  X  10“^  at  1.06  /am  was  observed 
at  a  pump  intensity  of  1.0  W/cm‘.  This  large  nonlinear  effect  is  attributed  to  free 
electron-hole  pairs  produced  by  the  dissociation  of  excitons  generated  by  the  short  wavelength 
beam.  A  carrier  lifetime  of  (17±  1)  /is  which  determines  the  modulator  switching  time 
is  in  good  agreement  with  theoretical  predictions.  To  our  knowledge,  this  is  the  first 
observation  of  free-carrier-induced  inde.x  modulation  in  crystalline  organic  waveguides. 


For  more  than  two  decades,  organic  materials  have 
been  investigated  for  their  potential  usefulness  to  nonlinear 
optics  (NLO).’  Most  research  has  been  focused  on  non¬ 
resonant  nonlinear  effects  since  the  losses  are  smaller  and 
the  response  times  are  shorter  than  at  resonance.  Resonant 
NLO  effects,  however,  are  interesting  because  they  can  be 
induced  at  very  small  optical  powers.^  Resonant  NLO  ef¬ 
fects  also  clarify  the  nature  of  optical  transitions  and  dy¬ 
namic  processes  which  are  not  yet  completely  understood 
for  organic  materials.  In  previous  work,^  the  linear  dielec¬ 
tric  properties  of  the  crystalline  organic  compound,  3,  4,  9, 
10  perylenetetracarboxylic  dianhydride  (PTCDA)  were 
reported.  In  this  letter,  we  report  on  the  demonstration  of 
resonant  NLO  effects  associated  with  free  carriers 
(plasma)  generated  by  dissociated  excitons  in  PTCDA  rib 
waveguides.  To  our  knowledge,  this  is  the  first  observation 
of  light-generated  plasma  effects  in  organic  crystals. 

Charge  transport  in  molecular  semiconductors  is 
thought  to  be  initiated  via  excitonic  transitions.^  Here,  free 
electron-hole  pairs  are  generated  in  a  second-order  process 
via  dissociation  of  photogenerated  e.xcitons  at  surfaces  or 
bulk  material  impurities.  That  is,  given  a  molecule  in 
ground  state  Sq,  a  dominant  channel  for  free-carrier  pro¬ 
duction  due  to  absorption  of  light  of  energy  ficj,  is 
5o-f  5*,  followed  by  S*-hM^e~  +  h'^ .  Here  5*  is 
the  excited  singlet  exciton  state  of  the  molecule,  M  repre¬ 
sents  an  impurity  or  interface,  and  e~  and  h  are  the  free 
electron  and  hole,  respectively.  The  absorbed  light  must 
have  a  photon  energy  greater  than  the  exciton  transition 
energy  (typically  1.5-3  eV).  Since  the  process  is  second 
order,  the  quantum  efficiency  (t;)  for  electron-hole  gener¬ 
ation  is  usually  <1%.  Assuming  that  there  is  a  refractive 
index  change  (An)  associated  with  the  optical  generation 
of  free  carriers  (via  excitons),  the  Drude  model  gives 


^rt  = 


-e^S 


(I) 


where  e  is  the  electronic  charge.  .V  is  the  free-carrier  den¬ 
sity  induced  by  the  optical  beam,  n  is  the  refractive  index  in 
the  dark,  m*  is  the  carrier  effective  mass,  fg  is  the  permit¬ 
tivity  in  vacuum,  and  to  is  the  light  frequency.  A  steady- 
state  solution  of  the  continuity  equation  for  A’  gives 


^n  =  n-,I  = 


—  erarT]! 
2nm*eofuo-  ' 


(2) 


Here,  a  is  the  absorption  coefficient,  I  is  the  optical  inten¬ 
sity,  T  is  the  free-carrier  lifetime,  and  ft  is  Planck's  constant 
divided  by  2n-.  Note  that  the  NLO  coefficient,  n->.  given  by 
Eq.  (2),  is  only  due  to  the  free-carrier  concentration  (not 
the  exciton  population),  and  is  intensity  dependent. 

To  fabricate  the  waveguide  devices  in  which  rii  was 
measured,  a  l-/xm-thick.  AZ  1400  photoresist  layer  (with 
/ip  =  1.61 )  was  first  spun  onto  the  surface  of  a  cleaned  and 
polished  ( 100)  InP  substrate.  Next,  a  series  of  2-/im-wide 
strips  were  patterned  onto  the  photoresist  along  the  (110) 
direction  by  standard  photolithographic  techniques.  A  1- 
/rm-thick  crystalline  PTCDA  layer  was  then  deposited 
onto  both  the  InP  substrate  and  the  photoresist  ridges  us¬ 
ing  the  ultrahigh  vacuum  ('lO"’  Torr)  process  of  or¬ 
ganic  molecular-beam  deposition.^  Deposition  proceeds  by- 
heating  a  prepurified  source  of  PTCDA  to  450  °C  to 
achieve  a  deposition  rate  of  approximately  3  A/s.  while  the 
substrate  temperature  is  maintained  at  90  K.  It  has  previ¬ 
ously  been  found  that  thin-film  deposition  under  these  con¬ 
ditions  results  in  single-crystalline  growth  of  PTCDA  on 
the  photoresist  strips  into  low-loss  (  <2.5  dB),  TE-mode 
waveguides.^  To  form  waveguide  facets,  the  InP  substrate 
was  cleaved  along  the  (110)  direction.  The  facet  surfaces 
are  sufficiently  smooth  to  allow  for  observation  of  Fabry- 
Perot  resonances  as  the  index  of  the  guide  is  thermally  or 
optically  varied. 

In  the  experiments,  1.06-/um  wavelength  light  from  a 
YAG  laser  was  coupled  into  and  out  of  the  waveguide 
using  two  microscope  objective  lenses  (see  inset.  Fig.  1 ).  A 
high-speed  Ino  53Gag47As  avalanche  photodetector  (APD) 
was  used  to  detect  the  light  signal.  To  ensure  a  constant 
coupling  efficiency,  a  T'V  camera  monitored  the  beam  spot 
position  using  a  beam  splitter  placed  in  front  of  the  APD 
The  waveguide  was  placed  on  a  copper  block  mounted  on 
a  thermoelectric  cooler.  Using  a  thermocouple  and  a  tem¬ 
perature  controller,  the  temperature  of  the  waveguide  de¬ 
vice  was  controlled  to  within  ±0.05  "C. 

The  2-/im-wide  PTCDA  nb  waveguide  allowed  for 
propagation  of  a  single  TEgg  mode.  The  waveguide,  with  its 
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FIG.  1.  Temperature-dependent  intensity  measured  in  a  PTCDA  Fabry- 
Perot  waveguide  resonator  at  A  =:  1.06  nm  with  a  length  of  1.62  mm  and 
a  width  of  2  /im.  Inset:  The  experimental  setup.  1;  argon  laser  (A  =  0.114 
/xm);  2:  AO  modulator:  3:  ring  YAG  laser  (A  =  1.06  /am);  4.  Glan- 
Thompson  polarizer,  5:  thermoelectric  cooler  and  thermocouple;  6.  TV 
camera:  7;  monitor,  8;  iris;  9;  APD;  10:  oscilloscope. 

cleaved  facets,  acts  as  a  lossy  Fabry-Perot  resonator  with  a 
transmission  function  of 

^^^^“(l-/?e-‘^)^  +  4i?e-“^sin-(5/2)  ’ 

where  a  is  the  absorption  coefficient,  L  is  the  device  length 
(1.62  mm),  R  is  the  reflection  coefficient  at  the  waveguide 
facets  which  is  calculated  to  equal  0.11  considering  only 
Fresnel  reflections,  and  5  is  the  phase  shift  defined  as 
5  =  5o  +  A(5  =  (4trZ./A)(n-fAn).  When  the  temperature 
T  is  varied,  the  transmitted  intensity  varied  periodically 
(see  Fig.  1 )  as  is  characteristic  of  a  Fabry-Perot  resonator. 
A  complete  fringe  (two  adjacent  maxima)  corresponding 
to  a  phase  difference  of  A5  =  2-  was  observed  for  a  tem¬ 
perature  change  of  2  'C,  from  which  we  obtain 
(1/«)A«/Ar=  (8.1  ±  0.5)  X  lO-'K"'.  Using  Eq.  (3). 
however,  a  maximum  transmission  modulation  of  ~509c 
was  expected,  compared  to  an  observed  modulation  of  only 
10*:^:.  This  smaller-than-expected  modulation  is  possibly 
due  to  imperfect  waveguide  facets  which  lower  the  finesse 
of  the  resonator. 

To  measure  the  all-optical  modulation  characteristics 
of  the  guide,  a  pump  light  beam  of  0.514-)um  wavelength 
with  a  diameter  L'  =  I  mm  directly  illuminated  the  surface 
of  the  PTCDA  waveguide,  while  the  1.06-/im  beam  was 
end-fire  coupled  into  the  waveguide,  as  in  the  case  of  the 
In/ST  measurements.  In  order  to  avoid  thermally  in¬ 
duced  intensity  changes  from  the  pump,  the  device  tem¬ 
perature  was  precisely  mainuined  at  22.1  ‘C  using  the 
thermoelectric  cooler.  The  pump-induced  intensity 
changes  are  shown  in  Fig.  2.  Here,  the  transmitted  inten¬ 
sity  is  decreased  nearly  linearly  proponional  to  the  pump 
intensity  at  a  temperature  of  22. 1  ’C.  In  contrast,  the  trans¬ 
mitted  intensity  due  to  a  temperature  nse  at  22.1  *C  is 
increased  ( Fig.  1 ) . 
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FIG.  2.  Intensity  of  guided  light  (A  =  1.06  nm)  vs  pump  lighi  iniensi!> 
(A  =  0.514  /im)  in  a  PTCDA  waveguide  with  the  same  parameters  in 
Fig.  1. 

While  we  suggest  that  the  pump  power  dependence  of 
P{8)  is  due  to  an  optically  induced  index  change,  it  is 
nevertheless  important  to  consider  other  competing  pro¬ 
cesses.  The  most  significant  contributions  to  An  are  ther¬ 
mally  induced  index  changes  and  optically  induced  absorp¬ 
tion  effects.  To  differentiate  between  these  various 
mechanisms,  we  have  investigated  the  modulation  behavior 
at  very  low  pump  light  intensities  ( <  1  W/cm^)  over  a 
wide  temperature  range  ( 16-37  *C).  Three  observations 
exclude  the  possibility  of  thermal  effects:  First,  the  changes 
in  transmitted  intensity  under  a  constant  pump  power 
could  either  be  increased  or  decreased,  depending  on  tem¬ 
perature.  However,  the  sign  of  the  optically  induced 
change  of  transmitted  intensity  was  always  opposite  to  the 
thermally  induced  changes  over  the  entire  temperature 
range  explored.  For  instance,  at  a  cenain  temperature,  if 
the  transmitted  intensity  was  decreased  due  to  pump  light 
illumination,  it  was  increased  due  to  a  temperature  rise. 
This  is  similar  to  observations  for  inorganic  semi¬ 
conductors*  in  which  nonlinear  effects  were  attributed  to 
index  changes  due  to  excitons. 

A  further  observation  ruling  out  thermal  effects  is  that 
the  temperature  increase  is  calculated  to  be  —0.03  ’C  un¬ 
der  a  maximum  pump  light  intensity  of  1  W/cm^  This 
calculation  assumes  a  steady-state  thermal  conductivity  of 
0.1,  0.007,  and  0.8  W/cm  K  for  PTCD.A..  AZ1400,  and 
InP,  respectively.  Such  a  small  temperature  rise  does  not 
significantly  contribute  to  the  intensity  modulation,  as  in¬ 
ferred  from  Fig.  1.  Finally,  the  transient  thermal  rela.\ation 
time  was  calculated  to  be  >13  ms.  which  is  three  orders  of 
magnitude  larger  than  the  measured  modulation  time  of  17 
ps.  as  shown  below.  A  very  low  amplitude  "tail"  in  the 
modulation  response  with  a  time  constant  of  —6  ms  was 
observed  at  the  highest  pump  intensities.  We  attribute  this 
latter  time  constant  to  be  characteristic  of  very  small  ther¬ 
mal  effects,  whereas  the  shoner  time  constant  results  from 
significantly  larger  carrier  lifetime  effects. 

The  second  "parasitic"  effect  to  be  considered  is  pump 
induced  absorption  changes.  Aa.  Using  Eq.  (3;.  it  can  be 
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FIG.  3.  Rise  and  fall  limes  of  the  1.06-/im  wavelength  beam  due  to 
modulation  of  the  A  =  O.J14-jim  beam. 

shorn,  that  the  change  of  transmitted  Jight  intensity,  AP,  in 
the  limit  of  aZ.<l,  is  given  by  AP  =  -  A Aa,  where  is  a 
constant.  Since  the  waveguide  used  in  our  experiment  had 
aL  =  0.092,  this  expression  is  valid.  Thus,  if  the  optically 
induced  intensity  change  is  due  to  Aa,  we  should  see  a 
constant  decrease  of  the  transmitted  intensity  under  a  con¬ 
stant  pump  light  illumination,  independent  of  temperature. 
However,  as  noted  above,  the  sign  of  the  observed  change 
varied,  dependent  on  temperature.  In  addition,  the  absorp¬ 
tion  change  due  to  the  opticaUy  induced  plasma  is  given  by 
Aa/An  =  2/ct,  where  r  is  the  carrier  lifetime  and  c  is  the 
velocity  of  light.  For  PTCDA,  Aa/An  =  4.2x  10“*  cm~‘, 
and  absorption  changes  are  insignificant.  We  conclude, 
therefore,  that  the  modulation  is  due  to  an  intensity  depen¬ 
dent  refractive  index,  rather  than  a  change  in  the  absorp¬ 
tion  coeflBcient. 

Since  both  the  optically  and  the  thermally  induced 
modulation  is  caused  by  the  changes  in  refractive  index, 
the  thermal  measurements  can  be  used  Co  determine 
An  =  nil.^  That  is,  since  the  thermal  measurements  indi¬ 
cated  that  an  output  intensity  change  of  10%  corresponds 
to  a  ir-phase  change  (Fig.  1),  then  the  2%  decrease  in 
transmitted  intensity  induced  by  a  pump  power  of  /  =  0.98 
W/cm^  corresponds  to  A5  =  tt/5.  Using  the  relationship 
A<5  =  4iTin2l)L'/A,  we  get  n,  =  5.4 X  10~^  cm^/W,  where 
L'  =  I  mm  is  the  interaction  length.  Here,  22. 1  *C  was 
chosen  as  the  measurement  temperature  since  the  optically 
induced  intensity  change  is  nearly  linear  over  a  small  range 
of  pump  intensity  (  <  I  W/cm^).  Of  course,  the  measure¬ 
ment  can  be  made  equally  well  at  other  temperatures. 

By  comparing  the  measured  carrier  lifetime  to  theoret¬ 
ical  predictions,  we  verified  that  R2  indeed  due  to  an 
optically  generated  plasma.  Here,  t  was  directly  measured 
from  the  transient  response  time  of  the  l.Ob-^um  wave¬ 
length  beam  due  to  a  pulsed  0.5l4-/xm  pump  beam.  Using 
a  measurement  system  with  a  response  time  of  <  I  fis,  the 
rise  and  fall  times  of  the  l.06-/im  beam  were  found  to  be  18 
and  16  /xs,  respectively  (Fig.  3).  For  comparison,  from  the 
steady-state  measurements  of  R;,  using  Eq.  (2),  we  obtain 
T=  (14±  1)  /IS,  which  is  very  close  to  the  measurement. 
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The  constants  used  in  the  calculation  for  PTCDA  are 
a  =  4.6  X  10*  cm“'  at  A  =  0.514  pm  as  obtained  from  ab¬ 
sorbance  measurements,  a  hole  effective  mass  for  PTCD.A 
of  m*  =  0.1 8mo  (where  mg  is  the  electron  rest  mass),*  and 
a  quantum  efficiency  of  tj  =0.01.*  The  presence  of  elec¬ 
trons  was  ignored  since  their  effective  mass  is  at  least  an 
order  of  magnitude  larger  than  that  of  holes.*  In  this  case, 
Eq.  (2)  implies  that  their  effect  on  An  is  negligible. 

In  Table  I  we  list  the  resonant  nonlinear  propenies  of 
some  typical  inorganic  semiconductors,  along  with  one  or¬ 
ganic  polymer,  and  PTCDA.  From  this  Table,  r,  of 
PTCDA  is  one  of  the  highest  observed  to  date.  However, 
the  response  time  (r)  is  comparatively  long.  To  reduce  the 
carrier  lifetime  for  high-speed  applications,  organic  multi¬ 
ple-quantum-well  structures*  are  potentially  useful.  How¬ 
ever,  decreasing  the  carrier  lifetime  might  also  reduce  tj, 
thereby  ultimately  decreasing  r^.  Alternatively,  one  can 
reduce  r  by  electrically  injecting  holes  from  conucts,  and 
then  sweep  them  out  using  an  applied  electric  field. 

In  the  conclusion,  we  report  the  apparently  first  obser¬ 
vation  of  free-carricr-induced  nonlinear  effects  in  crystal¬ 
line  organic  semiconductors.  An  intensity-dependent  re¬ 
fractive  index  of  R2  =  5.4x10"*  cm~‘/W  has  been 
measured,  and  is  among  the  highest  values  reported  for 
practical  optical  device  materials. 
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We  directly  observe,  using  the  scanning  tunneling  microscope,  a  two-dimensional  crystal  of  the 
organic  compound  3,  4,  9,  10-peiylenetetracarboxylic  dianhydride  (PTCDA).  The  surface  unit 
cell  dimension  is  found  to  be  21.6^2.2  A  by  1S.2±  1.6  A,  or  approximately  20%  larger  than  the 
bulk  unit  cell.  Furthermore,  the  organic  lattice  is  oriented  with  respect  to  the  graphite  substrate 
even  though  the  two  lattices  are  incommensurate.  These  obse^ations  are  consistent  with 
reflection  high  energy  electron  diffraction  measurements,  and  energy  minimization  calculations, 
assuming  that  the  van  der  Waals  bond  is  the  predominant  intermolecular  force  which 
determines  the  equihbrium  crystal  structure.  The  combination  of  measurement  and  theory 
provides  the  first  step  in  developing  tools  for  predicting  the  conditions  which  lead  to 
quasiepitaxial  growth  of  these  technolo^cally  important  van  der  Waals  solids. 

The  archetype  crystalline  organic  seniiconductor  3,  4, 

9,  lO-perylenetetracarboxylic  dianhydride  (PTCDA  or 
C2406Hg),  has  many  practical  device  applications  such  as 
waveguides,’  couplers,^  and  organic/inorganic  semicon¬ 
ductor  heterojunction  devices.^  Since  the  organic  crystal  is 
bonded  by  weak  van  der  Waals  (vdW)  forces,  single¬ 
crystal  thin  films  can  be  deposited  on  inorganic  substrates 
without  regard  for  lattice  matching;  a  process  known  var¬ 
iously  as  “van  der  Waals  epitaxy”  or  “quasiepitaxy”.  Un¬ 
derstanding  this  new  growth  process  is  crucial  if  we  are  to 
fully  exploit  its  potential  for  use  in  a  wide  range  of  device 
structures. 

In  this  work  we  use  the  scanning  tunneling  microscope 
(STM)  to  study  the  crystal  structure  at  the  onset  of  film 
growth  by  examining  the  first  monolayer  of  PTCDA  de¬ 
posited  by  the  ultrahigh  vacuum  process  of  organic  molec¬ 
ular  beam  deposition  (OMBD)  on  highly  oriented  pyro¬ 
lytic  graphite  (HOPG)  substrates  held  at  low  temperature 
(90  K).  We  show  that  the  surface  unit  cell  is  expanded 
from  that  found  in  the  bulk,  and  is  also  distorted  by  rota¬ 
tions  of  individual  molecules  with  respect  to  the  unit  cell. 

Furthermore,  we  determine  an  orientational  relationship 
between  the  PTCDA  layer  and  the  graphite  substrate.  This 
is  supported  by  calculations  which  minimize  total  crystal 
energy  to  determine  the  equilibrium  positions  of  molecules 
within  the  two-dimensional  unit  cell  of  PTCDA,  and  is 
consistent  with  recent  observations  obtained  using  reflec¬ 
tion  high  energy  electron  diffraction  (RHEED).* 

A  perspective  view  along  the  c-axis  of  the  monoclinic 
bulk  structure  of  PTCDA  is  shown  in  Fig.  1.  The  material 
forms  planar  stacks  with  molecular  spacing  of  3.21  A,  and 
has  lattice  constants  of  a =3.72  A,  6=11.96  A,  and  c 
=  17.34  A,  and  a  lattice  angle  ^=98.8*.  Thick  films  (ca. 

1000  A)  of  PTCDA  have  long  range  crystalline  order 
when  deposited  at  high  rates,’  or  on  substrates  held  at  low 
temperature.'  Scanning  tunneling  microscopy  has  been 
shown  to  be  capable  of  molecular  resolution  of  organic 

“Current  address:  DeparUneni  of  Electrical  Engineering,  Princeton 
University.  Princeton.  NJ  08J44.  FIG.  1.  The  PTCDA  bulk  unit  cell  viewed  along  the  c-axis. 

2417  Appl.  Phys.  Lett  61  (20).  16  November  1992  0003-695 1/92/452417-03$03.00  ©  1992  American  Institute  of  Physics  2417 


molecules.*’’  Ludwig  et aL*  recently  presented  molecular 
resolution  images  of  PTCDA  deposited  on  HOPG  held  at 
340-400  K.  They  explained  their  results  in  terms  of  a  sur¬ 
face  unit  cell  similar  to  the  bulk  unit  cell,  but  S%-10% 
larger.  However,  their  films  were  deposited  under  condi¬ 
tions  far  from  those  shown  to  produce  single  crystalline 
layers  of  PTCDA.'-’  ’ 

The  OMBD  vacuum  system  used  to  prepare  our  sam¬ 
ples  is  described  elsewhere.’  A  1.4  cm’  HOPG  substrate 
was  cleaved  in  air  and  immediately  loaded  into  the  depo¬ 
sition  chamber  via  a  vacuum  load  lock.  The  PTCDA,  pre¬ 
purified  by  rqieated  temperature  gradient  sublimation,  was 
loaded  into  a  Knudsen  cell  and  thoroughly  degassed  at 
200 ’C  prior  to  film  deposition.  Approximately  4  A  (1 
monolayer)  of  PTCDA  was  deposited  at  a  rate  of  0. 1-0.2 
A/s,  measured  by  a  quartr  'nystal  thickness  monitor.  Such 
a  measured  thickness  is  consistent  with  the  disappearance 
of  the  substrate  streaks  in  the  RHEED  experiments.* 

The  deposited  film  was  imaged  in  air  using  a  commer¬ 
cially  available  STM'°  with  controlled  geometry,  Pt/Ir 
tips"  within  about  24  h  of  growth.  Current  images  (con¬ 
stant  height  mode)  apparently  showing  single  molecules  of 
PTCDA  were  readily  obtained  on  several  samples  at  tun- 
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FIG.  2.  (a)  False  color  filtered  STM  image  of  crystalline  PTCDA  on 
HOPG.  Tip  operating  parameters:  402.1  mV,  O.IO  nA  tip  positive,  (b) 
Unfiltered  STM  image  of  the  edge  of  an  island  of  PTCDA  on  HOPG.  Tip 
operating  parameters:  403.7  mV,  0.10  nA  tip  positive,  (c)  Detail  of  (a), 
showing  the  unit  cell  of  PTCDA  on  HOPG.  Superimposed  lines  are  a 
guide  to  the  eye. 

neling  currents  of  ~0.1  nA  and  a  tip  voltage  of  — 1-400 
mV  with  respect  to  the  substrate.  No  images  were  obtained 
with  the  tip  biased  negative  with  respect  to  the  substrate.  A 
typical  image,  filtered  to  reduce  noise,  is  shown  in  Fig. 
2(a).  Such  images  were  stable  for  —10-20  scans  before 


FIG.  3.  Proposed  structure  for  the  two-dimensional  unit  cell  of  PTCDA 
on  HOPG,  showing  the  parameters  measured  and  displayed  in  Table  1 

being  replaced  by  an  image  of  the  underlying  graphite  lat¬ 
tice.  Since  higher  tunneling  currents  result  in  a  higher  rate 
of  image  degradation,  this  is  attributed  to  the  weakly 
bound  organic  molecules  being  swept  away  by  the  tip. 

Figure  2(b)  shows  an  example  unfiltered  image  of  the 
edge  of  an  island  of  PTCDA;  the  individual  molecules  ap¬ 
pearing  as  broken  double  circles.  Coupled  with  the  high 
reproducibility  of  the  images  between  numerous  different 
samples  and  using  different  tips,  the  observation  of  island 
edges  provides  evidence  that  the  images  presented  are,  in 
fact,  of  a  monolayer  of  the  organic  thin  film  on  graphite. 

Based  on  the  STM  images,  we  propose  a  model  for  the 
two-dimensional  unit  cell  of  PTCDA  on  HOPG  shown  in 
Fig.  3,  with  the  measured  dimensions  and  angles  between 
the  molecules  listed  in  Table  I.  The  linear  dimensions,  a 
and  b,  are  measured  and  calibrated  relative  to  an  image  of 
the  HOPG  lattice.  The  angle  B  is  measured  directly  from 
the  images  taken  from  several  samples  studied,  whereas  the 
angle  a  is  measured  using  Fig.  2(c).  This  is  a  magnified 
detail  of  the  unit  cells  in  Fig.  2(a)  showing  the  fifth,  cen¬ 
tral,  molecule  of  the  unit  cell,  which  is  faintly  apparent  in 


TABLE  I.  PTCDA  crystal  dimensions  and  angles. 
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11  96 
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the  low  magnification  image  as  well.  The  positions  of  the 
long  molecular  axes  are  indicated  by  the  line  segments 
shown  in  the  figure. 

Lattice  dimensions  and  angles  measured  from  images 
of  different  samples  with  different  tips  are  reproducible 
within  the  limits  quoted  in  Table  1.  Two  PTCDA  mole¬ 
cules  define  a  rectangular  unit  cell  with  dimensions  a 
=  21.6±2.2  A  and  6=  !5.2±  1.6  A,  in  contrast  to  the  bulk 
unit  cell  with  dimensions  a  =  17.34  A  and  6  =  1 1.96  A.  The 
angles  6  and  a  are  55  ±2*  and  35  ±8°,  respectively.  These 
parameters  agree  with  those  recently  reported  using 
RHEED  (Table  I). 

The  apparent  removal  of  the  organic  molecules  by  the 
tip  allows  us  to  examine  the  orienutional  relationship  be¬ 
tween  the  organic  film  and  the  graphite  substrate.  From 
these  images  we  find  the  angle  between  the  graphite  and 
PTCDA  lattices,  shown  as  in  Fig.  3,  to  be  17 ±4°. 

Some  larger  area  scans  show  adjacent  islands  of 
PTCDA.  The  angle  between  equivalent  PTCDA  molecules 
in  different  islands  is  always  an  integer  multiple  of  60“. 
Since  the  graphite  substrate  has  sixfold  rotational  symme¬ 
try,  nucleating  islands  of  PTCDA  are  energetically  identi¬ 
cal  at  rotations  of  60“  with  respect  to  the  substrate.  The 
absence  of  any  other  angles  in  the  observed  STM  images 
strongly  suggests  that  even  though  the  two  lattices  are 
highly  incommensurate,  they  align  in  a  consistent  fashion 
with  respect  to  each  other;  striking  evidence  for  quasiepi- 
taxial  growth.  This  represents  a  fundamental  limitation  of 
quasiepitaxy:  if  the  symmetries  of  the  two  lattices  are  not 
identical,  the  deposited  layer  must  contain  high  angle  grain 
boundaries  between  islands  which  have  separately  nucle¬ 
ated  at  different  lattice  sites. 

To  interpret  the  STM  images,  we  have  calculated  the 
crystal  energy  of  a  monolayer  of  PTCDA  by  assuming  that 
the  dominant  intermolecular  and  substrate-molecule  bond 
is  due  to  the  vdW  interaction  ( true  for  nonpolar  molecules 
such  as  PTCDA  which  lie  flat  on  the  substrate.  Thus,  the 
interaction  potential,  u^,  between  the  ith  atom  in  one  mol¬ 
ecule  and  the  yth  atom  in  its  neighbor  is  derived  from  the 
Buckingham  potential,  i.e.; 

“i/ 

u,j=  —r+l3,je\p( 

where  is  the  distance  between  the  ith  and  yth  atoms,  and 
a,j,  Pij,  y,j  are  the  vdW  coefficients  between  atom  pairs 
involved  in  bonding  PTCDA  to  its  neighbors  and  to  the 
substrate.'^  ’"’  Periodic  boundary  conditions  are  assumed  so 
that  only  potentials  within  the  unit  cell  are  considered.  The 
total  interaction  potential  is  Uj-,  obtained  by  summing  u,j 
over  all  atoms  in  the  unit  cell  and  substrate.  The  favored 
crystal  configuration  is  then  given  by 

^[C/r(x,)]=0, 

where  ranges  over  the  independent  parameters,  a,  b,  6, 
and  d>  (Fig.  3).  Because  of  the  short  range  of  the  van  der 
Waals  potential  (  —  l/r**),  distant  molecular  interactions 
are  negligible  and  the  model  is  rapidly  tractable. 


Unit  cell  parameters  of  the  bulk  crystal  calculated  us¬ 
ing  this  model  agree  with  those  measured  by  X-ray  diffrac¬ 
tion  to  within  2%.  In  addition,  the  calculated  crystal  sub¬ 
limation  energy  of  ~30  kcal/mol  is  consistent  with 
observations  for  PTCDA  and  other  anthracene-like  crys¬ 
tals.  Therefore,  application  of  the  model  to  two- 
dimensional  surface  cells  appears  to  be  justified. 

Results  of  the  calculations  applied  to  the  PTCDA  sur¬ 
face  unit  ceil  are  also  given  in  Table  I.  The  calculated 
equilibrium  distance  of  the  PTCDA  layer  to  the  graphite 
substrate  is  3.21  A;  a  value  equal  to  the  intermolecular 
spacing  in  bulk  PTCDA.  The  close  agreement  between  a, 
b,  9,  and  a  predicted  by  the  model  and  those  measured 
using  STM  and  RHEED  clearly  supports  the  structure 
shown  in  Fig.  3.  For  example,  the  calculated  rotational 
orientation  of  <^=11“  is  energetically  favored  by  ~100 
meV  over  other  angles. 

From  these  model  calculations,  we  infer  that  the  en¬ 
larged  reconstructed  surface  unit  cell  is  a  direct  result  of 
the  lack  of  vdW  bonds  extending  above  the  top  layer,  al¬ 
lowing  the  cell  to  expand  slightly  to  minimize  the  repulsive 
core  energy  between  adjacent  atoms.  We  note,  however 
that  this  surface  reconstruction  does  not  depend  strongly 
on  the  substrate,  but  rather  is  a  result  of  the  comparatively 
strong  intermolecular  forces  within  the  cell  itself.  A  more 
detailed  discussion  of  these  results  will  be  presented  else¬ 
where. 

In  summary,  the  STM  images  of  the  surface  unit  cell 
are  in  close  agreement  both  with  in  situ  RHEED  data 
taken  for  films  grown  under  simttar  conditions  and  with 
crystal  energy  minimization  calculations.  Such  close  agree¬ 
ment  implies  that  the  calculations  should  be  predictive  of 
the  optimal  conditions  under  which  quasiepitaxial  growth 
can  be  obtained  for  a  wide  range  of  van  der  Waals-bonded 
materials  with  a  broad  range  of  applications  to  optoelec¬ 
tronics. 
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A  new  class  of  engineered  materials  based  on  crystalline  organic  semiconductors  grown  by  the  ultra-high 
vacuum  process  of  organic  molecular  beam  deposition  is  described.  We  discuss  the  conditions  which  lead 
to  ordered  growth  between  highly  lattice-mismatched  van  der  Waals  crystals  into  novel  ‘quasi-epitaxial’ 
structures.  Furthermore,  a  model  is  presented  which  describes  the  physical  mechanisms  underlying  quasi¬ 
epitaxy.  The  elearical  and  optical  properties  of  both  single  and  multiple  heterojunction  structures  grown 
using  alternating  layers  of  two  different,  lattice-mismatched  organic  molecules  are  considered  in  detail. 
It  is  found  that  the  layers  can  be  grown  sufficiently  thin  (10  A)  to  result  in  exciton  confinement  in  organic 
multiple  quantum  wells.  These  and  other  quasi-epitaxiaJly  grown  structures  can  lead  to  an  entirely  new 
family  of  optoelectronic  devices. 


Keywords:  molecular  beam  epitaxy,  nonlinear  optics,  quantum  well 


INTRODUCTION 

Recent  work  in  many  laboratories  worldwide''^  has  demonstrated  a  new  class  of 
engineered  materials  with  application  to  photonic  devices.  These  materials  are 
variously  known  as  Van  der  Waals  solids’,  'layered  materials’  or  ‘quasi-epitaxial 
materials’.  The  property  these  materials  have  in  common  is  that  the  cohesive  force 
which  bonds  the  various  material  layers  together  is  the  relatively  weak  van  der  Waals 
force.  This  results  in  an  ability  to  layer  materials  which  are  highly  lattice-mismatched 
without  inducing  defects  in  the  crystal  structure.  One  particularly  interesting  class 
of  van  der  Waals  solids  with  applications  to  a  very  broad  range  of  optical  and  elec¬ 
tronic  devices  is  based  on  multilayered  structures  consisting  of  crystalline  organic 
semiconductors.'  These  materials  have  been  found  to  have  excellent  electrical  and 
optical  characteristics  when  grown  by  the  ultra-high  vacuum  process  of  organic 
molecular  beam  deposition  (OMBD)  onto  other  organic  semiconductors,  or  even 
onto  the  surfaces  of  inorganic  semiconductors  such  as  Si,  GaAs  or  InP.'-’  For 
example,  high  bandwidth  photodetectors  consisting  of  layers  of  perylene-based 
compounds  deposited  on  Si  substrates  have  been  demonstrated,'  as  have  Field-effect 
transistors  based  on  organic/lnP  heterojunctions. Furthermore,  waveguides  and 
other  optical  devices  have  been  demonstrated  using  these  quasi-epitaxially  grown 
layers." 

In  this  work,  we  discuss  the  structural,  dielectric,  conductive,  and  optical  proper¬ 
ties  of  this  new  and  exciting  class  of  engineered  materials.  In  particular,  we  have 
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FIGURE  1  Two  views  of  a  unit  cell  of  PTCDA. 


grown  which  extended  several  centimetres  across  the  substrate  surface.  Follow¬ 
ing  that  initial  work,  Debe  et  al.'^  demonstrated  that  other  perylene  derivatives 
similar  to  PTCDA  could  attain  single  crystallinity  when  deposited  more  slowly 
onto  substrates  held  well  below  room  temperature.  More  recently,  we  have  found 
that  PTCDA,  and  an  analogous  compound,  3, 4, 7,8  naphthalenetetracarboxylic 
dianhydride  (NTCDA),  or  C,40tH,,  can  both  be  grown  quasi-epitaxially  at  a  rate 
of  1-3  A /s  when  deposited  onto  glass  substrates  held  at  approximately  SK)  K.'  Low 
substrate  temperature  during  growth  is  required  since  the  thermal  energy  of  the 
molecules  on  the  surface  needs  to  be  brought  well  below  the  intermolecular  bond 
energy  for  ordered  growth  to  be  established. 

Multilayer  stacks  of  alternating  layers  of  PTCDA  and  NTCDA  also  have  recently 
been  found  to  form  ordered  crystalline  structures.  This  order  is  achieved  even 
though  the  crystal  structures  of  PTCDA  and  NTCDA  are  highly  mismatched. 
In  Figure  2  we  show  a  unit  cell  of  NTCDA  (cf.  Figure  1).  Here,  it  is  seen  that  rather 
than  forming  planar  stacks  as  in  the  case  of  PTCDA,  NTCDA  molecules  stack  in  a 
herringbone  manner,  with  an  intermolecular  spacing  of  3.S  A .  In  the  herringbone 
structure,  adjacent  stacking  axes  are  positioned  at  90°  with  respect  to  eacL  other, 
once  again  leading  to  anisotropies  in  the  conductive  and  dielectric  properties  of  these 
crystals.  However,  the  asymmetries  are  much  less  pronounced  than  those  found  in 
PTCDA. 

In  growing  the  multilayer  stacks  of  PTCDA  and  NTCDA  onto  glass  substrates, 
pre-purified  powder  samples  of  the  materials  were  loaded  into  separate  suL'imation 
cells  in  the  OMBD  growth  chamber  which  had  a  base  pressure  of  lO'"*  torr.  The 
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focused  beam,  this  birefringence  measurement  determines  the  local  orientation 
of  the  optic  axis  of  PTCDA  (NTCDA  has  only  very  weak  birefringence  and  hence 
the  probe  is  not  sensitive  to  its  presence  in  the  stacks).  This  technique,  therefore,  is 
sensitive  to  the  degree  of  crystalline  order  over  the  region  covered  by  the  beam 
diameter.  The  finding  that  the  optic  ax-'  of  PTCDA  is  orientated  in  the  same 
direction  across  the  wafer  surface,  and  that  the  birefringence  of  a  multilayer  stack 
is  equal  in  magnitude  to  a  single  PTCDA  film  of  equal  total  thickness,  indicates  that 
the  crystalline  orientation  of  one  PTCDA  layer  is  the  same  as  all  other  PTCDA  layers 
in  a  multilayer  structure.  By  inference,  the  NTCDA  layers  must  be  likewise  aligned. 
This  ordering,  both  across  wafer  surfaces  as  well  as  between  layers  in  the  stacks,  is 
evidence  for  quasi-epitaxy  in  the  multilayer  samples. 

The  process  of  q-e  can  be  quantitatively  understood  if  we  require  that  a  monolayer 
of  one  material  must  find  a  minimum  energy  configuration  with  the  underlying 
crystal  during  growth.  That  is,  when  growth  of  material  B  is  initiated  on  the  surface 
of  material  A,  the  first  few  molecules  B  align  themselves  to  A  to  minimize  their 
individual  molecular  interaction  energies.  As  the  growth  of  layer  B  approaches  a 
full  monolayer,  different  islands  of  material  B  arising  from  nucleation  at  random 
sites  across  the  surface  of  A  will  be  out  of  register  with  each  other,  although  they 
may  all  be  approximately  aligned.  As  the  spacing  between  islands  is  filled  in  with 
additional  deposition  of  material  B,  the  islands  must  shift  slightly  (but  in  a  rigid 
manner)  to  accommodate  these  last,  interstitial  molecules.  This  process  must  occur 
without  increasing  the  overall  energy  of  the  layer  to  the  extent  that  a  dislocation 
is  generated.  Note  that  this  picture  requires  that  the  substrate  temperature  be  high 
enough  to  allow  for  some  surface  mobility  of  the  islands  and  molecules,  yet  not  so 
high  that  the  individual  molecules  are  only  very  loosely  bound  to  the  surface.  For 
this  reason,  deposition  of  q-e  grown  materials  occurs  on  substrates  cooled  below 
room  temperature. 

One  condition  which  must  be  fulfilled  for  the  above  scenario  to  occur  is  that  the 
binding  energy  minimum  between  molecules  A  and  B  must  be  very  broad  -  i.e.  there 
is  a  range  of  positions  which  molecule  B  can  have  when  aligned  to  molecule  A  which 
does  not  significantly  change  the  total  energy  of  the  system.  In  this  way,  small  shifts 
in  position  can  be  accommodated  without  significantly  increasing  the  total  crystal 
energy.  We  have  tested  this  hypothesis  by  calculating  the  vdW  bond  energy  between 
PTCDA  (molecule  A)  and  NTCDA  (molecule  B).  These  calculations  are  done  using 
the  atom-atom  potential  method.'*  In  this  method,  we  calculate  (and  minimize) 
the  vdW  energy  between  each  pair  of  atoms  in  both  the  PTCDA  and  NTCDA 
molecules.  The  total  bond  energy  is  simply: 

#  =  (1) 

where  4>,j  is  the  atom-atom  potential  between  the  /th  and  y  th  atoms  in  molecules 
A  and  B,  respectively.  This,  in  turn,  is  calculated  using  the  Buckingham  potential; 

<i>o  =  «/'•«  +  l8cxp{-7r„)  (2) 

Here,  r,j  is  the  distance  between  atoms  /  in  molecule  A  and  j  in  molecule  B,  and  a, 
/3  and  7  are  vdW  atom-atom  potential  constants  for  the  constituent  atoms  in  the  two 
molecules  (C,  O  and  H  for  PTCDA  and  NTCDA,  for  example).  The  constants  used 
in  the  case  of  PTCDA  and  NTCDA  are  given  in  Table  1. 

The  accuracy  of  this  method  has  been  tested  by  calculating  the  minimum  energy 
configuration  of  bulk  PTCDA,  and  comparing  the  results  to  the  measured  crystal 
structure.  Under  this  test,  the  model  works  surprisingly  well.  For  example,  the 
calculated  intermolecular  stacking  distance  is  3.26  A  as  compared  with  the  actual 
value  of  3.21  A.  Furthermore,  the  model  predicts  that  two  molecules  stack  with 
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From  the  energy  contour  plot,  we  see  that  the  energy  minimum  between  PTCDA 
and  NTCDA  is  indeed  broad,  allowing  for  nearly  a  2  A  translation  of  one  molecule 
with  respect  to  another  without  any  significant  increase  in  #.  Thus,  we  can  conclude 
from  this  first  order,  somewhat  simplified,  calculation  that  our  ‘picture’  of  the 
process  of  q-e  is  plausible.  That  is,  we  expect  that  once  the  islands  are  nucleated 
around  the  surface  of  the  substrate,  they  are  still  free  to  translate  rigidly  to 
accommodate  the  arrival  of  the  final  molecules  which  fill  in  the  interstices  between 
neighbouring  islands.  This  shift  occurs  without  inducing  strain  or  polycrystalline 
growth. 

While  this  model  is  still  unsophisticated  in  many  ways,  it  suggests  the  properties 
which  molecules  must  possess  before  they  can  be  grown  into  q-e  structures.  From 
these  initial  results,  we  infer  that  the  process  is,  in  fact,  a  very  general  property  of 
vdW  solid  growth. 

One  consequence  of  being  able  to  reproduce  the  anisotropic  crystal  structure  in 
a  thin  film  is  that  the  films  exhibit  enormous  anisotropies  in  both  their  conductive 
and  dielectric  properties.  For  example,  the  in-plane  conductivity  of  PTCDA  films  is 
found  to  be  at  least  six  orders  of  magnitude  lower  than  the  conductivity  perpendicular 
to  the  film  plane. It  is  also  expected  that  the  dielectric  properties  of  highly  ordered 
films  would  have  a  considerable  degree  of  anisotropy  along  different  crystalline 
directions.  The  off-resonance  dielectric  constant  of  a  material  at  frequency,  w,  is 
given  by: 


(3) 


where  N  is  the  molecular  number  density,  Uj  is  the  dipole  transition  frequency,  q  is 
the  electronic  charge,  and  m*  is  the  electron  effective  mass.  The  oscillator  strength 
of  the  dipole  transition  from  level  o  to  j  is  given  by  fj  =  2/n*/tw  |  |  VA’,  where 

is  the  electron  position  expectation  value,  and  h  is  Planck’s  constant  divided  by 
It.  Assuming  that  the  largest  contribution  to  the  dipole  moment  of  the  molecule 
is  due  to  the  r-orbitals,  and  that  the  electron  within  a  given  orbital  is  completely 
delocalized,  then  we  can  make  the  approximation  that  * 

<P-/L^.  For  the  case  of  PTCDA,  d  ^  i.l  k  is  the  extent  of  the  x-orbital  system 
perpendicular  to  the  molecular  plane,  and  L  is  its  extent  in  the  plane,  which  is 
approximately  equal  to  the  length  of  the  perylene  molecular  core  of  PTCDA  (i.e. 
L  =  6.9  A).  Hence,  provided  that  there  is  perfect  crystalline  alignment  throughout 
the  q-e  film,  we  can  expect  an  anisotropy  in  e  between  directions  perpendicular  and 
parallel  to  the  thin  film  plane  (and  hence  approximately  perpendicular  and  parallel 
to  the  molecular  stacking  axis)  of  (e^  -  l)/(e,|  -  1)  =  d^/L‘  »  0.22.  This  value 
is  only  an  approximation,  since  we  are  assuming  that  m*  is  isotropic, “  and  that  the 
dipole  moment  is  due  to  a  completely  delocallized  electron  in  the  extended  x-system. 
Nevertheless,  it  has  been  found’  that  the  measured  value  of  -  O/le,  -  1)  for 
PTCDA  is  0.26,  which  is  quite  close  to  the  value  ‘predicted’  by  the  above  model. 
Furthermore,  the  index  of  refraction  measured  at  a  wavelength  of  A  =  1.064  ;xm 
in  the  direction  perpendicular  to  the  substrate  plane  is  =  1.36  ±  0.01,  whereas 
parallel  to  the  plane,  n,  =  2.017  ±  0.005,  resulting  in  an  index  difference  of 
An  =  0.66.  To  our  knowledge,  these  are  the  largest  index  anisotropies  ever  measured 
for  thin  films,  giving  further  evidence  for  the  nearly  perfect  crystalline  order  achieved 
in  q-e  growth. 


ORGANIC  HETEROJUNCTIONS 

It  is  well  known  that  multiple  quantum  wells  can  be  fabricated  by  alternately  layering 
two  ino'ganic  semiconductors  with  different  bandgaps.  Thus,  the  question  arises: 


QUASI-EPITAXIALLY  GROWN  CRYSTALLINE  ORGANIC  SEMICONDUCTORS  213 


the  hole  mobility  was  found  to  increase  linearly  at  low  electric  fields,  and  saturates 
at  high  fields.  The  high  field  carrier  velocity  saturation  is  due  to  phonon  emission, 
consistent  with  band  theory.  Further,  the  long  carrier  mean  free  path  (as  compared 
with  the  intermolccular  stacking  distance)  in  these  crystalline  organic  materials 
estimated  from  the  low  temperature  mobility  data  also  supports  the  band  description. 
Provided  that  we  can  use  this  model  to  describe  the  charge  transport  in  PTCDA  and 
CuPc,  the  transport  theory  for  inorganic  semiconductor  heteroj unctions  can  also 
be  applied  (with  caution)  to  organic  heteroj  unctions. 

Using  the  treatment  of  Chang  for  isotype  heterojunctions, the  saturation 
current  density  for  an  HJ  is  given  by: 

y,  =  exp [  -  9 Voi /kT]  (4) 

where  B  is  a  pre-factor  which  depends  on  the  density  of  acceptors  and  the  hole 
effective  mass  in  CuPc,  kT  is  the  Boltzmann  energy  at  temperature  T,  Kq,  is  the 
diffusion  potential  on  the  depleted,  CuPc  side  of  the  heterojunction,  and  q  is  the 
electron  charge.  The  valence  band  discontinuity  energy,  A£,,  is  then  given  by: 


^.  =  QyDi  +  Q^D2  (5) 

In  the  above  equation,  Fq,  is  the  diffusion  potential  on  the  PTCDA  side  of  the 
heterojunction.  For  simplicity,  we  assume  that  there  is  no  significant  difference  in 
the  effective  hole  densities  of  states  and  in  the  acceptor  concentrations  between  CuPc 
and  PTCDA. 

From  Equation  4,  it  can  be  seen  that  J,  is  thermally  activated  with  energy  qV^z- 
In  Figure  5,  we  show  the  dependence  of  log  (J,)  on  l/Tfor  a  CuPc/PTCDA  hetero¬ 
junction.  Here,  J,  is  obtained  by  extrapolating  both  the  forward  and  reverse  bias 
currents  to  an  applied  voltage  of  F  =  0.  The  activation  energies  for  both  the  for¬ 
ward  and  reverse  bias  are  the  same,  yielding  =  0.51  ±  0.05  V.  This  value  is 
consistent  with  that  obtained  from  capacitance-voltage  (C-V)  measurements,  where 
we  find  that  Vp:  =  0.45  ±  0.15  V.  Following  the  treatment  of  Chang,  we  obtain 
Fdi  =  0.01  V,  and  hence  the  valence-band  discontinuity  energy  for  CuPc/PTCDA 
heterojunctions  is  A£.  =  F^,  -I-  F^j  =  0.52  ±  0.05  V. 

Based  on  the  electrical  characteristics  of  the  CuPc/PTCDA  heterojunctions,  a 


FIGURE  5  Saturation  current  versus  temperature  for  a  CuPc/PTCDA  heterojunction. 


QUASI-EPITAXIALLY  GROWN  CRYSTALLINE  ORGANIC  SEMICONDUCTORS  2P 


to  quantum  confinement  of  excitons  in  a  quantum  well  structure.  We  now  discuss 
both  models  in  detail. 

The  change  of  polarization  energy  of  PTCDA  due  to  the  proximity  of  NTCDA 
molecules  in  closely  spaced  layers  can  be  calculated  using  perturbation  theory.  A 
derivation  of  the  perturbation  of  the  exciton  energy  due  to  polarization  effects  is 
given  in  the  Appendix.  The  results  of  the  calculation  show  that  the  first  order  energy 
shift  of  the  polarization  energy  increases  linearly  with  the  layer  thickness,  i.  This 
effect,  while  being  very  short  range,  can  actually  increase  the  exciton  energy  with 
increasing  well  width,  contrary  to  our  observations.  Hence,  we  can  rule  out  this  first 
order  ‘energy  offset’  effect.  On  the  other  hand,  the  second  order  polarization  energy 
correction  is  proportional  to  t~*.  As  shown  in  Figure  10,  the  data  are  not  consistent 
with  the  results  of  the  calculation  since  their  energy  dependence  is  considerably 
smaller  than  even  /"^  Hence,  we  conclude  that  the  polarization  model  does  not 
adequately  explain  the  observed  blue  shift  in  the  exciton  line. 

Quantum  confinement  provides  an  alternative  explanation  for  the  observed  blue 
shift  in  the  absorption  spectrum.'-^  Here,  we  use  the  variational  method  to  evaluate 
the  binding  energy  of  an  exciton  in  a  quantum  well.  The  Schrodinger  equation  for 
an  exciton  can  be  solved  using  the  trial  wavefunction; 

=  (6) 

where  ^i,(p,  z)  is  the  hydrogenic  Is  wavefunction,  p  =  +  y‘,  and  are 
coordinates  in  the  z-direction  for  the  electron  and  hole,  respectively,  and 

t^h(Zh)  arc  the  exact  solutions  to  the  finite  square  well  problem,  and  z  =  z,  -  Zh  -  In 
Equation  6,  the  hydrogenic  state  has  the  form: 

1^,.  =  exp(-(pVa^  +  zV|3^)''^]  (7) 

Here,  a  is  a  variational  parameter  in  the  trial  solution,  and  0  is  set  equal  to  a^, 
where  Oq  >s  the  exciton  Bohr  radius  in  the  bulk  crystal.  For  crystalline  organic 
MQW  structures,  typically  aa<  L^,  where  Z..  is  the  quantum  well  width.  In  this 
case,  therefore,  it  is  reasonable  to  assume  0  -  a^,  with  a  being  the  only  variation 
parameter. 

The  parameters  used  in  this  model  to  fit  the  data  in  Figure  10  are  listed  in  Table 
2.  Here,  the  sum  of  V,  and  is  equal  to  the  difference  in  energy  gaps  of  PTCDA 
and  NTCDA.  Of  all  the  parameters  used  in  the  calculation,  the  results  are  most 
sensitive  to  the  choice  of  rn„  and  the  ratio,  V^/Vf,.  Using  these  parameters,  a  good 
fit  to  the  experimental  data  is  obtained  for  =  O.lSmo,  where  is  the  electron 
rest  mass.  A  small  value  of  and  a  relatively  large  value  of  m,  were  chosen  in 
the  calculation  since  PTCDA  is  predominantly  a  hole-transporting  material.  The 
effect  of  using  different  values  of  are  also  shown  for  comparison  in  Figure  10 
indicating  the  sensitivity  of  the  fit  to  this  parameter.  Note  that  our  calculation 
implies  that  K,  ►  This  is  consistent  with  electrical  measurements  made  for 
PTCDA/NTCDA  heteroj unctions  which  do  not  exhibit  significant  rectification  of 
the  hole  current  at  temperatures  ranging  from  room  temperature  to  about  90  K.  This 


TABLE  2 

Variational  calculation  parameters 


Parameter 

Symbol 

Unit 

Value 

Hole  mass 

"•h 

">0 

0.18 

Electron  mass 

">0 

10.0 

Electron  potential 

y. 

meV 

900 

Hole  potential 

Relative  dielectric 

meV 

50 

constant  of  PTCDA 

3.6 
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exhibit  electrical  and  optical  properties  analogous  in  many  ways  to  fully  inorganic 
heterojunctions. 

5.  Evidence  for  quantum  confinement  in  very  thin  layer  organic  multiple  quantum 
well  structures  has  furthered  our  understanding  of  the  nature  of  excitons  in 
crystalline  organic  materials.  Given  that  excitons  in  these  materials  have  been  the 
focus  of  study  for  over  40  years, in  many  ways  they  are  the  key  to  our 
understanding  the  physics  of  vdW  solids  in  general.  In  addition,  the  observation 
of  quantum  confinement  suggests  that  many  devices  based  on  organic  MQWs 
are  becoming  a  realistic  possibility.  One  example  of  such  a  device  would  be 
an  organic  MQW  optical  modulator  which  could  be  made  to  cover  all  regions  oi 
the  ultra-violet,  visible  or  near  infra-red  spectral  regions,  depending  on  the  choice 
of  organic  compounds  employed. 

6.  The  large  exciton-phonon  coupling  strength  characteristic  of  organic  semicon¬ 
ductors  can  lead  to  exciting,  very  large  nonlinear  optical  effects.  For  example, 
Lam  and  co-workers^^  predicted  that  organic  MQWs  have  large  effects,  and 
are  therefore  potentially  useful  for  a  broad  range  of  applications  which  are 
unattainable  using  the  polymers  and  photorefractive  materials  studied  today. 
They  have  shown  that  quasi-epitaxially  grown  crystalline  organic  thin  films 
and  MQWs  with  large  Franck-Condon  (FC)  shift  energies  exhibit  an  intensity- 
dependent  absorption  and  index  of  refraction  at  relatively  low  input  powers. 
Here,  the  FC  shift  is  the  energy  difference  between  the  lowest  energy  absorption 
and  the  highest  energy  luminescence  spectral  peaks,  and  is  proportional  to  the 
square  of  the  exciton-phonon  coupling  strength.  As  shown  by  the  arrows  in 
Figure  10,  the  PTCDA/NTCD.A  MQW  structure  with  40  A  layer  thickness  has 
FC  =  4050  cm"',  which  to  our  knowledge  is  the  largest  yet  observed  for  organic 
semiconductors. 

To  understand  the  nonlinear  optical  properties  of  organic  MQWs,  the  Hamil¬ 
tonian  for  organic  MQWs  under  illumination  can  be  written^*: 

3C  =  -  it)i)a*a  —  hkQa*a  -t-  hugb*b  -  fia* E  -  fiaE*  (8) 

where  Wj,,  w,,  and  are  the  exciton,  incident  light  and  phonon  frequencies,  respec¬ 
tively,  and  A  is  the  exciton-phonon  coupling  constant,  n  =  qa„  is  the  cxciton  dipole 
moment,  E  is  the  external  electrical  field,  and  a*a  and  b*b  are  the  exciton  and 
phonon  populations,  respectively.  Also,  Q  =  b*  -b  6  is  the  phonon  amplitude. 

The  origin  of  the  optical  nonlinearities  in  organic  MQWs  lies  in  the  second 
term  of  Equation  8,  which  is  proportional  to  the  exciton-phonon  coupling  strength, 
A,  given  by: 

A=  (9) 

The  second  term  in  Equation  8  ‘renormalizes’  the  exciton  energy,  by  —hXQ.  As 
the  incident  optical  field  intensity  increases,  the  exciton  population  increases,  result¬ 
ing  in  a  decrease  in  the  energy  of  the  exciton  absorption  peak.  Such  an  effect  leads 
to  an  intensity  dependent  absorption,  and  hence  an  intensity  dependent  index  of 
refraction  through  the  Kramers-Kroenig  effect.  Note  that  if  the  incident  light  beam 
is  detuned  from  the  exciton  peak  on  the  low  side  of  by  a  factor  of  w,  -  ,  then 

the  absorption  peak  energy  is  shifted  toward  the  pump  light  energy,  which  decreases 
the  detuning.  This  increases  the  absorption  with  increasing  incident  light  intensity, 
leading  to  optical  bistability. 

The  two-wave  mixing  gain  for  a  weak  probe  beam  interacting  with  a  strong 
pump  beam  for  PTCDA/NTCDA  MQW  structures,  predicted  by  the  time  depen¬ 
dent  solution  of  Schrodinger’s  Equation  and  using  Equation  8,  is  shown  in  Figure  1 1 . 
The  horizontal  axis  corresponds  to  the  (Rabi  frequency)^  which  is  proportional  to 
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APPENDIX:  POLARIZATION  EFFECTS  IN  ORGANIC  MQWs 

Consider  an  exciton  with  a  Bohr  radius  a„  generated  at  point  P  in  a  PTCDA  layer 
of  thickness  t  sandwiched  between  two  NTCDA  layers  with  the  same  thickness,  as 
shown  in  Figure  Al.  The  change  in  exciton  potential  due  to  the  presence  of  an 
NTCDA  molecule  at  point  N  on  the  NTCDA  surface  is: 

Ay)  =  AiiCOsS'Ap  ji?  -  r|^  (A-1) 

where  =  pp  —  mn.  with  Mp(Mn)  the  dipole  moment  per  unit  volume  of  PTCDA 
(NTCDA),  tp  is  the  dielectric  constant  of  PTCDA,  i?  is  a  vector  between  points  P 
and  N,  d '  is  the  angle  between  R  and  the  normal  to  the  NTCDA  plane,  and  the  charge 
density  of  the  exciton  is  described  by  the  spherical  coordinates:  r,  6,  and  <t>.  Using 
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Current  research  in  the  optical  properties  of  organic 
materials  has  been  directed  toward  the  elucidation  of  the 
dominant  mechanism  that  gives  rise  to  their  nonlinear 
optical  behavior.'  In  spite  of  the  many  experiments  per¬ 
formed  over  the  past  decade,*  little  depth  of  understand¬ 
ing  has  been  achieved.  The  complications  arise  from  the 
competing  effects  between  the  delocalization  of  the  pho- 
togencrated  charge  carriers,*  and  formation  of  the  exci¬ 
tations  or  quasiparticles.'*  Recent  nonlinear  optical- 
absorption  experiments*  performed  on  the  quasi-lD 
semiconductor  polydiacetylene  appeared  to  confirm  the 
concept*  that  excitons  are  responsible  for  the  nonlinear 
optical  behavior  of  1 D  organic  materials. 

In  this  work,  we  have  extended  the  studies  of  semicon¬ 
ductor  multiple  quantum  wells  (MQWs)  to  the  case  of 
crystalline  organic  MQWs  (CO-MQWs)  (Ref.  7)  and 
have  found  that  the  response  of  the  quantum-confined 
charge-transfer  (CT)  excitons  to  external  fields  have 
novel  nonlinear  optical  properties.  CT  excitons  are 
known  to  exist  in  molecular  crystals*  and  their  electronic 
structure  can  be  described  by  means  of  the  Wannier  pic¬ 
ture*  with  an  appropriate  static  dielectric  constant. 
However,  their  interaction  with  the  lattice  is  significantly 
different  from  the  Wannier  excitons  found  in  inorganic 
semiconducting  hosts.  The  CT  exciton  binding  energy 
lies  in  the  few-eV  range,  making  them  less  susceptible  to 
phonon-induced  ionization  as  compared  to  Wannier  exci¬ 
tons  in  inorganic  semiconductors. 

In  previous  work,  the  linear  optical  properties  of  CO- 
MQWs  were  measured  in  some  detail.*  It  was  found 
that  the  lowest-encrgy  CT  exciton  absorption  line  was 
blueshifted  with  decreasing  well  width.  This  observation 
is  consistent  with  quantum  confinement  of  the  CT  exci¬ 
ton  by  energy  wells  formed  in  one  of  the  two  MQW  lay¬ 
ers  (consisting  of  3,4,9,  lO-perylenetetracarboxylic  di¬ 
anhydride  or  PTCDA),  sandwiched  between  layers  of  a 
second  material  (3,4,7,8-naphthalene  tetracarboxylic  di¬ 
anhydride  or  NTCDA)  forming  energy  barriers.  A  vari¬ 
ational  study  of  the  dependence  of  exciton  energy  on  well 
width  indicates  that  the  exciton  radius  is  approximately 
15  A.  This  number  is  significantly  smaller  than  that 


found  in  Ill-V  semiconducting  compounds  because  of  the 
smaller  static  dielectric  constant  of  organic  crystals.  The 
strong  Coulomb  binding  energy  (as  compared  with  the 
kinetic  energy  of  each  charge  carrier)  in  the  CO-MQW 
materials  implies  that  the  quantization  due  to  the  well 
thickness  L  is  determined  by  the  center-of-mass  motion 
of  the  exciton.  Hence  the  exciton  binding  energy  is  given 
by  hcDjc  “hwB  +  {hJtn)^llML^  for  an  infinitely  deep  po¬ 
tential  well,  where  q>b  is  the  bulk  exciton  binding  energy, 
n  is  an  integer,  and  M  is  the  total  mass  of  the  exciton. 
This  expression  provides  a  qualitative  explanation  of  the 
observed  blueshift  that  was  reported  in  the  linear  absorp¬ 
tion  measurements.  Since  the  CT  exciton  represents  a 
correlated  electron-hole  pair  between  nearly  adjacent 
molecules  in  a  stack,  we  can  consider  this  radius  to  be 
the  spatial  dimension  in  an  electric  dipole  moment.  Such 
a  large  dipole  moment  should,  in  turn,  lead  to  large  opti¬ 
cal  noniinearities  in  these  materials.  These  optical  non- 
linearities  are  the  focus  of  this  study. 

The  starting  point  of  our  analysis  of  the  optical  prop¬ 
erties  of  CT  excitons  in  crystalline  organic  materials  is 
the  Davydov  Hamiltonian  wi.ich  describes  the  interac¬ 
tion  of  excitons  with  phonons  and  external  radiation 
fields.  Assuming  the  rotating-wave  approximation,  ex¬ 
istence  of  one  phonon  mode,  and  keeping  only  the  linear 
exciton-phonon  interaction,  the  Hamiltonian  is 

H  “h{u^  —  Q))a*a  +  hQ)ob*b 

~ bka*aQ  —  j pa^E  —  j paE*  ,  (1) 

where  hcoj  and  Acuo  are  the  quantum-confined  exciton 
binding  and  phonon  energies,  respectively.  Also,  X  is  the 
exciton-phonon  coupling  constant,  p  is  the  electric  dipole 
moment  of  the  exciton,  a* a  and  b^b  are  the  exciton  and 
phonon  populations,  respectively.  Further,  Q^b  +  b*  is 
the  phonon  amplitude,  E  is  the  slowly  varying  envelope 
of  the  external  field,  and  a  and  b  are  the  exciton  coher¬ 
ence  and  phonon  annihilation  operators,  respectively. 

Careful  interpretation  of  the  constants  X  and  p  must 
be  considered.  Current  measurements  appear  to  be  in¬ 
conclusive  concerning  the  effects  of  quantization  on  these 
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constants.  From  a  theoretical  point  of  view,  quantum 
confinement  should  play  a  role  since  the  envelope  wave 
functions  are  quantized  and  they  enter  in  the  computa¬ 
tion  of  the  matrix  elements  of  the  observables.  That  is, 
the  quantum  confinement  modifies  the  bulk  values  by  a 
factor  that  reflects  the  overlap  of  the  spatial  wave  func¬ 
tions  for  the  electrons  and  holes. 

The  temporal  evolution  of  the  exciton  coherence  a  and 
phonon  amplitude  Q  are  determined  by  the  Heisenberg 
equation  of  motion,  and  are  given  by 

■^  + [/(<u»  —  <u)  + yla  .  (2a) 

at  2n 

— — ^ -l-r~^  +  (cuq) ™2aJoXfl  a,  (2b) 

dr  at 

where  y  and  T  are  phenomenological  cxciton  dephasing 
and  phonon  decay  rates,  respectively.  Equations  (2)  pro¬ 
vide  insight  into  the  nonlinear  optical  behavior  of  exci- 
ton-phonon-couplcd  systems.  The  term  ikQa  in  Eq.  (2a) 
is  a  renormalization  of  the  exciton  frequency  due  to  its 
coupling  to  the  phonon  structure  of  the  material.  Since 
it  depends  on  the  phonon  amplitude  Q,  the  renormaliza¬ 
tion  factor  can  be  seen  to  be  proportional  to  the  exciton 
density  from  the  steady-state  solution  of  Eq.  (2b).  This 
implies  that  the  effective  exciton  frequency  is  a  function 
of  the  population  of  photogenerated  excitons,  which  is 
proportional  to  the  light  intensity.  Hence,  optical  non- 
linearities  in  these  materials  have  their  origin  in  an 
exciton-phonon-induced  frequency  shift, ^  in  a  mar  er 
similar  to  the  dynamic  Stark  shift  in  polaritons. 

The  nonlinear  evolution  of  coupled  waves  is  deter¬ 
mined  by  the  Maxwell  equations.  In  the  slowly  varying 
envelope  approximation,  they  are  given  as 

(3a) 

cz  c 

k  / 

where  the  nonlinear  polarization  density  P,  is  defined  ’)y 
Pa“Nfiia).  (3b) 

Here,  N  is  the  number  density  of  CT  excitons  and  (o'  is 
the  expectation  value  of  a.  TTie  subscript  o  denotes  the 
radiation  field  oscillating  at  frequency  <0,. 

Equations  (l)-(3)  have  exact,  closed-form  analytical 
solutions  in  the  steady-state  regime.  We  shall  consider 
three  important  cases.  The  first  involves  the  linear 
response  of  the  medium  to  an  optical  radiation  field.  A 
comparison  of  the  theoretical  model  to  the  available  ex¬ 
perimental  data  will  provide  an  estimate  of  the  coupling 
parameters.  Second,  we  will  •  explore  the  nonlinear 
response  of  the  medium  by  obtaining  an  exact  solution  to 
the  CT  exciton  population  (c^a>.  Finally,  we  will  use 
the  results  to  understand  the  process  of  two-beam  cou¬ 
pling'  in  these  materials.  The  latter  involves  the  non¬ 
linear  coupling  of  strong  and  weak  radiation  fields. 

For  the  case  of  a  single  input  optical  wave,  the  analyti¬ 
cal  solution  of  Eq.  (2)  is  obtained  under  the  condition  of 


factorization  of  the  respective  variables,  a  and  Q.  First, 
in  the  low-intensity  regime,  the  population  of  photogen¬ 
erated  e.xcitons  is  proportional  to  the  intensity  of  the  op¬ 
tical  wave,  and  the  polarization  density  is  determined  by 
the  steady-state  small-signal  solution  of  Eqs.  (2).  That 
is. 


P 


(4) 


Figure  I  depicts  the  linear  absorption  (dark  solid  line) 
measurement  and  the  theoretical  fit  (light  solid  line) 
which  assumes  the  existence  of  two  excitonic  lines  in  the 
So  to  Si  electronic  manifold.*  The  theoretical  result  is 
fitted  against  the  experimental  data  in  the  following 
manner.  The  wavelengths  of  the  two  ma.xima  located  in 
the  right-hand  side  of  the  experimental  absorption  profile 
are  chosen  to  set  the  horizontal  scale.  The  vertical  scale 
is  set  by  the  ratio  of  the  two  peaks  in  the  right-hand  side 
of  the  data,  assuming  that  the  height  of  the  largest  peak 
is  equal  to  unity.  In  order  to  obtain  reasonably  good 
quantitative  agreement,  the  imaginary  part  of  the  polar¬ 
ization  density  was  averaged  over  a  Maxwellian  distribu¬ 
tion  with  two  distinct  widths,  which  are  chosen  from  the 
experimental  data.  The  use  of  the  Maxwellian  distribu¬ 
tion  is  consistent  with  the  fast  phonon-induced  relaxation 
processes  that  exist  in  these  materials.  The  result  shown 
in  Fig.  1  indicates  that  the  experimental  data  contain 
more  than  two  excitonic  lines,  which  is  the  source  of 
discrepancy  between  the  theoretical  model  and  experi¬ 
mental  data.  In  spite  of  the  complexity,  the  theoretical 
model  gives  an  adequate  understanding  of  the  origin  of 
the  linear  absorption  spectra  of  these  materials. 

Second,  in  the  fully  nonlinear  regime,  the  solution  for 
the  population  (a*o>  of  the  CT  exciton  is  given  by  the 


FIG.  1.  Linear  absorption  coefficient  of  PTCDA/NTCDA 
MQWs.  The  experimental  data  (dark  solid  line)  contain  two 
additional  sidebands  on  the  left-hand  side  due  to  the  presence 
of  the  NTCDA.  The  theoretical  fit  (light  solid  line)  from  the 
solution  of  the  exciton-phonon  equations  is  based  on  thermali- 
zation  by  phonons  and  the  assumption  of  two  exciton  reso¬ 
nances. 
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cubic  equation 

/  t  _  (S) 

where  A"(tfx~<u  is  the  detuning  from  the  cxciton  reso- 

•  nance  and  n  -p£/2A  is  the  Rabi  frequency.  Figure  2 
shows  the  solution  of  this  cubic  equation  as  a  function  of 
the  Rabi  frequency  for  different  values  of  the  detuning 
parameter  A.  In  this  plot,  all  physical  variables  have 
been  normalized  to  2A.Va>o-  A  transition  to  multivalue 
behavior  is  observed  for  a  sufficiently  large  value  of  the 

#  detuning  parameter  A.  This  behavior  can  be  understood 
in  the  following  manner.  Multivalue  behavior  of  Eq.  (5) 
is  achieved  if  the  derivative  of  the  Rabi  frequency  with 
respect  to  the  cxciton  population  changes  sign.  A  simple 
calculation  of  this  criterion  asserts  that  bistability  is 
present  provided  that 

*  A>>/3y.  (6) 


This  condition  is  valid  even  in  the  absence  of  an  optical 
cavity.  Hence,  the  coupled  exciton-phonon  system 
posseses  the  property  of  intrinsic  bistable  behavior  which 
arises  from  the  renormalization  of  the  exciton  frequency 
mediated  by  the  exciton-phonon  interaction. 

Finally,  we  consider  the  interaction  between  a  strong 
wave  £o.  oscillating  at  frequency  o»,  and  a  weak  wave 
£i,  oscillating  at  frequency  aj  +  5.  Their  nonlinear  cou¬ 
pling  yields  a  coherent  traveling-wave  excitation  in  the 
medium  oscillating  at  frequency  S.  The  scattering  of  the 
strong  wave  from  the  coherent  excitation  changes  the  ab¬ 
sorption  coefficient  and  the  index  of  refraction  experi¬ 
enced  by  the  weak  wave.  A  calculation  of  the  optical 


response  function  in  the  undepleicd  pump  approximation 
gives  the  following  expression  for  the  spatial  evolution  of 
the  weak  wave; 


I  dE\ 
£,  dz 


a)  +  5  /Vti^cuo  ,  . 

— - ria-Tf/j) 

2n,  eohlX^ 


(7a) 


where  the  dimensionless  (all  physical  parameters  are 
normalized  to  2A.Va»o)  nonlinear  absorption  coefficient  a 
is 


(a3i-S^)C-5rD 

C^+D^ 


(7b) 


and  the  dimensionless  nonlinear  index  of  refraction  is 
given  by 


P 


{ai-s^)D+src 

C^+D^ 


(7c) 


with  the  following  expressions  for  C  and  D: 
C  -  yCoxj  -  5-)  +  5r(A  -5)~  5r{a "a) , 


D  “(ci}i  —  S~)(A  —  S)  —Sr /  —  (Iciii  —  S^)(a*a) . 


Figure  3  shows  the  behavior  of  the  nonlinear  absorp¬ 
tion  coefficient,  a  +  i'P,  of  the  weak  wave  as  a  function  of 
the  Rabi  frequency  induced  by  the  strong  wave  for 
different  values  of  the  quantum-well  dimension.  A  tran¬ 
sition  to  bistable  behavior,  accompanied  by  gain  (nega¬ 
tive  values  of  the  nonlinear  absorption  coefficient)  of  the 
weak  wave  at  the  expense  of  the  strong  wave,  is  observed 
for  a  critical  value  of  the  normalized  Rabi  frequency  and 
a  small  enough  value  of  the  quantum-well  dimension. 
The  dimension  of  the  quantum  well  plays  a  key  role  in 


FIG.  2.  Dependence  of  the  exciton  population  on  the  nor¬ 
malized  Rabi  frequency  to  the  second  power.  The  linewidth 
and  the  detuning  are  also  normalized  to  2XV<uo,  which  has  the 
unit  of  frequency.  Curves  are  shown  for  normalized  detuning 
A  “3,  A  •■5.19,  which  corresponds  to  the  transition  region  for 
bistability  to  begin  taking  place,  and  A -8. 


FIG.  3.  Two-wave-mixing  gain  (negative  value)  and  absorp¬ 
tion  (positive  value)  coefficient  for  three  different  values  of  the 
quantum-well  size.  Prediction  of  optical  bistability  and  energy 
transfer  is  given  for  £,  “62  A.  The  normalized  probe-pump  de¬ 
tuning  parameter  S  is  set  equal  to  -  1. 
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Ihc  detuning  parameter  A.  For  large  enough  detuning  or 
small  enough  well  si^e,  the  value  of  the  detuning  param¬ 
eter  satisfies  the  bistability  condition  (6).  Hence,  a 
coherent  energy  transfer  from  the  strong  to  the  weak 
wave  takes  place  with  a  threshold  behavior.  This  phe¬ 
nomenon  can  be  thought  of  as  a  coherent  bistable  optical 
switch.  That  is.  the  energy  transfer  takes  place  from  the 
strong  to  the  weak  optical  beams  when  a  certain  thresh¬ 
old  is  achieved.  The  bistable  behavior  is  a  reflection  of 
the  nonlinear  functional  dependence  of  the  exciton  popu¬ 
lation  on  the  pump  intensity. 

The  phenomena  discussed  in  the  previous  paragraphs 
provide  insight  into  subtle  effects  that  appear  in  the 
Davidov  Hamiltonian.  Estimates  of  the  physical  param¬ 
eters  such  as  X  and  cuo  are  crucial  to  the  understanding 
of  the  materials  growth  conditions  as  well  as  to  the  fu¬ 
ture  applications  of  these  novel  materials  for  optoelec¬ 
tronics.  The  beauty  and  simplicity  of  our  results  are 
contained  in  one  single  physical  parameter,  the  Franck- 
Condon  (FC)  shift.  This  frequency  shift  is  related  to  X 
and  <uo  by  the  following  expression: 

FC-xV<uo. 

and  is  a  measure  of  the  degree  of  reduction  of  the  poten¬ 
tial  energy  of  the  material  due  to  the  cxciton-phonon 
coupling.  From  measured  values  in  aromatic  mole¬ 
cules,"  X  is  approximately  equal  to  wq.  These  numbers 
imply  a  FC  shift  of  700  cm  for  the  case  of  naphtalene 
compounds. 

It  is  interesting  to  calculate  the  range  of  laser  powers 
required  to  observe  the  onset  of  optical  bistability  and 
two-beam  energy  transfer.  However,  the  nonlinear  ab¬ 
sorption  coefficient  is  a  sensitive  function  of  the  material 
parameters.  For  example,  if  one  assumes  a  FC  shift  of 
approximately  1000  cm”‘  expected  for  large  molecules 
such  as  PTCDA,'^  an  oscillator  strength  of  unity  (con¬ 
sistent  with  the  observed  large  linear  absorption  coeffi¬ 
cients  of  10*  cm”'),  a  normalized  exciton  dephasing 
rate  ranging  from  10”*  to  1,  one  finds  that  the  power 
density  ranges  from  100  to  10*  W/cm*  for  A“I0”*FC. 


A  more  accurate  value  must  wait  for  a  detailed  measure¬ 
ment  of  the  FC  shift,  and  the  exciton  dephasing  rate. 

In  summary,  we  presented  new  results  on  the  optical 
behavior  of  these  newly  discovered  materials.  We  pre¬ 
dict  the  existence  of  intrinsic  bistability  as  well  a 
coherent  energy  transfer  from  a  strong  wave  to  a  weak 
wave.  The  energy  e.xchange  occurs  under  bistable  condi¬ 
tions,  leading  to  the  possibility  of  novel  optical  devices 
using  these  materials. 
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We  have  measured  the  indices  of  refraction  and  dielectric  constants  along  different  directions 
in  thin  films  of  the  crystalline  organic  semiconductor  compound  3,  4,  9,  10 
perylenetetracarboxylic  dianhydhde  (PTCDA).  The  films  were  deposited  via  organic 
molecular  beam  deposition,  resulting  in  single-crystalline,  “quasi-epitaxial"  films.  Due  to 
inherent  asymmetries  in  the  molecular  crystal  structure,  film  ordering  results  in  giant 
anisotropies  in  their  dielectric  properties.  For  example,  the  index  of  refraction  measured  at  a 
wavelength  of  A  =  1.064  pm  in  the  direction  perpendicular  to  the  substrate  plane  is 
=  1.36  ±  0.01,  whereas  parallel  to  the  plane,  =  2.017±0.(X)S,  resulting  in  an  index 
difference  of  An =0.66.  Furthermore,  the  low-frequency  dielectric  constant  of  the 
films  is  =  1.9  ±  0. 1  and  e||  =  4.5  ±  0.2.  To  our  knowledge,  these  are  the  largest  anisotropies 
ever  measured  for  thin  films.  We  discuss  a  guided  wave  polarization-selective  device 
which  takes  advantage  of  the  large  dielectric  anisotropies  characteristic  of  the  thin  organic 
films. 

Crystalline  organic  semiconductors  exhibit  a  diversity 
of  optical  and  electronic  properties  which  have  numerous 
device  applications.*'’  Recently,  we  have  demonstrated 
that  under  certain  deposition  conditions,  single-crystalline 
thin  films  of  the  organic  compound  3,  4,  9,  10  perylenetet- 
racarboyxlic  dianhydride  (PTCDA)  can  be  deposited  onto 
a  variety  of  substrate  materials  including  glass,  metal,  or 
semiconductors.  It  has  been  found  that  the  conductivity  of 
such  films  is  extremely  anisotropic*  due  to  jaymmetries  in 
the  crystalline  structure  of  the  organic  compounds  investi¬ 
gated.  For  example,  the  in-plane  conductivity  of  PTCDA 
films  is  found  to  be  at  least  six  orders  of  magnitude  lower 
than  the  conductivity  perpendicular  to  the  film  plane. 

One  would  expect  that  the  dielectric  properties  of  or¬ 
dered  films  would  also  be  anisotropic  along  different  crys¬ 
talline  directions.  A  unit  cell  of  PTCDA  is  shown  in  the 
inset  of  Fig.  1.  The  molecules  form  tilted  planar  stacks, 
with  a  distance  of  3.21  A  between  molecules  in  the  stack. 

The  long  axes  of  molecules  in  adjacent  stacks  are  rotated 
by  90°  to  the  long  axis  of  a  molecule  in  the  stack  at  the 
center  of  the  unit  cell.  Now,  the  off-resonance  dielectric 
constant  for  a  material  is  given  by 

e  —  l=4irNq^//[m*(cjj  —  0)^)],  (I) 

where  N  is  the  number  density  of  the  material,  <uy  is  the 
dipole  transition  frequency,  a  is  the  frequency  of  the  inci¬ 
dent  electromagnetic  radiation,  q  is  the  electronic  charge, 
and  m*  is  the  electron  effective  mass.  The  oscillator 
strength  of  the  dipole  transition  from  level  o  to  j  is  given  by 
fj  =  where  xqj  is  the  electron  position  ex¬ 

pectation  value,  and  is  Planck’s  constant  divided  by  2n. 

Assuming  the  largest  contribution  to  the  dipole  moment  of 
the  molecule  is  due  to  the  ir  orbitals,  and  that  the  electron 
within  a  given  orbital  is  completely  delocalized,  then  we 
can  make  the  approximation  //x///ii  ~1 
~  S/V-.  Here,  das 3.2  A  is  the  extent  of  the  rr  orbital 
perpendicular  to  the  molecular  plane,  and  L  is  its  extent  in 
the  plane,  which  is  taken  as  the  length  of  the  perylene 


molecular  core  of  PTCDA  (i.e.,  I.=s6.9  A).  Due  to  the 
weak  interaction  between  molecules  in  crystalline  organic 
semiconductors,  this  single  molecule  approximation  also 
can  be  accurately  applied  to  bulk  crystals.^  Hence,  pro¬ 
vided  that  there  is  perfect  crystalline  alignment  throughout 
the  thin  film,  we  can  expect  an  anisotropy  in  e  between 
directions  perpendicular  and  parallel  to  the  thin-film  plane 
of(Ci  —  l)/(eii  —  1)  ^  dVZ.^  ~  0.22.  This  value  is  only  an 
approximation,  since  we  are  assuming  that  m*  is  isotropic, 
and  that  the  dipole  moment  is  due  to  a  completely  delo¬ 
calized  electron  in  the  extended  ir  system.  Nevertheless,  it 
is  shown  below  that  this  value  is  close  to  the  measured 
anisotropy  in  f  for  quasi-epitaxial  films  of  PTCDA. 

To  measure  the  dielectric  constant  along  different  film 
directions,  capacitors  oriented  along  the  various  film  axes 
were  fabricated  on  quartz  substrates.  For  measuring  the 
capacitance  along  the  in-plane  direction,  interdigitated  pat¬ 
terns  of  Cr-Au  contacts  were  employed.  Here,  a  thin  layer 
(  —  150  A)  of  Cr  was  deposited,  followed  by  a  0.7-//m- 
thick  Au  layer  which  forms  a  current-blocking  contact  to 


FIG.  l.  Capacitance  measured  vs  frequency  for  a  horizontal,  interdigi¬ 
tated  capacitor  with  PTCDA  serving  as  the  dielectric.  Inset:  Two  views  of 
a  unit  cell  of  PTCDA. 
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PTCDA  ‘  Next,  the  contacts  were  patterned  into  2-nun- 
long  interdigitated  electrodes  with  2  /xm  spacing  and 
width.  Some  electrode  sets  were  positioned  perpoidicular 
to  others  to  determine  if  there  were  significant  asymmetries 
in  the  capacitance  in  the  two  in-plane  directions.  Next, 
purified  PTCDA  was  deposited  to  a  thickness  of  0.9  urn  via 
organic  molecular  beam  deposition  (OMBD)  onto  the 
quartz  substrate  which  was  maintained  at  a  temperature  of 
approximately  90  K  during  growth.  As  has  been  discussed 
previously,*  deposition  under  these  conditions  ensures  that 
the  resulting  film  is  single  crystalline,  independent  of  the 
substrate  material  employed.  The  capacitance  measured 
between  100  Hz  and  10  MHz  is  shown  in  Fig.  1.  The 
capacitance  is  independent  of  frequency  over  the  measure¬ 
ment  range,  indicating  that  the  films  are  free  of  traps. 

It  can  be  shown  that  the  ratio  of  capacitances  of  the 
interdigitated  contact  pattern  measured  before  (Cy^,)  and 
after  (C„)  deposition  of  the  PTCDA  thin  film  is  simply’ 
=  (c,  +  l)/(6,  +  e||),  where  e,  is  the  dielectric  con¬ 
stant  for  quartz.  We  note  that  the  small  amount  of  electric 
field  penetration  above  the  PTCDA  film  can  lead  to  devi¬ 
ations  from  this  expression.  However,  since  the  PTCDA 
thickness  was  larger  than  that  of  the  contact  huge's,  this 
error  is  expected  to  lead  to  an  insignificant  underestimate 
of  Cy.  By  measuring  and  C„  for  these  samples,  we 
obtain  in-plane  dielectric  constants  for  PTCDA  of  ey 
=  4.6  ±  0.2  and  4.3±0.2  for  the  two  perpendicular  orien¬ 
tations  of  the  patterns.  Since  these  values  are  within  each 
other's  error  limits,  we  assume  that  the  dielectric  constants 
along  the  two  directions  are  not  significantly  different, 
leading  to  fy  =  4.S  ^  0.2. 

To  measure  a  parallel-plate  capacitor  consisting  of  a 
PTCDA  film  (0.9  /xm  thick)  sandwiched  between  two 
Cr-Au  contacts  was  used.  Using  this  geometry,  it  is  found 
that  €i  =  1.9  ^  0.1  over  the  same  frequency  range  as  that 
used  to  measure  ey.  From  these  data  we  obtain  (ex 
—  1 )/(«!(  —  1 )  =  0.26  which  is  close  to  our  estimate  of 
0.22  for  p^ectly  oriented  crystalline  films.  In  earlier  mea¬ 
surements  of  ex  for  PTCDA  deposited  under  conditions 
leading  to  polycrystxdline  growth,  it  was  found  that*  ex 
=  3.6  ±  0.4.  Assuming  a  homogeneous  medium  was  being 
analyzed  in  that  case,  then  the  spatially  averaged  dielectric 
constant  is  e  ai  (2ey  -I-  ex)/3.  Using  the  values  of  e  ob¬ 
tained  here,  we  find  ew3.6,  in  agreement  with  the  former 
value.  This  indicates  that  the  films  deposited  on  room- 
temperature  substrates  are  structurally  randomized, 
whereas  deposition  via  OMBD  on  cold  substrates  results  in 
nearly  perfect  crystalline  order  across  broad  areas. 

Since  =  R^  where  n  is  the  refractive  index,  and  e.  is 
the  thin-film  dielectric  constant  at  optical  frequencies,  we 
also  expect  to  observe  asymmetries  in  n  as  a  result  of  the 
ordering  of  the  PTCDA  films.  Previously*  we  measured  an 
in-plane  refractive  index  of  Ry  =  2.017  ±  0.005.  In  this 
work,  therefore,  we  measure  only 

Determination  of  Rx  proceeds  by  measuring  the  reflec¬ 
tivity  from  the  thin  film  as  a  function  of  beam  incident 
angle  using  several  film  thicknesses  and  light  wavelengths. 
Considering  an  isotropic  homogeneous  thin  film  sand- 
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FIG  2.  Reflectance  at  2  a:  1.064  /im  vs  incident  beam  angle  for  a  I  .2S- 
ftm-diick  PTCDA  film  on  a  quaitz  substrate.  Data  are  shown  as  closed 
circles,  and  theory  is  shonm  as  a  solid  line.  Inset:  Index  of  refraction 
ellipsoid  for  PTCDA  showing  the  relationship  between  the  indices  along 
different  thin-film  directions. 

wiched  between  two  semi-infinite  media,  the  reflectivity  is 
given  by'° 

+  ^3  -t-  2r|2r„  ecK  2)8 
*  +  12^23 

where  p  =  2vn2t  cos  Bj/X.  Here,  r^  is  the  refractive  index 
measured  along  the  refractive  angle  62  in  the  thin  film  of 
thickness  t.  Also,  r,2  (r2})  is  the  reflectivity  at  the  interface 
between  medium  1  (3)  and  medium  2.  The  reflectivity 
maxima  and  minima  for  a  given  index,  R2,  occur  at  angles 
$2  which  satisfy  fl2/  =  mX/4  cos  $2,  where  m  is  the  integer 
order  of  the  extremum.  If  the  film  thickness  is  accurately 
determined,  the  refractive  index  can  thence  be  obtained  by 
measuring  the  angles  corresponding  to  reflectivity  extrema 
as  predicted  by  these  equations. 

To  extend  this  technique  to  anisotropic  thin  films,  the 
incident  light  polarization  is  made  parallel  to  the  plane  of 
incidence,  i.e.,  it  is  a  TM-polarized  wave.  Furthermore, 
using  birefringence  measurements,*  the  optical  b  axis  is 
located,  and  is  also  oriented  parallel  to  the  incident  light 
polarization  vector.  When  the  incident  angle  (0i)  of  the 
beam  to  the  film  is  changed,  this  routes  62  in  an  index 
ellipsoid  with  axes  of  length  Ry  =  r^  and  Rx,  as  shown  in  the 
inset  in  Fig.  2.  Here,  r»  is  the  index  along  the  b  axis.  From 
this  figure,  it  follows  that  the  index  rI  =  Ry  +  (1 
—  Ry/nJ)sin  ^  0|,  from  which  is  extracted. 

A  typical  measurement  result  made  at  A  =  1 .064  pm 
for  a  1.2S-;xm-thick  PTCDA  film  together  with  a  theoret¬ 
ical  calculation  using  Eq.  (2)  is  shown  in  Fig.  2.  The  ex¬ 
perimental  setup  consists  of  a  YAG  laser  providing  the 
polarized  incident  beam.  A  large  area  detector  is  posi¬ 
tioned  near  the  top  film  surface  to  measure  the  reflected 
light  intensity.  The  PTCDA  is  deposited  on  a  quartz  sub¬ 
strate  at  low  temperatures.  The  bwk  surface  of  the  quartz 
is  frosted  to  minimize  reflections  from  the  quartz/air  in¬ 
terface.  The  positions  of  the  theoretical  extrema  (particu¬ 
larly  the  minima)  can  be  made  to  match  those  obuined 
experimenully,  although  the  magnitude  of  the  reflected 
signal  is  sometimes  different  from  that  predicted  due  to 
unwanted  reflections  from  surface  imperfections  and  the 
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Table  I.  Reflectivity  data  for  PTCDA  Uun  films. 


Mistn) 

Sample 

thickness 

Incident 

angle* 

n,  (Cal). 

1.064 

1.13 

29* 

I.3S4  (m»8) 

1.319 

1.25 

4r 

1.364  Im^t) 

1.064 

1.13 

41* 

1.343  (in»6) 

1.319 

1.25 

55.5* 

1.346  (m»6) 

‘Coriesfiondiog  to  the  reflectivity  tninima. 


quartz/air  interface.  Similar  fits  are  obtained  for  a  film 
thickness  of  1.13  fua.  and  at  A  =  1.3  fim,  and  the  data  are 
listed  in  Table  I.  From  these  data,  we  obtain  =  1.36 
±  0.01  at  A  =  1.064 /xm,  giving  An=/i||  —  =  0.66,  which 
apparently  is  the  largest  value  of  Ait  reported  for  thin  films 
far  from  their  absorption  edge  (which  for  PTCDA  lies  at 
Aai6000  A‘).  Very  little  dispersion  is  observed  in  n  mea¬ 
sured  at  A  =s  1.064  and  1.3  fim,  as  indicated  in  Table  I. 

It  is  interesting  to  compare  the  asymmetries  measured 
for  the  dielectric  constant  and  the  index  of  refraction.  Ig¬ 
noring  dispersion  between  the  frequencies  at  which  these 
parameters  were  determined,  we  find  that  =  2.4 
»  ( nil/ni  =  2.2.  This  is  a  remarkable  agreement  given  the 
widely  different  conditions  under  which  the  various  mea¬ 
surements  were  made.  Indeed,  we  see  that  the  low-fre¬ 
quency  dielectric  constant  is  approximately  equal  to  for 
each  of  the  various  film  directions. 

A  waveguide  polarizer  was  fabricated  to  take  advan¬ 
tage  of  the  large  birefringence  characteristic  of  OMBD- 
grown  PTCDA  thin  films.  Thus,  PTCDA  rib  waveguides 
were  made  of  lengths  of  1.5  and  10  mm.  The  guides  were 
fabricated  by  first  spinning  a  l-/tm-thick  film  of  the  pho¬ 
toresist  AZ  1400  (with  «pR  =  1.61 )  onto  the  ( 100)  surface 
of  an  InP  wafer  (with  n,  =  3.27  at  A  =  1.064  /xm).  Next, 
2-/xm-wide  ridges  were  photolithographically  patterned  in 
the  AZ  1400,  followed  by  deposition  of  approximately  1 
/xm  of  PTCDA.  The  regions  where  PTCDA  was  deposited 
onto  the  photoresist  ridges  form  the  guides  to  TE  modes 
(where  ii||  >  Hpn),  but  not  to  TM  modes  (where 
<  Rpr).  The  substrates  were  cleaved  along  (110)  planes, 
allowing  for  coupling  into  the  end  facets.  Similar  guides 
have  been  measured  to  have  a  loss  of  <  2.5  dB/cm. 

Due  to  the  selective  guiding  of  TE  modes,  the 
waveguides  form  a  strongly  polarizing  medium.  To  mea¬ 
sure  the  TM/TE  extinction  ratio,  A  =  1.064  /xm  light  from 
a  diode-pumped  YAG  laser  was  coupled  into  and  out  of 
the  guides  using  microscope  objective  lenses.  The  input 
light  polarization  was  continuously  rotated  through  360* 
using  a  quarter  wave  plate  placed  between  the  laser  and  the 
sample.  In  addition,  a  Glan-Thomson  polarizer  (GTP) 
with  an  extinction  ratio  >  50  dB  was  placed  in  front  of  the 
laser.  This  setup,  shown  in  the  inset  of  Fig.  3,  ensures  a 
constant  coupling  efficiency  of  a  polarized  beam  along  any 
direction.  The  outcoupled  light  was  passed  through  a  sec¬ 
ond  GTP  which  analyzed  the  light  prior  to  being  detected 
using  a  Si  CCD  camera. 

Measurements  for  the  10-mm-long  device  are  shown  in 
Fig.  3.  Here,  both  the  TE  and  TM  components  of  the 
output  light  beam  are  plotted  versus  the  input  light  polar- 


FIG.  3.  Transmitted  light  intensity  for  TE  and  TM  waves  in  a  lO-mm- 
long,  2.^m-wide  PTCDA  rib  waveguide  as  a  function  of  incident  angle  of 
polarization.  Measurements  are  made  at  2  =  1.064  ^m.  Inset:  Experimen¬ 
tal  setup  used  to  characterize  the  PTCDA  waveguide  polarizer  under 
study. 

ization  angle.  An  extinction  ratio  of  TM/TE  <  —  48  dB 
was  obtained,  where  the  measurement  accuracy  was  lim¬ 
ited  by  the  GTP.  Similar  results  were  obtained  for  the  1.5 
mm  guide,  where  an  extinction  ratio  of  TM/TE<  —  25  dB 
was  measured.  The  accuracy  here  was  limited  by  light  scat¬ 
tered  over  the  top  of  the  short  guide.  The  propagation  loss 
of  the  short  guide  was  found  to  be  <0.4  dB. 

In  summary,  extremely  large  anisotropies  in  dielectric 
properties  have  been  observed  in  crystalline  organic  thin 
films  of  the  compound,  PTCDA.  The  anisotropies  in  both 
the  dielectric  constant  and  index  of  refraction  for  these 
films  is  a  result  of  their  nearly  perfect  crystalline  order 
achieved  via  OMBD  on  cold  substrates.  Similar  results  are 
expected  for  other  anisotropic  organic  semiconductor  crys¬ 
tals  deposited  using  the  OMBD  technique.  Furthermore,  a 
TE-pass  waveguide  polarizer  with  a  high  TM  polarization 
extinction  and  low  TE  propagation  loss  was  demonstrated. 
Since  such  polarizing  guides  can  be  deposited  on  semicon¬ 
ductor  substrates,  these  films  have  uses  such  as  optical 
isolators  integrated  with  semiconductor  lasers  ot  ampli¬ 
fiers. 
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